CG
Version 25
|
#include "Cgsm.h"
#include "display.h"
#include "CompositeGridOperators.h"
#include "ParallelUtility.h"
#include "ParallelOverlappingGridInterpolator.h"
#include "SmParameters.h"
#include "GridFunctionFilter.h"
#include "GridMaterialProperties.h"
Macros | |
#define | advSM EXTERN_C_NAME(advsm) |
#define | FD4_2D(u, i1, i2, i3, c) ( -( u(i1-2,i2,i3,c)+u(i1+2,i2,i3,c)+u(i1,i2-2,i3,c)+u(i1,i2+2,i3,c) ) +4.*( u(i1-1,i2,i3,c)+u(i1+1,i2,i3,c)+u(i1,i2-1,i3,c)+u(i1,i2+1,i3,c) ) -12.*u(i1,i2,i3,c) ) |
#define | FD4_3D(u, i1, i2, i3, c) ( -( u(i1-2,i2,i3,c)+u(i1+2,i2,i3,c)+u(i1,i2-2,i3,c)+u(i1,i2+2,i3,c)+u(i1,i2,i3-2,c)+u(i1,i2,i3+2,c) ) +4.*( u(i1-1,i2,i3,c)+u(i1+1,i2,i3,c)+u(i1,i2-1,i3,c)+u(i1,i2+1,i3,c)+u(i1,i2,i3-1,c)+u(i1,i2,i3+1,c) ) -18.*u(i1,i2,i3,c) ) |
#define | FOR_3(i1, i2, i3, I1, I2, I3) for( i3=I3.getBase(); i3<=I3.getBound(); i3++ ) for( i2=I2.getBase(); i2<=I2.getBound(); i2++ ) for( i1=I1.getBase(); i1<=I1.getBound(); i1++ ) |
#define | FN(m) fn[m+numberOfFunctions*(grid)] |
Functions | |
void | advSM (const int &nd, const int &n1a, const int &n1b, const int &n2a, const int &n2b, const int &n3a, const int &n3b, const int &nd1a, const int &nd1b, const int &nd2a, const int &nd2b, const int &nd3a, const int &nd3b, const int &nd4a, const int &nd4b, const int &mask, const real &rx, const real &xy, const real &um, const real &u, real &un, const real &f, const int &ndMatProp, const int &matIndex, const real &matValpc, const real &matVal, const int &bc, const real &dis, const real &varDis, const int &ipar, const real &rpar, int &ierr) |
#define advSM EXTERN_C_NAME(advsm) |
Referenced by Cgsm::advanceSOS().
#define FD4_2D | ( | u, | |
i1, | |||
i2, | |||
i3, | |||
c | |||
) | ( -( u(i1-2,i2,i3,c)+u(i1+2,i2,i3,c)+u(i1,i2-2,i3,c)+u(i1,i2+2,i3,c) ) +4.*( u(i1-1,i2,i3,c)+u(i1+1,i2,i3,c)+u(i1,i2-1,i3,c)+u(i1,i2+1,i3,c) ) -12.*u(i1,i2,i3,c) ) |
#define FD4_3D | ( | u, | |
i1, | |||
i2, | |||
i3, | |||
c | |||
) | ( -( u(i1-2,i2,i3,c)+u(i1+2,i2,i3,c)+u(i1,i2-2,i3,c)+u(i1,i2+2,i3,c)+u(i1,i2,i3-2,c)+u(i1,i2,i3+2,c) ) +4.*( u(i1-1,i2,i3,c)+u(i1+1,i2,i3,c)+u(i1,i2-1,i3,c)+u(i1,i2+1,i3,c)+u(i1,i2,i3-1,c)+u(i1,i2,i3+1,c) ) -18.*u(i1,i2,i3,c) ) |
#define FOR_3 | ( | i1, | |
i2, | |||
i3, | |||
I1, | |||
I2, | |||
I3 | |||
) | for( i3=I3.getBase(); i3<=I3.getBound(); i3++ ) for( i2=I2.getBase(); i2<=I2.getBound(); i2++ ) for( i1=I1.getBase(); i1<=I1.getBound(); i1++ ) |
void advSM | ( | const int & | nd, |
const int & | n1a, | ||
const int & | n1b, | ||
const int & | n2a, | ||
const int & | n2b, | ||
const int & | n3a, | ||
const int & | n3b, | ||
const int & | nd1a, | ||
const int & | nd1b, | ||
const int & | nd2a, | ||
const int & | nd2b, | ||
const int & | nd3a, | ||
const int & | nd3b, | ||
const int & | nd4a, | ||
const int & | nd4b, | ||
const int & | mask, | ||
const real & | rx, | ||
const real & | xy, | ||
const real & | um, | ||
const real & | u, | ||
real & | un, | ||
const real & | f, | ||
const int & | ndMatProp, | ||
const int & | matIndex, | ||
const real & | matValpc, | ||
const real & | matVal, | ||
const int & | bc, | ||
const real & | dis, | ||
const real & | varDis, | ||
const int & | ipar, | ||
const real & | rpar, | ||
int & | ierr | ||
) |