94 if ( (
dzpuis == 0) || (tmax(nzm1) ==
double(0)) ) {
96 double s_max =
max( tmax ) ;
98 for (
int l=0; l<nz; l++) {
99 trel.
set(l) = tdiff(l) / s_max ;
106 for (
int l=0; l<nzm1; l++) {
107 s_max = (tmax(l) > s_max) ? tmax(l) : s_max ;
110 for (
int l=0; l<nzm1; l++) {
111 trel.
set(l) = tdiff(l) / s_max ;
114 trel.
set(nzm1) = tdiff(nzm1) / tmax(nzm1) ;
125 for(
int l=0; l<nz; l++) {
126 err.
set(0, l) = trel(l) ;
127 err.
set(1, l) = tdiff(l) ;
128 err.
set(2, l) = tmax(l) ;
135 ostr <<
"Max. source :" ;
136 for (
int l=0; l<nz; l++) {
137 ostr <<
" " << err(2, l) ;
140 ostr << endl <<
"Abs. error : " ;
141 for (
int l=0; l<nz; l++) {
142 ostr <<
" " << err(1, l) ;
146 ostr << endl <<
"Rel. error : " ;
147 for (
int l=0; l<nz; l++) {
148 ostr <<
" " << err(0, l) ;
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
const Scalar & laplacian(int ced_mult_r=4) const
Returns the Laplacian of *this.
double & set(int i)
Read/write of a particular element (index i) (1D case)
Tensor field of valence 0 (or component of a tensorial field).
friend Tbl max(const Scalar &)
Maximum values of a Scalar in each domain.
int dzpuis
Power of r by which the quantity represented by this must be divided in the compactified external do...
Tbl test_poisson(const Scalar &uu, ostream &ostr, bool detail=false) const
Checks if a Poisson equation with *this as a source has been correctly solved.
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
int get_nzone() const
Returns the number of domains.
bool check_dzpuis(int dzi) const
Returns false if the last domain is compactified and *this is not zero in this domain and dzpuis is n...
friend Scalar abs(const Scalar &)
Absolute value.
const Map *const mp
Mapping on which the numerical values at the grid points are defined.
const Map & get_mp() const
Returns the mapping.