Rheolef  7.2
an efficient C++ finite element environment
 
Loading...
Searching...
No Matches
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
int main()
Definition field2bb.cc:58
check_macro(expr1.have_homogeneous_space(Xh1), "dual(expr1,expr2); expr1 should have homogeneous space. HINT: use dual(interpolate(Xh, expr1),expr2)")
This file is part of Rheolef.
rheolef - reference manual