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