param_elliptic_2d.C
/*
* Copyright (c) 2004 Philippe Grandclement
*
* 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 version 2
* as published by the Free Software Foundation.
*
* 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 param_elliptic_2d_C[] = "$Header: /cvsroot/Lorene/C++/Source/Param_elliptic/param_elliptic_2d.C,v 1.3 2014/10/13 08:53:37 j_novak Exp $" ;
/*
* $Id: param_elliptic_2d.C,v 1.3 2014/10/13 08:53:37 j_novak Exp $
* $Log: param_elliptic_2d.C,v $
* Revision 1.3 2014/10/13 08:53:37 j_novak
* Lorene classes and functions now belong to the namespace Lorene.
*
* Revision 1.2 2014/10/06 15:13:15 j_novak
* Modified #include directives to use c++ syntax.
*
* Revision 1.1 2004/08/24 09:14:49 p_grandclement
* Addition of some new operators, like Poisson in 2d... It now requieres the
* GSL library to work.
*
* Also, the way a variable change is stored by a Param_elliptic is changed and
* no longer uses Change_var but rather 2 Scalars. The codes using that feature
* will requiere some modification. (It should concern only the ones about monopoles)
*
*
* $Header: /cvsroot/Lorene/C++/Source/Param_elliptic/param_elliptic_2d.C,v 1.3 2014/10/13 08:53:37 j_novak Exp $
*
*/
#include "headcpp.h"
#include <cmath>
#include <cstdlib>
#include "base_val.h"
#include "map.h"
#include "ope_elementary.h"
#include "param_elliptic.h"
#include "change_var.h"
#include "scalar.h"
namespace Lorene {
void Param_elliptic::set_poisson_2d(const Scalar& source, bool indic) {
int dzpuis = source.get_dzpuis() ;
if (type_map != MAP_AFF) {
cout << "set_poisson_2d only defined for an affine mapping..." << endl ;
abort() ;
}
else {
int nz = get_mp().get_mg()->get_nzone() ;
int nr ;
double alpha, beta ;
int m_quant, l_quant, base_r_1d ;
int conte = 0 ;
for (int l=0 ; l<nz ; l++) {
nr = get_mp().get_mg()->get_nr(l) ;
alpha = get_alpha (l) ;
beta = get_beta (l) ;
for (int k=0 ; k<get_mp().get_mg()->get_np(l)+1 ; k++)
for (int j=0 ; j<get_mp().get_mg()->get_nt(l) ; j++) {
if (operateurs[conte] != 0x0)
delete operateurs[conte] ;
source.get_spectral_va().base.give_quant_numbers(l, k, j, m_quant, l_quant, base_r_1d) ;
if (k!=1)
if ((indic) || ((!indic) && (l_quant !=0)))
operateurs[conte] = new Ope_poisson_2d (nr, base_r_1d, alpha, beta, l_quant, dzpuis) ;
else
operateurs[conte] = 0x0 ;
else
operateurs[conte] = 0x0 ;
conte ++ ;
}
}
}
}
void Param_elliptic::set_helmholtz_minus_2d(int zone, double masse, const Scalar& source) {
int dzpuis = source.get_dzpuis() ;
assert (masse > 0) ;
if (type_map != MAP_AFF) {
cout << "set_helmholtz_minus_2d only defined for an affine mapping..." << endl ;
abort() ;
}
else {
int nz = get_mp().get_mg()->get_nzone() ;
if (zone == nz-1)
source.check_dzpuis(2) ;
int nr ;
double alpha, beta ;
int m_quant, l_quant, base_r_1d ;
int conte = 0 ;
for (int l=0 ; l<nz ; l++) {
nr = get_mp().get_mg()->get_nr(l) ;
alpha = get_alpha (l) ;
beta = get_beta (l) ;
for (int k=0 ; k<get_mp().get_mg()->get_np(l)+1 ; k++)
for (int j=0 ; j<get_mp().get_mg()->get_nt(l) ; j++) {
if (l==zone) {
if (operateurs[conte] != 0x0) {
delete operateurs[conte] ;
source.get_spectral_va().base.give_quant_numbers
(l, k, j, m_quant, l_quant, base_r_1d) ;
operateurs[conte] = new Ope_helmholtz_minus_2d (nr, base_r_1d, alpha, beta, l_quant, masse, dzpuis) ;
}
}
conte ++ ;
}
}
}
}
}
Lorene::Base_val::give_quant_numbers
void give_quant_numbers(int, int, int, int &, int &, int &) const
Computes the various quantum numbers and 1d radial base.
Definition base_val_quantum.C:81
Lorene::Ope_helmholtz_minus_2d
Class for the operator of the Helmholtz equation in 2D.
Definition ope_elementary.h:789
Lorene::Ope_poisson_2d
Class for the operator of the Poisson equation in 2D.
Definition ope_elementary.h:723
Lorene::Param_elliptic::get_mp
const Map_radial & get_mp() const
Returns the mapping.
Definition param_elliptic.C:118
Lorene::Param_elliptic::operateurs
Ope_elementary ** operateurs
Array on the elementary operators.
Definition param_elliptic.h:142
Lorene::Param_elliptic::set_poisson_2d
void set_poisson_2d(const Scalar &, bool indic=false)
Set everything to do a 2d-Poisson, with or without l=0 (not put by default...).
Definition param_elliptic_2d.C:59
Lorene::Param_elliptic::type_map
int type_map
Type of mapping either MAP_AFF or MAP_LOG.
Definition param_elliptic.h:138
Lorene::Param_elliptic::set_helmholtz_minus_2d
void set_helmholtz_minus_2d(int zone, double mas, const Scalar &)
Set the 2D Helmholtz operator (with minus sign).
Definition param_elliptic_2d.C:100
Lorene::Scalar
Tensor field of valence 0 (or component of a tensorial field).
Definition scalar.h:387
Lorene::Scalar::get_dzpuis
int get_dzpuis() const
Returns dzpuis.
Definition scalar.h:557
Lorene::Scalar::check_dzpuis
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...
Definition scalar.C:873
Lorene::Scalar::get_spectral_va
const Valeur & get_spectral_va() const
Returns va (read only version).
Definition scalar.h:601
Lorene::Valeur::base
Base_val base
Bases on which the spectral expansion is performed.
Definition valeur.h:305
Lorene
Lorene prototypes.
Definition app_hor.h:64