static char help[] =
"testing mesh cut test\n\n";
int main(
int argc,
char *argv[]) {
try {
elem_coords(0, 0) = -1;
elem_coords(0, 1) = -1;
elem_coords(0, 2) = -1;
elem_coords(1, 0) = 2;
elem_coords(1, 1) = 0;
elem_coords(1, 2) = 0;
elem_coords(2, 0) = 0;
elem_coords(2, 1) = 1;
elem_coords(2, 2) = 0;
elem_coords(3, 0) = 0;
elem_coords(3, 1) = 0;
elem_coords(3, 2) = 1;
init_local_coords(0, 0) = 0;
init_local_coords(0, 1) = 0;
init_local_coords(0, 2) = 0;
init_local_coords(1, 0) = 0.5;
init_local_coords(1, 1) = 0;
init_local_coords(1, 2) = 0;
init_local_coords(2, 0) = 0;
init_local_coords(2, 1) = 0.5;
init_local_coords(2, 2) = 0;
init_local_coords(3, 0) = 0;
init_local_coords(3, 1) = 0;
init_local_coords(3, 2) = 0.5;
init_local_coords(4, 0) = 1. / 3.;
init_local_coords(4, 1) = 1. / 3.;
init_local_coords(4, 2) = 1. / 3.;
CHKERR Tools::shapeFunMBTET<3>(&shape(0, 0), &init_local_coords(0, 0),
&init_local_coords(0, 1),
&init_local_coords(0, 2), 5);
&elem_coords(0, 0), &global_coords(0, 0), init_local_coords.size1(),
&local_coords(0, 0));
std::cout << residual << std::endl;
for (
auto v : residual.data())
"Should be zer, but is v = %3.4e",
v);
}
return 0;
}
#define CATCH_ERRORS
Catch errors.
@ MOFEM_ATOM_TEST_INVALID
#define CHKERR
Inline error check.
int main(int argc, char *argv[])
double v
phase velocity of light in medium (cm/ns)
UBlasMatrix< double > MatrixDouble
implementation of Data Operators for Forces and Sources
static MoFEMErrorCode Initialize(int *argc, char ***args, const char file[], const char help[])
Initializes the MoFEM database PETSc, MOAB and MPI.
static MoFEMErrorCode Finalize()
Checks for options to be called at the conclusion of the program.