Rheolef  7.1
an efficient C++ finite element environment
vortex_position.cc
Go to the documentation of this file.
1 #include "rheolef.h"
26 using namespace rheolef;
27 int 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 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:459
see the environment page for the full documentation
Definition: environment.h:115
idiststream din
see the diststream page for the full documentation
Definition: diststream.h:427
odiststream dout(cout)
see the diststream page for the full documentation
Definition: diststream.h:430
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)