I have a question about the initialization of a petsc_matrix.
The init call is
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?
Thanks!
Jens
------------------------------------------------------------------------------
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
