space.h

Go to the documentation of this file.

/* -----------------------------------------------------------------------------------------
  $Id: space.h 573 2006-11-20 00:25:28Z svacek $
  $Author: svacek $
   ----------------------------------------------------------------------------------------- */
#ifndef __SPACE_SVH
#define __SPACE_SVH
/* -------------------------------------------------------------------------- */
/*@{*/
/** \addtogroup FEM_SPACE_DEFINITION */
/* function SetElement()  - nastaveni paramatru pro ruzne typy elementu  */
void SetElement();                                           /**< Set type of GLOBAL element. Global element should be only the maximal one */
/* function ConstructSpace(),DestructSpace() */

/** \addtogroup SPACE_IO 
 *   \ingroup    FEM_SPACE_DEFINITION
 */
/*@{*/
void SaveProgramState(char *FileName);                       /**< Saves all what is needed for computation for future computation restart. */
void LoadProgramState(char *FileName);                       /**< Loads all what is needed for computation for computation restart. */
void SaveProgramStateMin(char *FileName);                    /**< Saves minimal program state. */
void LoadProgramStateMin(char *FileName);                    /**< Loads minimal program state. */
void LoadNodeApproximation(const char *filename,vector *_load); /**< Load node (P1) approximation of solution */


void SaveAllOutputs(const long iter,const vector pom);       /**< Save all outputs (according to INI file settings) */
void GMVOutput(char filename[],vector pom);                  /**< Save approximation pom (UVP) to filename[] in GMV format*/
void TecplotOutput(char filename[],int delka,char *nazvy, vector seznam[]); /**< not documented (NOT IN CURRENT USE) \todo Remove this function */
void TecplotOutputR(const char fname[], vector pom,int How);                /**< Save pom(UVP) to TECPLOT format with How=0 - only P1 appr. or higher with How=1 */
void ElasticOutput(const char fname[],const vector X,const vector XDOT,const vector XDOTDOT);                          /**< Save Elastic Deformation */
void ElasticOutput1(const char fname[],const vector X,const vector XDOT,const vector XDOTDOT);                          /**< Save Elastic Deformation */
void ResultsOutput(char filename[], vector pom);             /**< Output in columns X Y U V P ...., usable for MATLAB, ANGENER mesh refinement, etc. */

void SaveBndrPressureInterpolate(const char fname[],const vector u);                 /**< Save PRESSURE on surface. */
double ControlStep(const vector u1,const vector u2,const long ino,const double ctm); /**< Control step - residuals, forces control etc. */
/*@}*/

/** 
 *   \defgroup   SPACE_INITIALIZATION Space Initialization
 *   \ingroup    FEM_SPACE_DEFINITION
 *   \addtogroup SPACE_INITIALIZATION 
 */
/*@{*/
void ConstructSpace(double xvel(double,double),double yvel(double,double));    /**< Allocate and fills all data of FEM space definition (velocity space). */
void DestructSpace();                                        /**< De-allocate all data of FEM space definition. */
void ConstructPressureSpace();                               /**< Allocate and fills all data of FEM space definition (pressure space). */

void ConstructBoundaryConditions(double xvel(double,double),double yvel(double,double)); /**< P1 approximation of boundary condition prescribed by functions */
void DestructBoundaryConditions();                                                       /**< destruct boundary condition representation */

/* velocity space */
short GetMark(long SideIdx);                                /**< Gets Mark for a side (velocity) */
void  ComputeBaseIndices();                                 /**< Maps LOCAL(element) base functions to GLOBAL(grid) - velocity */
void  ComputeBaseValues();                                  /**< Computes values on reference element (at Quadrature nodes)*/
void  ComputeGradBaseValues();                              /**< Computes GRADIENT values on reference element (at Quadrature nodes 2D) */
void  ComputeHesBaseValues();                               /**< Computes HESSIAN values on reference element (at Quadrature nodes 2D) */
void  ComputeSideBaseValues();                              /**< Computes base values on reference side (at Quadrature nodes 1D) */ 
void  ComputeExportValues();                                /**< Computes values for refined output, used in ::TecplotOuputR */
void  DestructExportValues();                               /**< Deallocate Export Values */
/* pressure space */
short GetPressureMark(long SideIdx);                        /**< Gets Mark for a side (pressure) */
void  ComputePressureBaseIndices();                         /**< Maps LOCAL(element) base functions to GLOBAL(grid) - pressure */
long  SelectFixedPressureIdx();                             /**< select a point where pressure should be fixed, \todo Remove this - not needed! */


void InitialVelocities();                                    /**< Fills up first approximation of the solution. */
void ConstructVectors();                                     /**< Constructs all vectors needed for computation. */
void DestructVectors();                                      /**< Destructs all vectors needed for computation. */
void ExtrapolateVelocities();                                /**< Extrapolate Fluid Velocities */
/*@}*/
void Project2HDIV(vector ww);                                /**< Projection of ww to H(div) space */
void HDivProjection(sparse *S,vector b,const vector ww);     /**< Assembly of the linear system for H(div) projection S*ww = b */
/* kvadratura */
void QuadratureFormula1d(double *W,double *X,int q1);        /**< Fill Up 1D quadr. formula with degree of 2*q1-1 */                    
int  QuadratureFormulae2dTriGauss(double **w,double ***L,int q); /**< Fill Up 2D quadr. formula with degree of 2*q */                    
int  QuadratureFormulae2dTriangle(double **w,double ***L,int q); /**< Fill Up 2D quadr. formula with degree of q \todo Remove this function */                    
void InitQuadratureFormulae();                               /**< Init Q.F */
void InitQuadratureFormulaePrecise();                        /**< Init Q.F.P. */
/* konverze */
void Pressure2Velocity(vector out,const vector press);       /**< Convert pressure to velocity */
void Velocity2Pressure(vector out,const vector vel);         /**< Convert velocity to pressure */
void InitPressureInterpolate();

/* ------------------------------------------------------------------------------------- */
/* norms etc. */
/** \addtogroup FORMS_EVALUATION
 *   \ingroup    FEM_SPACE_DEFINITION
 */
/*@{*/
void InitNorms();
void PrepareH1norm(sparse *H1);
void PrepareL2norm(sparse *L2);
void PreparePressureL2norm(sparse *L2);
void PrepareBoundaryForm( vector beta,double _phi(double,double));
void PrepareForms(vector dzeta,double _eta(double,double));
void DestructNorms();
/*@}*/
/* ------------------------------------------------------------------------------------- */
void ControlMarks();
void CheckBoundaryConditions(const char fname[]);
void add2triplet(long _i,long _j,double _val);                /**< add to triplet value (_i,_j,_val) */
double LinearFunctionDX(const long iELE,const double LinF[]); /**< Linear Function x-derivative */
double LinearFunctionDY(const long iELE,const double LinF[]); /**< Linear Function y-derivative */
double LinearFunctionValue(const double LinF[],const int m);  /**< Linear Function Value */
/* ------------------------------------------------------------------------------------- */
/** \addtogroup COARSE_TRIANG 
 *   \ingroup    FEM_SPACE_DEFINITION
 *   \todo Correct that.
 */
/*@{*/
void ConstructCoarse();
void FillCoarseSystem(const sparse S);
void DestructCoarse();
void ConstructCoarsePrecond();
void DestructCoarsePrecond();
void PrecondCoarse(vector x);
/*@}*/
/*@}*/
void InitOriginalCfg();
double MaxLinearError(const vector _new,const vector _old);
double MaxLinearErrorVelocity(const vector _new,const vector _old);
double LinearError(const vector _new,const vector _old);
double LinearErrorVelocity(const vector _new,const vector _old);
double QuadraticError(const vector _new,const vector _old,double *out1,double *out2,double *outp);
double QuadraticErrorComponent(const vector _new,const vector _old);
/* ------------------------------------------------------------------------------------- */
#endif
/* ------------------------------------------------------------------------------------- */

---