[Getfem-users] Print(u) print Test2_u

[Torquil Macdonald Sørensen]

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