On Thu, 10 Mar 2011, John Peterson wrote: > The attached patch makes PetscMatrix::print_personal() print to file. > Would you try it out and see if that works for you?
Your patch makes it even more difficult to write to a file in parallel, because my trick to open the file on one processor and use std::cout on the others confuses your method in that way that you create your custom PetscViewer on the first processor but call MatView(...,PETSC_VIEWER_STDOUT_SELF) on all other processors. This causes PETSc to crash. Anyway, I somehow dislike the idea of a temporary file. Also, I still think that a method doing the thing independently of the solver package would be desirable. I have meanwhile implemented my print_sparse() method, see attachment. Changing it to the other API (with bool argument) would be very easy. I agree that the complex version of print() (in the header file) is wrong. Since I have never used libMesh with complex numbers and don't plan to do this at the moment, I will leave this issue to somebody else. Best Regards, Tim -- Dr. Tim Kroeger CeVis -- Center of Complex Systems and Visualization University of Bremen [email protected] Universitaetsallee 29 [email protected] D-28359 Bremen Phone +49-421-218-59246 Germany Fax +49-421-218-59277 ------------------------------------------------------------------------------ Colocation vs. Managed Hosting A question and answer guide to determining the best fit for your organization - today and in the future. http://p.sf.net/sfu/internap-sfd2d _______________________________________________ Libmesh-devel mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-devel
