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.
Is there a demo or a piece of code showing how to use
MeshQuality::radius_ratio_matplolib_histogram ?
Jan
_______________________________________________
fenics mailing list
[email protected]
http://fenicsproject.org/mailman/listinfo/fenics
