On Fri, Aug 27, 2010 at 11:36 AM, John Peterson <peter...@cfdlab.ae.utexas.edu> wrote: > On Fri, Aug 27, 2010 at 11:30 AM, Roy Stogner <royst...@ices.utexas.edu> > wrote: >> >> On Fri, 27 Aug 2010, Johannes A Huber wrote: >> >>> I was wondering, is it might be possible petsc_matric.C on line 72 (and >>> perhaps in other lines as well, like 139, 142, ...) >>> The matrix that is created in these lines is always quadratic >>> (n_global*n_global). I guess, the second argument to MatCreateSeqAIJ should >>> be >>> m_global instead of n_global. >> >> Let me run this by libmesh-devel (Ben Kirk is the one who wrote that >> code, and he surely understands it better than I do), but it looks to >> me like you're right. And since libMesh (and, I'd be willing to bet, >> all its primary developers) only creates square sparse matrices, we'd >> never have noticed a bug that breaks rectangular sparse matrix >> creation. > > I think you're right, especially considering the call to > MatCreateMPIAIJ a few lines below which does use m_global. > > Nice catch!
And fix checked in. -- John ------------------------------------------------------------------------------ Sell apps to millions through the Intel(R) Atom(Tm) Developer Program Be part of this innovative community and reach millions of netbook users worldwide. Take advantage of special opportunities to increase revenue and speed time-to-market. Join now, and jumpstart your future. http://p.sf.net/sfu/intel-atom-d2d _______________________________________________ Libmesh-devel mailing list Libmesh-devel@lists.sourceforge.net https://lists.sourceforge.net/lists/listinfo/libmesh-devel
