Hi!
I have used the expression Print(u) in the high-level assembly, but the
output is:
Print term Test2_u on Gauss point 0/1 of element 0: sizes (3, 1)
vector(3) [ 0.333333, 0.333333, 0.333333 ]
Should it not print "u" instead of "Test2_u"? Or have I misunderstood
something?
Here is the contents of main():
GETFEM_MPI_INIT(argc, argv);
getfem::mesh mesh;
std::vector<getfem::size_type> ind(3);
ind[0] = mesh.add_point(getfem::base_node(0.0, 0.0));
ind[1] = mesh.add_point(getfem::base_node(1.0, 0.0));
ind[2] = mesh.add_point(getfem::base_node(0.0, 1.0));
mesh.add_simplex(2, ind.begin());
getfem::mesh_fem mesh_fem(mesh);
mesh_fem.set_finite_element(getfem::fem_descriptor("FEM_PK(2,1)"));
getfem::mesh_im mesh_im(mesh);
mesh_im.set_integration_method(getfem::int_method_descriptor("IM_TRIANGLE(1)"));
getfem::ga_workspace workspace;
getfem::base_vector U(mesh_fem.nb_dof());
workspace.add_fem_variable("u", mesh_fem, gmm::sub_interval(0,
mesh_fem.nb_dof()), U);
workspace.add_expression("Test_u*Print(u)", mesh_im);
workspace.assembly(2);
GETFEM_MPI_FINALIZE;
return(0);
Using an updated SVN GetFEM built with gcc (Debian 5.2.1-27) 5.2.1 20151129.
Best regards,
Torquil Sørensen
_______________________________________________
Getfem-users mailing list
[email protected]
https://mail.gna.org/listinfo/getfem-users
