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;
}
Float
see the Float page for the full documentation
communicator
see the communicator page for the full documentation
field
see the field page for the full documentation
point
see the point page for the full documentation
rheolef::disarray
see the disarray page for the full documentation
Definition disarray.h:497
rheolef::environment
see the environment page for the full documentation
Definition environment.h:121
main
int main()
Definition field2bb.cc:58
check_macro
check_macro(expr1.have_homogeneous_space(Xh1), "dual(expr1,expr2); expr1 should have homogeneous space. HINT: use dual(interpolate(Xh, expr1),expr2)")
rheolef
This file is part of Rheolef.
Definition compiler_eigen.h:39
rheolef.h
rheolef - reference manual