On 08/15/2012 05:38 PM, Roy Stogner wrote: > > On Wed, 15 Aug 2012, Jens Lohne Eftang wrote: > >> init (const unsigned int m, >> const unsigned int n, >> const unsigned int m_l, >> const unsigned int n_l, >> const std::vector<unsigned int>& n_nz, >> const std::vector<unsigned int>& n_oz); >> >> where n_nz is the number of non-zeros per row. Hence n_nz.size() = m_l >> on each processor, no? >> >> However there's an assert in petsc_matrix.C, line 110: >> >> libmesh_assert (n_nz.size() == n_l); >> >> What am I missing here? > > It looks to me like you caught a bug. Most of our apps only test > symmetric matrices. Looks like David Knezevic added this code, > though; I'll Cc: him before fixing it.
I agree, looks like a bug. Apologies! David ------------------------------------------------------------------------------ Live Security Virtual Conference Exclusive live event will cover all the ways today's security and threat landscape has changed and how IT managers can respond. Discussions will include endpoint security, mobile security and the latest in malware threats. http://www.accelacomm.com/jaw/sfrnl04242012/114/50122263/ _______________________________________________ Libmesh-users mailing list [email protected] https://lists.sourceforge.net/lists/listinfo/libmesh-users
