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

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;
}
rheolef::space_numbering::ndof
size_type ndof(const basis_basic< T > &b, const geo_size &gs, size_type map_dim)
Definition: space_numbering.cc:28
check_macro
check_macro(expr1.have_homogeneous_space(Xh1), "dual(expr1,expr2); expr1 should have homogeneous space. HINT: use dual(interpolate(Xh, expr1),expr2)")
main
int main(int argc, char **argv)
Definition: vortex_position.cc:27
field
see the field page for the full documentation
rheolef.h
rheolef - reference manual
rheolef::din
idiststream din
see the diststream page for the full documentation
Definition: diststream.h:427
rheolef::environment
see the environment page for the full documentation
Definition: environment.h:115
rheolef
This file is part of Rheolef.
Definition: compiler_eigen.h:37
Float
see the Float 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:459
rheolef::dout
odiststream dout(cout)
see the diststream page for the full documentation
Definition: diststream.h:430