hi everybody
i'm working with the fourth example, yet i tried to use instead of the
standard LAGRANGE bases the HIERARCHICAL ones (which happen to be C^0
continuous).
i made a slight modification on the part of code in the "main" function
where there is an IF statement ...*
if* ((family == *"LAGRANGE"*) && (order == *"FIRST"*))
and i enforced for the mesh to be constructed through the following command
only...
*MeshTools*::Generation::build_cube (mesh,
ps,
(dim>1) ? ps : 0,
(dim>2) ? ps : 0,
-1., 1.,
-halfwidth, halfwidth,
-halfheight, halfheight,
(dim==1) ? EDGE2 :
((dim == 2) ? QUAD4 : HEX8));
while the rest of code in the ELSE statement is commented!
in case i use the following arguments in the command line:
./ex4-opt -d 2 -n 10 -f HIERARCHIC -o THIRD
then the program aborts. is there any explanation for that??? :)
many thanks in advance
regards,
Vasilis
------------------------------------------------------------------------------
Everyone hates slow websites. So do we.
Make your web apps faster with AppDynamics
Download AppDynamics Lite for free today:
http://p.sf.net/sfu/appdyn_d2d_nov
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
