On 2013-11-19 11:38, Jan Blechta wrote:
Hi Garth,
would you check line 56 in mesh/MeshQuality.cpp:
44 std::pair<double, double>
MeshQuality::radius_ratio_min_max(const Mesh& mesh)
45 {
46 CellIterator cell(mesh);
47 double qmin = cell->radius_ratio();
48 double qmax = cell->radius_ratio();
49 for (; !cell.end(); ++cell)
50 {
51 qmin = std::min(qmin, cell->radius_ratio());
52 qmax = std::max(qmax, cell->radius_ratio());
53 }
54
55 qmin = MPI::min(qmin);
56 qmax = MPI::min(qmax);
57 return std::make_pair(qmin, qmax);
58 }
This seems like a typo.
Thanks Jan. I've pushed a fix.
Is there a demo or a piece of code showing how to use
MeshQuality::radius_ratio_matplolib_histogram ?
Take a look at the function
test_radius_ratio_matplotlib(self):
in
test/unit/mesh/python/MeshQuality.py
Garth
Jan
_______________________________________________
fenics mailing list
[email protected]
http://fenicsproject.org/mailman/listinfo/fenics
