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);
------------------------------------------------------------------------------
This SF email is sponsosred by:
Try Windows Azure free for 90 days Click Here
http://p.sf.net/sfu/sfd2d-msazure
_______________________________________________
Libmesh-users mailing list
[email protected]
https://lists.sourceforge.net/lists/listinfo/libmesh-users
