Rheolef  7.2
an efficient C++ finite element environment
vortex_position.cc

The stream function minima and its position.

The stream function minima and its position

#include "rheolef.h"
using namespace rheolef;
int main (int argc, char** argv) {
environment rheolef (argc, argv);
check_macro (communicator().size() == 1, "please, use sequentially");
field psi_h;
din >> psi_h;
size_t idof_min = 0;
Float psi_min = std::numeric_limits<Float>::max();
for (size_t idof = 0, ndof = psi_h.ndof(); idof < ndof; idof++) {
if (psi_h.dof(idof) >= psi_min) continue;
psi_min = psi_h.dof(idof);
idof_min = idof;
}
const disarray<point>& xdof = psi_h.get_space().get_xdofs();
point xmin = xdof [idof_min];
dout << "xc\t\tyc\t\tpsi" << std::endl
<< xmin[0] << "\t" << xmin[1] << "\t" << psi_min << std::endl;
}
see the Float page for the full documentation
see the communicator page for the full documentation
see the field page for the full documentation
see the point page for the full documentation
see the disarray page for the full documentation
Definition: disarray.h:497
see the environment page for the full documentation
Definition: environment.h:121
check_macro(expr1.have_homogeneous_space(Xh1), "dual(expr1,expr2); expr1 should have homogeneous space. HINT: use dual(interpolate(Xh, expr1),expr2)")
size_type ndof(const basis_basic< T > &b, const geo_size &gs, size_type map_dim)
This file is part of Rheolef.
rheolef - reference manual
int main(int argc, char **argv)