scalar_asymptot.C
/*
* Function Scalar::asymptot
*
*/
/*
* Copyright (c) 2003 Eric Gourgoulhon & Jerome Novak
*
* Copyright (c) 1999-2002 Eric Gourgoulhon (for previous class Cmp)
* Copyright (c) 1999-2001 Philippe Grandclement (for previous class Cmp)
*
* This file is part of LORENE.
*
* LORENE is free software; you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation; either version 2 of the License, or
* (at your option) any later version.
*
* LORENE is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with LORENE; if not, write to the Free Software
* Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
*
*/
char scalar_asymptot_C[] = "$Header: /cvsroot/Lorene/C++/Source/Tensor/Scalar/scalar_asymptot.C,v 1.5 2014/10/13 08:53:46 j_novak Exp $" ;
/*
* $Id: scalar_asymptot.C,v 1.5 2014/10/13 08:53:46 j_novak Exp $
* $Log: scalar_asymptot.C,v $
* Revision 1.5 2014/10/13 08:53:46 j_novak
* Lorene classes and functions now belong to the namespace Lorene.
*
* Revision 1.4 2014/10/06 15:16:15 j_novak
* Modified #include directives to use c++ syntax.
*
* Revision 1.3 2004/02/21 17:05:14 e_gourgoulhon
* Method Scalar::point renamed Scalar::val_grid_point.
* Method Scalar::set_point renamed Scalar::set_grid_point.
*
* Revision 1.2 2003/10/08 14:24:09 j_novak
* replaced mult_r_zec with mult_r_ced
*
* Revision 1.1 2003/09/25 07:18:00 j_novak
* Method asymptot implemented.
*
* $Header: /cvsroot/Lorene/C++/Source/Tensor/Scalar/scalar_asymptot.C,v 1.5 2014/10/13 08:53:46 j_novak Exp $
*
*/
// Headers C
#include <cmath>
// Headers Lorene
#include "tensor.h"
namespace Lorene {
Valeur** Scalar::asymptot(int n0, const int flag) const {
assert(n0 >= 0) ;
const Mg3d& mg = *(mp->get_mg()) ;
int nz = mg.get_nzone() ;
int nzm1 = nz-1 ;
assert(mg.get_type_r(nzm1) == UNSURR) ;
int np = mg.get_np(nzm1) ;
int nt = mg.get_nt(nzm1) ;
int nr = mg.get_nr(nzm1) ;
int nrm1 = nr-1 ;
Valeur** vu = new Valeur*[n0+1] ;
for (int h=0; h<=n0; h++) {
vu[h] = new Valeur(mg.get_angu()) ;
}
Scalar uu = *this ;
int precis = 2 ;
// The terms 1/r^h with h < dzpuis are null :
// -----------------------------------------
for (int h=0; h<dzpuis; h++) {
vu[h]->set_etat_zero() ;
}
// Terms 1/r^h with h >= dzpuis :
// -----------------------------
for (int h=dzpuis; h<=n0; h++) {
// Value on the sphere S^2 at infinity
// -----------------------------------
vu[h]->set_etat_c_qcq() ;
vu[h]->c->set_etat_qcq() ;
for (int l=0; l<nzm1; l++) {
vu[h]->c->t[l]->set_etat_zero() ;
}
vu[h]->c->t[nzm1]->set_etat_qcq() ;
for (int k=0; k<np; k++) {
for (int j=0; j<nt; j++) {
vu[h]->set(nzm1, k, j, 0) = uu.val_grid_point(nzm1, k, j, nrm1) ;
}
}
vu[h]->set_base( uu.va.base ) ;
// Printing
// --------
if (flag != 0) {
cout << "Term in 1/r^" << h << endl ;
cout << "-------------" << endl ;
double ndec = 0 ;
double vmin = (*vu[h])(nzm1, 0, 0, 0) ;
double vmax = vmin ;
cout << " Values at the point (phi_k, theta_j) : " << endl ;
cout.precision(precis) ;
cout.setf(ios::showpoint) ;
for (int k=0; k<np; k++) {
cout << " k=" << k << " : " ;
for (int j=0; j<nt; j++) {
double xx = (*vu[h])(nzm1, k, j, 0) ;
cout << " " << setw(precis) << xx ;
ndec += fabs(xx) ;
vmin = ( xx < vmin ) ? xx : vmin ;
vmax = ( xx > vmax ) ? xx : vmax ;
}
cout << endl;
}
ndec /= np*nt ;
cout << "Minimum value on S^2 : " << vmin << endl ;
cout << "Maximum value on S^2 : " << vmax << endl ;
cout << "L^1 norm on S^2 : " << ndec << endl ;
}
// The value at infinity is substracted
// ------------------------------------
for (int k=0; k<np; k++) {
for (int j=0; j<nt; j++) {
double v_inf = (*vu[h])(nzm1, k, j, 0) ;
for (int i=0; i<nr; i++) {
uu.set_grid_point(nzm1, k, j, i) -= v_inf ;
}
}
}
// Mutliplication by r
// -------------------
uu.mult_r_ced() ;
} // End of loop on h (development order)
return vu ;
}
}
Lorene::Mg3d
Multi-domain grid.
Definition grilles.h:273
Lorene::Mg3d::get_angu
const Mg3d * get_angu() const
Returns the pointer on the associated angular grid.
Definition mg3d.C:473
Lorene::Mg3d::get_np
int get_np(int l) const
Returns the number of points in the azimuthal direction ( ) in domain no. l.
Definition grilles.h:462
Lorene::Mg3d::get_nt
int get_nt(int l) const
Returns the number of points in the co-latitude direction ( ) in domain no. l.
Definition grilles.h:457
Lorene::Mg3d::get_nzone
int get_nzone() const
Returns the number of domains.
Definition grilles.h:448
Lorene::Mg3d::get_nr
int get_nr(int l) const
Returns the number of points in the radial direction ( ) in domain no. l.
Definition grilles.h:452
Lorene::Mg3d::get_type_r
int get_type_r(int l) const
Returns the type of sampling in the radial direction in domain no.
Definition grilles.h:474
Lorene::Mtbl::t
Tbl ** t
Array (size nzone ) of pointers on the Tbl 's.
Definition mtbl.h:132
Lorene::Mtbl::set_etat_qcq
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
Definition mtbl.C:299
Lorene::Scalar::val_grid_point
double val_grid_point(int l, int k, int j, int i) const
Returns the value of the field at a specified grid point.
Definition scalar.h:637
Lorene::Scalar::Scalar
Scalar(const Map &mpi)
Constructor from mapping.
Definition scalar.C:204
Lorene::Scalar::mult_r_ced
void mult_r_ced()
Multiplication by r in the compactified external domain (CED), the dzpuis flag is not changed.
Definition scalar_r_manip.C:269
Lorene::Scalar::set_grid_point
double & set_grid_point(int l, int k, int j, int i)
Setting the value of the field at a given grid point.
Definition scalar.h:684
Lorene::Scalar::va
Valeur va
The numerical value of the Scalar.
Definition scalar.h:405
Lorene::Scalar::asymptot
Valeur ** asymptot(int n, const int flag=0) const
Asymptotic expansion at r = infinity.
Definition scalar_asymptot.C:63
Lorene::Scalar::dzpuis
int dzpuis
Power of r by which the quantity represented by this must be divided in the compactified external d...
Definition scalar.h:403
Lorene::Tbl::set_etat_zero
void set_etat_zero()
Sets the logical state to ETATZERO (zero).
Definition tbl.C:347
Lorene::Tbl::set_etat_qcq
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
Definition tbl.C:361
Lorene::Valeur
Values and coefficients of a (real-value) function.
Definition valeur.h:287
Lorene::Valeur::set_etat_c_qcq
void set_etat_c_qcq()
Sets the logical state to ETATQCQ (ordinary state) for values in the configuration space (Mtbl c ).
Definition valeur.C:701
Lorene::Valeur::set_etat_zero
void set_etat_zero()
Sets the logical state to ETATZERO (zero).
Definition valeur.C:689
Lorene::Valeur::set_base
void set_base(const Base_val &)
Sets the bases for spectral expansions (member base ).
Definition valeur.C:810
Lorene::Valeur::c
Mtbl * c
Values of the function at the points of the multi-grid.
Definition valeur.h:299
Lorene::Valeur::set
Tbl & set(int l)
Read/write of the value in a given domain (configuration space).
Definition valeur.h:363
Lorene::Valeur::base
Base_val base
Bases on which the spectral expansion is performed.
Definition valeur.h:305
Lorene::Tensor::mp
const Map *const mp
Mapping on which the numerical values at the grid points are defined.
Definition tensor.h:295
Lorene
Lorene prototypes.
Definition app_hor.h:64