scamac_generator.h

Go to the documentation of this file.

#ifndef SCAMAC_GENERATOR_H
#define SCAMAC_GENERATOR_H

#ifndef __cplusplus
#include <complex.h>
#endif


#include "scamac_error.h"
#include "scamac_inttypes.h"
#include "scamac_defs.h"
#include <stdbool.h>

typedef struct scamac_generator_st ScamacGenerator;
typedef struct scamac_workspace_st ScamacWorkspace;

ScamacErrorCode scamac_generator_check(const ScamacGenerator * gen, char ** desc);
ScamacErrorCode scamac_generator_finalize(ScamacGenerator * gen);
ScamacErrorCode scamac_generator_destroy(ScamacGenerator * gen);

ScamacErrorCode scamac_workspace_alloc(const ScamacGenerator * gen, ScamacWorkspace ** ws);
ScamacErrorCode scamac_workspace_free (ScamacWorkspace * ws);

ScamacErrorCode scamac_generate_row(const ScamacGenerator * gen, ScamacWorkspace * ws, ScamacIdx irow, ScamacFlag flag,
                                    ScamacIdx * nzr, ScamacIdx * cind, double * val);

/* more specialized versions */
ScamacErrorCode scamac_generate_row_real    (const ScamacGenerator * gen, ScamacWorkspace * ws, ScamacIdx irow, ScamacFlag flag, ScamacIdx * nzr, ScamacIdx * cind, double * val);
#ifndef __cplusplus
ScamacErrorCode scamac_generate_row_cplx    (const ScamacGenerator * gen, ScamacWorkspace * ws, ScamacIdx irow, ScamacFlag flag, ScamacIdx * nzr, ScamacIdx * cind, double complex * val);
#endif
ScamacErrorCode scamac_generate_row_int     (const ScamacGenerator * gen, ScamacWorkspace * ws, ScamacIdx irow, ScamacFlag flag, int * nzr, int * cind, double * val);
ScamacErrorCode scamac_generate_row_int_real(const ScamacGenerator * gen, ScamacWorkspace * ws, ScamacIdx irow, ScamacFlag flag, int * nzr, int * cind, double * val);
#ifndef __cplusplus
ScamacErrorCode scamac_generate_row_int_cplx(const ScamacGenerator * gen, ScamacWorkspace * ws, ScamacIdx irow, ScamacFlag flag, int * nzr, int * cind, double complex * val);
#endif


ScamacErrorCode scamac_alloc_cind_val(const ScamacGenerator * gen, ScamacFlag flag, ScamacIdx ** cind, double ** val);

const char * scamac_generator_query_name (const ScamacGenerator * gen);

ScamacIdx scamac_generator_query_nrow    (const ScamacGenerator * gen);
ScamacIdx scamac_generator_query_ncol    (const ScamacGenerator * gen);
ScamacIdx scamac_generator_query_maxnzrow(const ScamacGenerator * gen);
ScamacIdx scamac_generator_query_maxnzcol(const ScamacGenerator * gen);
ScamacIdx scamac_generator_query_maxnz   (const ScamacGenerator * gen);
ScamacIdx scamac_generator_query_valtype (const ScamacGenerator * gen);
ScamacIdx scamac_generator_query_symmetry(const ScamacGenerator * gen);

 int scamac_generator_query_coorddim(const ScamacGenerator * gen);
bool scamac_generator_get_coord     (const ScamacGenerator * gen, ScamacWorkspace * ws, ScamacIdx idx, double * x, double * y, double * z, int * ilayer, ScamacIdx * iblock);
//bool scamac_generator_get_coord     (const ScamacGenerator * gen, ScamacIdx idx, double * x, double * y, double * z, int * ilayer, ScamacIdx * iblock);
#endif /* SCAMAC_GENERATOR_H */