From: Andrew E S. <and...@gm...> - 2012-07-16 15:13:46
|
I would like to view the the solution vector of a TransientExplicitSystem to test that some of my calculations are working correctly. I have tried the following, but get a segmentation fault. AutoPtr<NumericVector<Number> > ptr = system.current_local_solution; ptr->print(std::cout); It prints 9 values, then faults (my test has a single QUAD4 element). The first value is correct for the first node of the first variable, the rest of the numbers are garbage. If I loop through the nodes and use point_value(...) things work fine. I am probably just not using the pointer correctly. I would appreciate any help. Thanks, Andrew -- Andrew E. Slaughter, PhD and...@gm... Materials Process Design and Control Laboratory Sibley School of Mechanical and Aerospace Engineering 169 Frank H. T. Rhodes Hall Cornell University Ithaca, NY 14853-3801 (607) 229-1829 http://aeslaughter.github.com/ |