int TestPeriodic(void) {
bool test=true;
field f;
init_empty_field(&f);
f.model.m=_INDEX_MAX; // num of conservative variables
f.vmax = _VMAX; // maximal wave speed
f.model.NumFlux=VlasovP_Lagrangian_NumFlux;
f.model.Source = NULL;
f.model.BoundaryFlux = TestPeriodic_BoundaryFlux;
f.model.InitData = TestPeriodic_InitData;
f.model.ImposedData = TestPeriodic_ImposedData;
f.varindex=GenericVarindex;
f.pre_dtfield=NULL;
f.post_dtfield=NULL;
f.update_after_rk=NULL;
f.model.cfl=0.05;
f.interp.interp_param[0]=f.model.m; // _M
f.interp.interp_param[1]=3; // x direction degree
f.interp.interp_param[2]=0; // y direction degree
f.interp.interp_param[3]=0; // z direction degree
f.interp.interp_param[4]=10; // x direction refinement
f.interp.interp_param[5]=1; // y direction refinement
f.interp.interp_param[6]=1; // z direction refinement
// read the gmsh file
ReadMacroMesh(&(f.macromesh), "test/testcube.msh");
// try to detect a 2d mesh
Detect1DMacroMesh(&(f.macromesh));
assert(f.macromesh.is1d);
// mesh preparation
f.macromesh.period[0]=1;
BuildConnectivity(&(f.macromesh));
PrintMacroMesh(&(f.macromesh));
//assert(1==2);
//AffineMapMacroMesh(&(f.macromesh));
// prepare the initial fields
Initfield(&f);
f.nb_diags = 0;
// prudence...
CheckMacroMesh(&(f.macromesh),f.interp.interp_param+1);
printf("cfl param =%f\n",f.hmin);
// time derivative
//dtField(&f);
//DisplayField(&f);
//assert(1==2);
// apply the DG scheme
// time integration by RK2 scheme
// up to final time = 1.
//RK2(&f,0.5,0.1);
f.vmax=_VMAX;
real dt = set_dt(&f);
RK2(&f,0.5, dt);
// save the results and the error
Plotfield(0,(1==0),&f,"sol","dgvisu.msh");
Plotfield(0,(1==1),&f,"error","dgerror.msh");
real dd=L2error(&f);
real dd_Kinetic=L2_Kinetic_error(&f);
printf("erreur kinetic L2=%lf\n",dd_Kinetic);
printf("erreur L2=%lf\n",dd);
test= test && (dd<3e-3);
//SolvePoisson(&f);
return test;
}
int Test_TransportVP()
{
bool test = true;
field f;
init_empty_field(&f);
int vec=1;
f.model.m=_INDEX_MAX; // num of conservative variables f(vi) for
// each vi, phi, E, rho, u, p, e (ou T)
f.model.NumFlux=VlasovP_Lagrangian_NumFlux;
//f.model.Source = NULL;
f.model.InitData = Test_TransportVP_InitData;
f.model.ImposedData = Test_TransportVP_ImposedData;
f.model.BoundaryFlux = Test_TransportVP_BoundaryFlux;
f.varindex = GenericVarindex;
f.interp.interp_param[0] = f.model.m; // _M
f.interp.interp_param[1] = 2; // x direction degree
f.interp.interp_param[2] = 0; // y direction degree
f.interp.interp_param[3] = 0; // z direction degree
f.interp.interp_param[4] = 16; // x direction refinement
f.interp.interp_param[5] = 1; // y direction refinement
f.interp.interp_param[6] = 1; // z direction refinement
// read the gmsh file
ReadMacroMesh(&(f.macromesh), "../test/testcube.msh");
// try to detect a 2d mesh
Detect1DMacroMesh(&(f.macromesh));
bool is1d = f.macromesh.is1d;
assert(is1d);
// mesh preparation
BuildConnectivity(&(f.macromesh));
//AffineMapMacroMesh(&(f.macromesh));
// prepare the initial fields
f.model.cfl = 0.05;
Initfield(&f);
f.vmax = _VMAX; // maximal wave speed
f.nb_diags = 3;
f.pre_dtfield = UpdateVlasovPoisson;
f.post_dtfield=NULL;
f.update_after_rk = PlotVlasovPoisson;
f.model.Source = VlasovP_Lagrangian_Source;
// prudence...
CheckMacroMesh(&(f.macromesh), f.interp.interp_param + 1);
printf("cfl param =%f\n", f.hmin);
real tmax = 0.03;
real dt = set_dt(&f);
RK2(&f, tmax, dt);
//RK2(&f,0.03,0.05);
// save the results and the error
int iel = 2 * _NB_ELEM_V / 3;
int iloc = _DEG_V;
printf("Trace vi=%f\n", -_VMAX + iel * _DV + _DV * glop(_DEG_V, iloc));
Plotfield(iloc + iel * _DEG_V, false, &f, "sol","dgvisu.msh");
Plotfield(iloc + iel * _DEG_V, true, &f, "error","dgerror.msh");
Plot_Energies(&f, dt);
real dd_Kinetic = L2_Kinetic_error(&f);
printf("erreur kinetic L2=%lf\n", dd_Kinetic);
test= test && (dd_Kinetic < 1e-2);
return test;
}
