mat.h

Go to the documentation of this file.

#ifndef _MAT_H_
#define _MAT_H_

#include <mc/mc_base.h>

#include <mc/mtool.h>
#include <mc/slu.h>


struct sMat {

    char   name[10];  
    Vmem   *vmem;     
    int    iMadeVmem;  

    MATformat format;
 
    MATstate state;    

    MATsym sym;        

    MATimpl impl;   

    int    numR;   
    int    numC;  
    int    numA;       
    int    numO;      
    int    numZ; 
    int    numBR;    
    int    numBC;      

    int    iMallocIJA;
    int    iMallocA;   
    int    *IJA;       
    double *A;    
    int    *BR;   
    int    *BC; 

    int    *IA;  
    int    *JA;   
    double *diag;  
    double *offU; 
    double *offL; 

    Slu    *slu;
    Vset   *lnkL;   
    Vset   *lnkU;
    void   *xln;   
    void   *xlnt;   

};

typedef struct sMat Mat;

/*
 * ***************************************************************************
 * Class Mat: Inlineable methods (mat.c)
 * ***************************************************************************
 */

#if !defined(VINLINE_BAM)
#else /* if defined(VINLINE_BAM) */
#endif /* if !defined(VINLINE_BAM) */

Mat* Mat_ctor(Vmem *vmem, const char *name, int pnumR, int pnumC);

void Mat_dtor(Mat **thee);

void Mat_initStructure(Mat *thee,
    MATformat frmt, MATsym sym, int numO, int *IJA, double *A);

void Mat_copyStructure(Mat *thee, Mat *model);

void Mat_killStructure(Mat *thee);

int Mat_numR(Mat *thee);

int Mat_numC(Mat *thee);

int Mat_numA(Mat *thee);

int Mat_numO(Mat *thee);

int Mat_numZ(Mat *thee);

MATformat Mat_format(Mat *thee);

MATsym Mat_sym(Mat *thee);

MATstate Mat_state(Mat *thee);

MATimpl Mat_impl(Mat *thee);

void Mat_setFormat(Mat *thee, MATformat format);

void Mat_setSym(Mat *thee, MATsym sym);

void Mat_setState(Mat *thee, MATstate state);

void Mat_setImpl(Mat *thee, MATimpl impl);

int Mat_sizeIJA(Mat *thee);

int Mat_sizeA(Mat *thee);

int *Mat_IJA(Mat *thee);

int *Mat_IA(Mat *thee);

int *Mat_JA(Mat *thee);

double *Mat_A(Mat *thee);

double *Mat_diag(Mat *thee);

double *Mat_offU(Mat *thee);

double *Mat_offL(Mat *thee);

void Mat_print(Mat *thee);

void Mat_printSp(Mat *thee, char *fname, int pflag);

void Mat_printNoD(Mat *thee);

void Mat_printSpNoD(Mat *thee, char *fname, int pflag);

void Mat_zero(Mat *thee);

void Mat_set(Mat *thee, int i, int j, double val);

void Mat_addTo(Mat *thee, int i, int j, double val);

void Mat_buildBRC(Mat *thee, int numBR, int numBC, int *BR, int *BC);

void Mat_zeroBRC(Mat *thee);

void Mat_diagBRC(Mat *thee);

void Mat_galerkin(Mat *thee, Mat *rmat, Mat *amat, Mat *pmat);

int Mat_sluDirect(Mat *thee);

void Mat_sluCreate(Mat *thee);

int Mat_sluFactor(Mat *thee);

int Mat_sluSolve(Mat *thee, int key, double *f, double *u);

void Mat_sluDestroy(Mat *thee);

void Mat_memChk(Mat *thee);

void Mat_copy(Mat *thee, Mat *smat);

void Mat_copyROW2COL(Mat *thee, Mat *smat);

void Mat_copyCOL2ROW(Mat *thee, Mat *smat);

void Mat_copyDRC2RLN(Mat *thee, Mat *smat);

void Mat_copyROW2RLN(Mat *thee, Mat *smat);

void Mat_copyCOL2RLN(Mat *thee, Mat *smat);

void Mat_copyRLN2ROW(Mat *thee, Mat *smat);

void Mat_copyDRC2CLN(Mat *thee, Mat *smat);

void Mat_copyROW2CLN(Mat *thee, Mat *smat);

void Mat_copyCOL2CLN(Mat *thee, Mat *smat);

void Mat_copyCLN2COL(Mat *thee, Mat *smat);

void Mat_copyCLN2RLN(Mat *thee, Mat *smat);

void Mat_copyRLN2CLN(Mat *thee, Mat *smat);

void Mat_copyDRC2XLN(Mat *thee, Mat *smat);

void Mat_copyROW2XLN(Mat *thee, Mat *smat);

void Mat_copyCOL2XLN(Mat *thee, Mat *smat);

void Mat_copyRLN2XLN(Mat *thee, Mat *smat);

void Mat_copyCLN2XLN(Mat *thee, Mat *smat);

void Mat_copyXLN2DRC(Mat *thee, Mat *smat);

void Mat_copyXLN2ROW(Mat *thee, Mat *smat);

void Mat_copyXLN2COL(Mat *thee, Mat *smat);

void Mat_copyXLN2RLN(Mat *thee, Mat *smat);

void Mat_copyXLN2CLN(Mat *thee, Mat *smat);

void Mat_squeezeBRC(Mat *thee, int key);

void Mat_copy2(Mat *Y, Mat *X);

void Mat_axpy(Mat *Y, Mat *X, double val, int key);


void Mat_initStructureLN(Mat *thee, MATformat frmt, MATsym sym);

void Mat_killStructureLN(Mat *thee);

LinkRC* Mat_accessXLN(Mat *thee, int idx, int key);

void Mat_contribXLN(Mat *thee, int key, int i, int j, double val);

void Mat_contribNSYMXLN(Mat *thee, int key, int i, int j, double val);

void Mat_contribSSYMXLN(Mat *thee, int key, int i, int j, double val);

void Mat_contribSYMXLN(Mat *thee, int key, int i, int j, double val);

void Mat_printLN(Mat *thee);

void Mat_printLNSp(Mat *thee, char *fname, int pflag);

#endif /* _MAT_H_ */

---