Hi,
Thank you for this fantastic library.
I an trying to write a DG programme, as your introduction, XYZ basis is
better,
but it seems that the function FEXYZ::compute_face_value
can not working here?
hints as follow:
"ERROR logging event compute_face_values()
Did you forget to start or restart it?
[0] src/utils/perf_log.C, line 476, compiled Jan 4 2008 at 14:21:31
p0_15814: p4_error: interrupt SIGx: 6
"
I change to MONOMIAL , it is OK but very slowly.
Code like this , just testing the reinit function:
// The main program.
int main (int argc, char** argv)
{
// Initialize libMesh.
libMesh::init (argc, argv);
{
// Check the input file
GetPot infile("euler.in");
const unsigned int dim=infile("dim", 2);
assert(dim==2);
//input mesh
Mesh mesh (dim);
MeshData mesh_data(mesh);
mesh_data.activate();
const std::string mesh_file = infile("mesh_file", "ideas.msh");
std::cout<<mesh_file;
mesh.read (mesh_file,&mesh_data);
EquationSystems equation_systems (mesh);
TecplotIO tecplot_io(mesh);
mesh.print_info();
//get FE order and family
std::string order;
order= infile("order", "SECOND");
std::string family;
family=infile("family","XYZ");
TransientLinearImplicitSystem & euler_system =
equation_systems.add_system<TransientLinearImplicitSystem> ("euler");
//add conservative variables to the system
euler_system.add_variable("rho",
Utility::string_to_enum<Order> (order),
Utility::string_to_enum<FEFamily>(family));
const unsigned int rho_var = euler_system.variable_number ("rho");
FEType fe_type = euler_system.variable_type(rho_var);
AutoPtr<FEBase> fe_elem_face(FEBase::build(dim, fe_type));
QGauss qface(dim-1, fe_type.default_quadrature_order());
fe_elem_face->attach_quadrature_rule(&qface);
MeshBase::const_element_iterator el = mesh.active_local_elements_begin();
const MeshBase::const_element_iterator end_el =
mesh.active_local_elements_end();
for ( ; el != end_el; ++el)
{
const Elem* elem = *el;
//element boundary quadrature
for (unsigned int side=0; side<elem->n_sides(); side++)
{
//Boundary condition, with Ke and Fe, no offset matrix
if (elem->neighbor(side) == NULL)
{
fe_elem_face->reinit(elem, side);
}
else
{
fe_elem_face->reinit(elem, side);
}
}
}
}//libmesh
// All done.
return libMesh:}
Best regards!
Luyi
2008-04-15
luyi06
2008-04-15
luyi06
-------------------------------------------------------------------------
This SF.net email is sponsored by the 2008 JavaOne(SM) Conference
Don't miss this year's exciting event. There's still time to save $100.
Use priority code J8TL2D2.
http://ad.doubleclick.net/clk;198757673;13503038;p?http://java.sun.com/javaone
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
