rheolef  6.5
reference_element2_tst.cc
Go to the documentation of this file.
1 #include "rheolef/reference_element.h"
2 //#include "rheolef/geo_element.h"
3 using namespace rheolef;
4 using namespace std;
5 
6 void
7 put_prism (ostream& out, size_t my_order)
8 {
10  typedef point_basic<size_type> ilat;
11  for (size_type k = 0; k <= my_order; k++) {
12  for (size_type j = 0; j <= my_order; j++) {
13  for (size_type i = 0; i+j <= my_order; i++) {
14  size_type loc_inod = reference_element_P::ilat2loc_inod (my_order, ilat(i, j, k));
15  out << " loc_inod("<<i<<","<<j<<","<<k<<") = " << loc_inod << endl;
16  }
17  }
18  }
19 }
20 int main() {
22  cout << "P1(P):"<<endl; put_prism (cout, 1);
23  cout << "P2(P):"<<endl; put_prism (cout, 2);
24  cout << "P3(P):"<<endl; put_prism (cout, 3);
25  return 0;
26 }
27 
28