36 #ifndef NEKTAR_LIB_LIBUTILITIES_DBUTILS_HPP
37 #define NEKTAR_LIB_LIBUTILITIES_DBUTILS_HPP
49 const int stop = StopDefault)
53 ASSERTL1(start < in.num_elements(),
"Start value is outside array range ");
55 if(stop == StopDefault)
58 for(i = start; i < in.num_elements(); ++i)
60 cout << in[i] << endl;
65 ASSERTL1(stop <= in.num_elements(),
"Stop value is outside array range ");
67 for(i = start; i < stop; ++i)
69 cout << in[i] << endl;
76 const int stop = StopDefault)
80 ASSERTL1(start < in.num_elements(),
"Start value is outside array range ");
82 ofstream ofile(outfile.c_str());
84 if(stop == StopDefault)
87 for(i = start; i < in.num_elements(); ++i)
89 ofile << in[i] << endl;
94 ASSERTL1(stop <= in.num_elements(),
"Stop value is outside array range ");
96 for(i = start; i < stop; ++i)
98 ofile << in[i] << endl;
106 const int stop = StopDefault)
110 ASSERTL1(start < in.num_elements(),
"Start value is outside array range ");
112 if(stop == StopDefault)
115 for(i = start; i < in.num_elements(); ++i)
117 ofile << in[i] << endl;
122 ASSERTL1(stop <= in.num_elements(),
"Stop value is outside array range ");
124 for(i = start; i < stop; ++i)
126 ofile << in[i] << endl;
134 const int stop = StopDefault)
140 if(stop == StopDefault)
145 ofile << in[i] << endl;
152 for(i = start; i < stop; ++i)
154 ofile << in[i] << endl;
169 out <<
"Norm: " << vExchange[0] << endl;
void NormGlobalVector(const int n, Array< OneD, const T > &in, std::ostream &out, MultiRegions::AssemblyMapCGSharedPtr &map)
unsigned int GetDimension() const
Returns the number of dimensions for the point.
void Output1DArray(const Array< OneD, const T > &in, const int start=0, const int stop=StopDefault)
T Dot2(int n, const T *w, const T *x, const int *y)
vvtvp (vector times vector times vector): z = w*x*y
boost::shared_ptr< AssemblyMapCG > AssemblyMapCGSharedPtr
#define ASSERTL1(condition, msg)
Assert Level 1 – Debugging which is used whether in FULLDEBUG or DEBUG compilation mode...