|
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 | ||
| ) |
1.8.3