[Libmesh-users] Grid Projections on Lonestar

[Andrea Hawkins-D. <and...@gm...>]

Hello-

I'm attempting to use a grid projection operation on lonestar. It has
worked previously on machines in ICES, (and a colleague, Ernesto, has
verified it works there as of yesterday), but on lonestar we are
getting errors.

It fails because it is looking for a point (0,0,0) that does not exist
on the meshes. I have checked that the solution I'm passing in to
project onto this system is on the domain [12.8,25.6]^2. Also, I svn
up-ed this morning, same issue.

Has anyone else had issues with this functionality on lonestar? Or any
suggestions for what I should check?
Thanks!
Andrea

The (what I think are) relevant lines of code are below:



//Signature for function to be used in the project_solution
Number GridProjection1(const Point& p, const Parameters& parameters,
            const std::string& sys_name,
            const std::string& unknown_name);

Gradient gptr1(const Point& p, const Parameters& parameters,
             const std::string& sys_name,
             const std::string& unknown_name);


....

//Create uniform mesh to project solutions onto
  MeshTools::Generation::build_square (mesh3, 70, 70, 12.8, 25.6,
12.8, 25.6, QUAD9);

//Create systems to hold projected systems
  EquationSystems equation_systems3 (mesh3);

...

function1 = AutoPtr<MeshFunction>
       (new MeshFunction(equation_systems1,
                *f1_soln,
                f1_system.get_dof_map(),
                f1_system.variable_number("u")));
    function1->init();

//Project Solutions
  f3_system.project_solution(GridProjection1,gptr1,equation_systems3.parameters);