Rheolef  7.2
an efficient C++ finite element environment
 
Loading...
Searching...
No Matches
vortex_position.cc
Go to the documentation of this file.
1
25#include "rheolef.h"
26using namespace rheolef;
27int main (int argc, char** argv) {
28 environment rheolef (argc, argv);
29 check_macro (communicator().size() == 1, "please, use sequentially");
30 field psi_h;
31 din >> psi_h;
32 size_t idof_min = 0;
33 Float psi_min = std::numeric_limits<Float>::max();
34 for (size_t idof = 0, ndof = psi_h.ndof(); idof < ndof; idof++) {
35 if (psi_h.dof(idof) >= psi_min) continue;
36 psi_min = psi_h.dof(idof);
37 idof_min = idof;
38 }
39 const disarray<point>& xdof = psi_h.get_space().get_xdofs();
40 point xmin = xdof [idof_min];
41 dout << "xc\t\tyc\t\tpsi" << std::endl
42 << xmin[0] << "\t" << xmin[1] << "\t" << psi_min << std::endl;
43}
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