am.h

Go to the documentation of this file.

#ifndef _AM_H_
#define _AM_H_

#include <mc/mc_base.h>

#include <mc/aprx.h>

/*
 * ***************************************************************************
 * Class AM: Parameters and datatypes
 * ***************************************************************************
 */

struct sAM {

    /* ------- Objects out of our control (always exist) ------------------ */

    Vmem       *vmem;        
    int        iMadeVmem;   

    Aprx       *aprx;       
    Bmat       *P;       

    /* ------- Objects under our control ---------------------------------- */

    int        algExist;    

    Bmat       *A;   
    Bmat       *M;   

    Bvec       *f;  
    Bvec       *u;     
    Bvec       *ud;  
    Bvec       *ui;   
    Bvec       *ut;   
    Bvec       *r; 
    Bvec       *w0;     

};

typedef struct sAM AM;

/*
 * ***************************************************************************
 * Class AM: Inlineable methods (am.c)
 * ***************************************************************************
 */

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

/*
 * ***************************************************************************
 * Class AM: Non-inlineable methods (am.c)
 * ***************************************************************************
 */

AM* AM_ctor(Vmem *vmem, Aprx *taprx);

void AM_dtor(AM **thee);

void AM_create(AM *thee);

void AM_destroy(AM *thee);

int AM_markRefine(AM *thee, int key, int color,
    int bkey, double elevel);

int AM_refine(AM *thee, int rkey, int bkey, int pkey);

int AM_unRefine(AM *thee, int rkey, int pkey);

int AM_deform(AM *thee);

int AM_read(AM *thee, int key, Vio *sock);

double AM_assem(AM *thee,
    int evalKey, int energyKey, int residKey, int tangKey, int massKey,
    int bumpKey,
    Bvec *u, Bvec *ud, Bvec *f, int ip[], double rp[]);

double AM_evalJ(AM *thee);

void AM_evalFunc(AM *thee,
    int number, int block, int numPts, double *pts,
    double *vals, int *marks);

void AM_bndIntegral(AM *thee);

double AM_evalError(AM *thee, int pcolor, int key);

void AM_applyDiriZero(AM *thee, Bvec *which);

void AM_iniGuess(AM *thee, Bvec *which);

int AM_part(AM *thee, int pkey, int pwht, int ppow);

int AM_partSet(AM *thee, int pcolor);

int AM_partSmooth(AM *thee);

void AM_printJ(AM *thee);

void AM_printA(AM *thee);

void AM_printAnoD(AM *thee);

void AM_printAsp(AM *thee, char *fname);

void AM_printAspNoD(AM *thee, char *fname);

void AM_printP(AM *thee);

void AM_printPsp(AM *thee, char *fname);

void AM_printV(AM *thee, int num);

void AM_printVsp(AM *thee, int num, char *fname);

void AM_writeGEOM(AM *thee, Vio *sock,
    int defKey, int colKey, int chartKey, double gluVal, int fkey,
    int number, char *format);

void AM_writeSOL(AM *thee, Vio *sock, int number, char *format);

void AM_memChk(AM *thee);

/*
 * ***************************************************************************
 * Class AM: Non-inlineable methods (lsolv.c)
 * ***************************************************************************
 */

void AM_lSolve(AM *thee, int prob,
    int meth, int itmax, double etol, int prec, int gues, int pjac);

void AM_hlSolve(AM *thee, int prob,
    int meth, int itmax, double etol, int prec, int gues, int pjac);
#if 0

void AM_hPbcg(AM *thee,
    int lev, int key, int flag, int itmax, double etol, int meth,
    int uu, int ff, int rr, int dd, int ut);
#endif

/*
 * ***************************************************************************
 * Class AM: Non-inlineable methods (nsolv.c)
 * ***************************************************************************
 */

void AM_nSolve(AM *thee,
    int meth, int itmax, double etol,
    int lmeth, int litmax, double letol, int lprec, int gues, int pjac);

void AM_newton(AM *thee,
    int itmax, double etol,
    int lmeth, int litmax, double letol, int lprec, int pjac, double loadParm);

VPUBLIC void AM_homotopy(AM *thee,
    int itmax, double etol,
    int lmeth, int litmax, double letol, int lprec, int pjac);

#endif /* _AM_H_ */

---