CG
Version 25
|
#include "Maxwell.h"
#include "CompositeGridOperators.h"
#include "display.h"
#include "UnstructuredMapping.h"
#include "OGPolyFunction.h"
#include "OGTrigFunction.h"
#include "OGPulseFunction.h"
#include "RadiationBoundaryCondition.h"
#include "ParallelUtility.h"
#include "ParallelGridUtility.h"
Macros | |
#define | interfaceMaxwell EXTERN_C_NAME(interfacemaxwell) |
#define | newInterfaceMaxwell EXTERN_C_NAME(newinterfacemaxwell) |
#define | interface3dMaxwell EXTERN_C_NAME(interface3dmaxwell) |
#define | FOR_3D(i1, i2, i3, I1, I2, I3) int I1Base =I1.getBase(), I2Base =I2.getBase(), I3Base =I3.getBase(); int I1Bound=I1.getBound(), I2Bound=I2.getBound(), I3Bound=I3.getBound(); for(i3=I3Base; i3<=I3Bound; i3++) for(i2=I2Base; i2<=I2Bound; i2++) for(i1=I1Base; i1<=I1Bound; i1++) |
#define | FOR_3(i1, i2, i3, I1, I2, I3) I1Base =I1.getBase(), I2Base =I2.getBase(), I3Base =I3.getBase(); I1Bound=I1.getBound(), I2Bound=I2.getBound(), I3Bound=I3.getBound(); for(i3=I3Base; i3<=I3Bound; i3++) for(i2=I2Base; i2<=I2Bound; i2++) for(i1=I1Base; i1<=I1Bound; i1++) |
#define | FOR_3IJD(i1, i2, i3, I1, I2, I3, j1, j2, j3, J1, J2, J3) int I1Base =I1.getBase(), I2Base =I2.getBase(), I3Base =I3.getBase(); int I1Bound=I1.getBound(), I2Bound=I2.getBound(), I3Bound=I3.getBound(); int J1Base =J1.getBase(), J2Base =J2.getBase(), J3Base =J3.getBase(); for(i3=I3Base,j3=J3Base; i3<=I3Bound; i3++,j3++) for(i2=I2Base,j2=J2Base; i2<=I2Bound; i2++,j2++) for(i1=I1Base,j1=J1Base; i1<=I1Bound; i1++,j1++) |
#define | FOR_4IJD(i1, i2, i3, i4, I1, I2, I3, I4, j1, j2, j3, j4, J1, J2, J3, J4) int I1Base =I1.getBase(), I2Base =I2.getBase(), I3Base =I3.getBase(), I4Base =I4.getBase(); int I1Bound=I1.getBound(), I2Bound=I2.getBound(), I3Bound=I3.getBound(), I4Bound=I4.getBound(); int J1Base =J1.getBase(), J2Base =J2.getBase(), J3Base =J3.getBase(), J4Base =J4.getBase(); for(i4=I4Base,j4=J4Base; i4<=I4Bound; i4++,j4++) for(i3=I3Base,j3=J3Base; i3<=I3Bound; i3++,j3++) for(i2=I2Base,j2=J2Base; i2<=I2Bound; i2++,j2++) for(i1=I1Base,j1=J1Base; i1<=I1Bound; i1++,j1++) |
#define | FOR_4IJ(i1, i2, i3, i4, I1, I2, I3, I4, j1, j2, j3, j4, J1, J2, J3, J4) I1Base =I1.getBase(), I2Base =I2.getBase(), I3Base =I3.getBase(), I4Base =I4.getBase(); I1Bound=I1.getBound(), I2Bound=I2.getBound(), I3Bound=I3.getBound(), I4Bound=I4.getBound(); J1Base =J1.getBase(), J2Base =J2.getBase(), J3Base =J3.getBase(), J4Base =J4.getBase(); for(i4=I4Base,j4=J4Base; i4<=I4Bound; i4++,j4++) for(i3=I3Base,j3=J3Base; i3<=I3Bound; i3++,j3++) for(i2=I2Base,j2=J2Base; i2<=I2Bound; i2++,j2++) for(i1=I1Base,j1=J1Base; i1<=I1Bound; i1++,j1++) |
#define | FOR_4IJD_WITH_STRIDE(i0, i1, i2, i3, indexi, j0, j1, j2, j3, indexj) for( int i3=indexi[3][0],j3=indexj[3][0]; i3<=indexi[3][1]; i3+=indexi[3][2], j3+=indexj[3][2] )for( int i2=indexi[2][0],j2=indexj[2][0]; i2<=indexi[2][1]; i2+=indexi[2][2], j2+=indexj[2][2] )for( int i1=indexi[1][0],j1=indexj[1][0]; i1<=indexi[1][1]; i1+=indexi[1][2], j1+=indexj[1][2] )for( int i0=indexi[0][0],j0=indexj[0][0]; i0<=indexi[0][1]; i0+=indexi[0][2], j0+=indexj[0][2] ) |
Functions | |
void | interfaceMaxwell (const int &nd, const int &nd1a, const int &nd1b, const int &nd2a, const int &nd2b, const int &nd3a, const int &nd3b, const int &gridIndexRange1, real &u1, const int &mask1, const real &rsxy1, const real &xy1, const int &boundaryCondition1, const int &md1a, const int &md1b, const int &md2a, const int &md2b, const int &md3a, const int &md3b, const int &gridIndexRange2, real &u2, const int &mask2, const real &rsxy2, const real &xy2, const int &boundaryCondition2, const int &ipar, const real &rpar, real &aa2, real &aa4, real &aa8, int &ipvt2, int &ipvt4, int &ipvt8, int &ierr) |
void | newInterfaceMaxwell (const int &nd, const int &nd1a, const int &nd1b, const int &nd2a, const int &nd2b, const int &nd3a, const int &nd3b, const int &gridIndexRange1, real &u1, const int &mask1, const real &rsxy1, const real &xy1, const int &boundaryCondition1, const int &md1a, const int &md1b, const int &md2a, const int &md2b, const int &md3a, const int &md3b, const int &gridIndexRange2, real &u2, const int &mask2, const real &rsxy2, const real &xy2, const int &boundaryCondition2, const int &ipar, const real &rpar, int &ierr) |
void | interface3dMaxwell (const int &nd, const int &nd1a, const int &nd1b, const int &nd2a, const int &nd2b, const int &nd3a, const int &nd3b, const int &gridIndexRange1, real &u1, const int &mask1, const real &rsxy1, const real &xy1, const int &boundaryCondition1, const int &md1a, const int &md1b, const int &md2a, const int &md2b, const int &md3a, const int &md3b, const int &gridIndexRange2, real &u2, const int &mask2, const real &rsxy2, const real &xy2, const int &boundaryCondition2, const int &ipar, const real &rpar, real &aa2, real &aa4, real &aa8, int &ipvt2, int &ipvt4, int &ipvt8, int &ierr) |
#define FOR_3 | ( | i1, | |
i2, | |||
i3, | |||
I1, | |||
I2, | |||
I3 | |||
) | I1Base =I1.getBase(), I2Base =I2.getBase(), I3Base =I3.getBase(); I1Bound=I1.getBound(), I2Bound=I2.getBound(), I3Bound=I3.getBound(); for(i3=I3Base; i3<=I3Bound; i3++) for(i2=I2Base; i2<=I2Bound; i2++) for(i1=I1Base; i1<=I1Bound; i1++) |
#define FOR_3D | ( | i1, | |
i2, | |||
i3, | |||
I1, | |||
I2, | |||
I3 | |||
) | int I1Base =I1.getBase(), I2Base =I2.getBase(), I3Base =I3.getBase(); int I1Bound=I1.getBound(), I2Bound=I2.getBound(), I3Bound=I3.getBound(); for(i3=I3Base; i3<=I3Bound; i3++) for(i2=I2Base; i2<=I2Bound; i2++) for(i1=I1Base; i1<=I1Bound; i1++) |
#define FOR_3IJD | ( | i1, | |
i2, | |||
i3, | |||
I1, | |||
I2, | |||
I3, | |||
j1, | |||
j2, | |||
j3, | |||
J1, | |||
J2, | |||
J3 | |||
) | int I1Base =I1.getBase(), I2Base =I2.getBase(), I3Base =I3.getBase(); int I1Bound=I1.getBound(), I2Bound=I2.getBound(), I3Bound=I3.getBound(); int J1Base =J1.getBase(), J2Base =J2.getBase(), J3Base =J3.getBase(); for(i3=I3Base,j3=J3Base; i3<=I3Bound; i3++,j3++) for(i2=I2Base,j2=J2Base; i2<=I2Bound; i2++,j2++) for(i1=I1Base,j1=J1Base; i1<=I1Bound; i1++,j1++) |
#define FOR_4IJ | ( | i1, | |
i2, | |||
i3, | |||
i4, | |||
I1, | |||
I2, | |||
I3, | |||
I4, | |||
j1, | |||
j2, | |||
j3, | |||
j4, | |||
J1, | |||
J2, | |||
J3, | |||
J4 | |||
) | I1Base =I1.getBase(), I2Base =I2.getBase(), I3Base =I3.getBase(), I4Base =I4.getBase(); I1Bound=I1.getBound(), I2Bound=I2.getBound(), I3Bound=I3.getBound(), I4Bound=I4.getBound(); J1Base =J1.getBase(), J2Base =J2.getBase(), J3Base =J3.getBase(), J4Base =J4.getBase(); for(i4=I4Base,j4=J4Base; i4<=I4Bound; i4++,j4++) for(i3=I3Base,j3=J3Base; i3<=I3Bound; i3++,j3++) for(i2=I2Base,j2=J2Base; i2<=I2Bound; i2++,j2++) for(i1=I1Base,j1=J1Base; i1<=I1Bound; i1++,j1++) |
#define FOR_4IJD | ( | i1, | |
i2, | |||
i3, | |||
i4, | |||
I1, | |||
I2, | |||
I3, | |||
I4, | |||
j1, | |||
j2, | |||
j3, | |||
j4, | |||
J1, | |||
J2, | |||
J3, | |||
J4 | |||
) | int I1Base =I1.getBase(), I2Base =I2.getBase(), I3Base =I3.getBase(), I4Base =I4.getBase(); int I1Bound=I1.getBound(), I2Bound=I2.getBound(), I3Bound=I3.getBound(), I4Bound=I4.getBound(); int J1Base =J1.getBase(), J2Base =J2.getBase(), J3Base =J3.getBase(), J4Base =J4.getBase(); for(i4=I4Base,j4=J4Base; i4<=I4Bound; i4++,j4++) for(i3=I3Base,j3=J3Base; i3<=I3Bound; i3++,j3++) for(i2=I2Base,j2=J2Base; i2<=I2Bound; i2++,j2++) for(i1=I1Base,j1=J1Base; i1<=I1Bound; i1++,j1++) |
#define FOR_4IJD_WITH_STRIDE | ( | i0, | |
i1, | |||
i2, | |||
i3, | |||
indexi, | |||
j0, | |||
j1, | |||
j2, | |||
j3, | |||
indexj | |||
) | for( int i3=indexi[3][0],j3=indexj[3][0]; i3<=indexi[3][1]; i3+=indexi[3][2], j3+=indexj[3][2] )for( int i2=indexi[2][0],j2=indexj[2][0]; i2<=indexi[2][1]; i2+=indexi[2][2], j2+=indexj[2][2] )for( int i1=indexi[1][0],j1=indexj[1][0]; i1<=indexi[1][1]; i1+=indexi[1][2], j1+=indexj[1][2] )for( int i0=indexi[0][0],j0=indexj[0][0]; i0<=indexi[0][1]; i0+=indexi[0][2], j0+=indexj[0][2] ) |
#define interface3dMaxwell EXTERN_C_NAME(interface3dmaxwell) |
Referenced by Maxwell::assignInterfaceBoundaryConditions().
#define interfaceMaxwell EXTERN_C_NAME(interfacemaxwell) |
Referenced by Maxwell::assignInterfaceBoundaryConditions().
#define newInterfaceMaxwell EXTERN_C_NAME(newinterfacemaxwell) |
Referenced by Maxwell::assignInterfaceBoundaryConditions().
void interface3dMaxwell | ( | const int & | nd, |
const int & | nd1a, | ||
const int & | nd1b, | ||
const int & | nd2a, | ||
const int & | nd2b, | ||
const int & | nd3a, | ||
const int & | nd3b, | ||
const int & | gridIndexRange1, | ||
real & | u1, | ||
const int & | mask1, | ||
const real & | rsxy1, | ||
const real & | xy1, | ||
const int & | boundaryCondition1, | ||
const int & | md1a, | ||
const int & | md1b, | ||
const int & | md2a, | ||
const int & | md2b, | ||
const int & | md3a, | ||
const int & | md3b, | ||
const int & | gridIndexRange2, | ||
real & | u2, | ||
const int & | mask2, | ||
const real & | rsxy2, | ||
const real & | xy2, | ||
const int & | boundaryCondition2, | ||
const int & | ipar, | ||
const real & | rpar, | ||
real & | aa2, | ||
real & | aa4, | ||
real & | aa8, | ||
int & | ipvt2, | ||
int & | ipvt4, | ||
int & | ipvt8, | ||
int & | ierr | ||
) |
void interfaceMaxwell | ( | const int & | nd, |
const int & | nd1a, | ||
const int & | nd1b, | ||
const int & | nd2a, | ||
const int & | nd2b, | ||
const int & | nd3a, | ||
const int & | nd3b, | ||
const int & | gridIndexRange1, | ||
real & | u1, | ||
const int & | mask1, | ||
const real & | rsxy1, | ||
const real & | xy1, | ||
const int & | boundaryCondition1, | ||
const int & | md1a, | ||
const int & | md1b, | ||
const int & | md2a, | ||
const int & | md2b, | ||
const int & | md3a, | ||
const int & | md3b, | ||
const int & | gridIndexRange2, | ||
real & | u2, | ||
const int & | mask2, | ||
const real & | rsxy2, | ||
const real & | xy2, | ||
const int & | boundaryCondition2, | ||
const int & | ipar, | ||
const real & | rpar, | ||
real & | aa2, | ||
real & | aa4, | ||
real & | aa8, | ||
int & | ipvt2, | ||
int & | ipvt4, | ||
int & | ipvt8, | ||
int & | ierr | ||
) |
void newInterfaceMaxwell | ( | const int & | nd, |
const int & | nd1a, | ||
const int & | nd1b, | ||
const int & | nd2a, | ||
const int & | nd2b, | ||
const int & | nd3a, | ||
const int & | nd3b, | ||
const int & | gridIndexRange1, | ||
real & | u1, | ||
const int & | mask1, | ||
const real & | rsxy1, | ||
const real & | xy1, | ||
const int & | boundaryCondition1, | ||
const int & | md1a, | ||
const int & | md1b, | ||
const int & | md2a, | ||
const int & | md2b, | ||
const int & | md3a, | ||
const int & | md3b, | ||
const int & | gridIndexRange2, | ||
real & | u2, | ||
const int & | mask2, | ||
const real & | rsxy2, | ||
const real & | xy2, | ||
const int & | boundaryCondition2, | ||
const int & | ipar, | ||
const real & | rpar, | ||
int & | ierr | ||
) |