ope_helmholtz_plus.C
/*
* Copyright (c) 2003 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 ope_helmholtz_plus_C[] = "$Header: /cvsroot/Lorene/C++/Source/Ope_elementary/ope_helmholtz_plus.C,v 1.7 2014/10/13 08:53:33 j_novak Exp $" ;
/*
* $Id: ope_helmholtz_plus.C,v 1.7 2014/10/13 08:53:33 j_novak Exp $
* $Log: ope_helmholtz_plus.C,v $
* Revision 1.7 2014/10/13 08:53:33 j_novak
* Lorene classes and functions now belong to the namespace Lorene.
*
* Revision 1.6 2014/10/06 15:13:15 j_novak
* Modified #include directives to use c++ syntax.
*
* Revision 1.5 2007/05/06 10:48:13 p_grandclement
* Modification of a few operators for the vorton project
*
* Revision 1.4 2004/01/15 09:15:38 p_grandclement
* Modification and addition of the Helmholtz operators
*
* Revision 1.3 2003/12/11 16:12:10 e_gourgoulhon
* Changed sqrt(2) to sqrt(double(2)).
*
* Revision 1.2 2003/12/11 15:57:27 p_grandclement
* include stdlib.h encore ...
*
* Revision 1.1 2003/12/11 14:48:50 p_grandclement
* Addition of ALL (and that is a lot !) the files needed for the general elliptic solver ... UNDER DEVELOPEMENT...
*
*
* $Header: /cvsroot/Lorene/C++/Source/Ope_elementary/ope_helmholtz_plus.C,v 1.7 2014/10/13 08:53:33 j_novak Exp $
*
*/
#include <cmath>
#include <cstdlib>
#include "proto.h"
#include "ope_elementary.h"
// Standard constructor :
namespace Lorene {
Ope_helmholtz_plus::Ope_helmholtz_plus (int nbr, int base, int lquant, double alf,
double bet, double mas):
Ope_elementary(nbr, base, alf, bet), lq(lquant), masse (mas) {
}
// Constructor by copy :
Ope_helmholtz_plus::Ope_helmholtz_plus (const Ope_helmholtz_plus& so) :
Ope_elementary(so), lq (so.lq), masse (so.masse) {
}
// Destructor :
Ope_helmholtz_plus::~Ope_helmholtz_plus() {}
// True functions :
void Ope_helmholtz_plus::do_ope_mat() const {
if (ope_mat != 0x0)
delete ope_mat ;
ope_mat = new Matrice
(helmholtz_plus_mat(nr, lq, alpha, beta, masse, base_r)) ;
}
void Ope_helmholtz_plus::do_ope_cl() const {
if (ope_mat == 0x0)
do_ope_mat() ;
if (ope_cl != 0x0)
delete ope_cl ;
ope_cl = new Matrice
(cl_helmholtz_plus(*ope_mat, base_r)) ;
}
void Ope_helmholtz_plus::do_non_dege() const {
if (ope_cl == 0x0)
do_ope_cl() ;
if (non_dege != 0x0)
delete non_dege ;
non_dege = new Matrice
(prepa_helmholtz_plus_nondege(*ope_cl, base_r)) ;
}
Tbl Ope_helmholtz_plus::get_solp (const Tbl& so) const {
if (non_dege == 0x0)
do_non_dege() ;
Tbl res(solp_helmholtz_plus (*ope_mat, *non_dege, so, alpha, beta, base_r));
Tbl valeurs (val_solp (res, alpha, base_r)) ;
sp_plus = valeurs(0) ;
sp_minus = valeurs(1) ;
dsp_plus = valeurs(2) ;
dsp_minus = valeurs(3) ;
return res ;
}
Tbl Ope_helmholtz_plus::get_solh() const {
Tbl res (solh_helmholtz_plus (nr, lq, alpha, beta, masse, base_r)) ;
// Un peu tricky...
if (res.get_ndim() == 1) {
Tbl val_lim (val_solp (res, alpha, base_r)) ;
s_one_plus = val_lim(0) ;
s_one_minus = val_lim(1) ;
ds_one_plus = val_lim(2) ;
ds_one_minus = val_lim(3) ;
}
else {
Tbl auxi (nr) ;
auxi.set_etat_qcq() ;
for (int i=0 ; i<nr ; i++)
auxi.set(i) = res(0,i) ;
Tbl val_one (val_solp (auxi, alpha, base_r)) ;
s_one_plus = val_one(0) ;
s_one_minus = val_one(1) ;
ds_one_plus = val_one(2) ;
ds_one_minus = val_one(3) ;
for (int i=0 ; i<nr ; i++)
auxi.set(i) = res(1,i) ;
Tbl val_two (val_solp (auxi, alpha, base_r)) ;
s_two_plus = val_two(0) ;
s_two_minus = val_two(1) ;
ds_two_plus = val_two(2) ;
ds_two_minus = val_two(3) ;
}
return res ;
}
void Ope_helmholtz_plus::inc_l_quant() {
cout << "inc_l_quant not implemented for Helmholtz operator." << endl ;
abort() ;
}
}
Lorene::Matrice
Matrix handling.
Definition matrice.h:152
Lorene::Ope_elementary::ope_mat
Matrice * ope_mat
Pointer on the matrix representation of the operator.
Definition ope_elementary.h:114
Lorene::Ope_elementary::ds_two_minus
double ds_two_minus
Value of the derivative of the second homogeneous solution at the inner boundary.
Definition ope_elementary.h:160
Lorene::Ope_elementary::s_two_plus
double s_two_plus
Value of the second homogeneous solution at the outer boundary.
Definition ope_elementary.h:146
Lorene::Ope_elementary::s_one_minus
double s_one_minus
Value of the first homogeneous solution at the inner boundary.
Definition ope_elementary.h:131
Lorene::Ope_elementary::beta
double beta
Parameter of the associated mapping.
Definition ope_elementary.h:109
Lorene::Ope_elementary::dsp_plus
double dsp_plus
Value of the derivative of the particular solution at the outer boundary.
Definition ope_elementary.h:177
Lorene::Ope_elementary::ds_one_plus
double ds_one_plus
Value of the derivative of the first homogeneous solution at the outer boundary.
Definition ope_elementary.h:136
Lorene::Ope_elementary::ds_one_minus
double ds_one_minus
Value of the derivative of the first homogeneous solution at the inner boundary.
Definition ope_elementary.h:141
Lorene::Ope_elementary::Ope_elementary
Ope_elementary(int nbr, int baser, double alf, double eta)
Standard constructor, protected because the class is an abstract one.
Definition ope_elementary.C:42
Lorene::Ope_elementary::alpha
double alpha
Parameter of the associated mapping.
Definition ope_elementary.h:108
Lorene::Ope_elementary::s_two_minus
double s_two_minus
Value of the second homogeneous solution at the inner boundary.
Definition ope_elementary.h:150
Lorene::Ope_elementary::base_r
int base_r
Radial basis of decomposition.
Definition ope_elementary.h:107
Lorene::Ope_elementary::s_one_plus
double s_one_plus
Value of the first homogeneous solution at the outer boundary.
Definition ope_elementary.h:127
Lorene::Ope_elementary::ope_cl
Matrice * ope_cl
Pointer on the banded-matrix of the operator.
Definition ope_elementary.h:118
Lorene::Ope_elementary::sp_minus
double sp_minus
Value of the particular solution at the inner boundary.
Definition ope_elementary.h:165
Lorene::Ope_elementary::ds_two_plus
double ds_two_plus
Value of the derivative of the second homogeneous solution at the outer boundary.
Definition ope_elementary.h:155
Lorene::Ope_elementary::non_dege
Matrice * non_dege
Pointer on the non-degenerated matrix of the operator.
Definition ope_elementary.h:122
Lorene::Ope_elementary::sp_plus
double sp_plus
Value of the particular solution at the outer boundary.
Definition ope_elementary.h:169
Lorene::Ope_elementary::nr
int nr
Number of radial points.
Definition ope_elementary.h:106
Lorene::Ope_elementary::dsp_minus
double dsp_minus
Value of the derivative of the particular solution at the inner boundary.
Definition ope_elementary.h:173
Lorene::Ope_helmholtz_plus::lq
int lq
The quantum number l.
Definition ope_elementary.h:449
Lorene::Ope_helmholtz_plus::do_non_dege
virtual void do_non_dege() const
Computes the non-degenerated matrix of the operator.
Definition ope_helmholtz_plus.C:92
Lorene::Ope_helmholtz_plus::get_solh
virtual Tbl get_solh() const
Computes the homogeneous solutions(s).
Definition ope_helmholtz_plus.C:119
Lorene::Ope_helmholtz_plus::~Ope_helmholtz_plus
virtual ~Ope_helmholtz_plus()
Destructor.
Definition ope_helmholtz_plus.C:70
Lorene::Ope_helmholtz_plus::do_ope_cl
virtual void do_ope_cl() const
Computes the banded-matrix of the operator.
Definition ope_helmholtz_plus.C:81
Lorene::Ope_helmholtz_plus::Ope_helmholtz_plus
Ope_helmholtz_plus(int nbr, int baser, int lq, double alf, double bet, double mas)
Standard constructor.
Definition ope_helmholtz_plus.C:59
Lorene::Ope_helmholtz_plus::masse
double masse
The mass parameter m .
Definition ope_elementary.h:450
Lorene::Ope_helmholtz_plus::do_ope_mat
virtual void do_ope_mat() const
Computes the matrix of the operator.
Definition ope_helmholtz_plus.C:73
Lorene::Ope_helmholtz_plus::get_solp
virtual Tbl get_solp(const Tbl &so) const
Computes the particular solution, given the source so .
Definition ope_helmholtz_plus.C:103
Lorene::Ope_helmholtz_plus::inc_l_quant
virtual void inc_l_quant()
Increases the quatum number l by one unit (CURRENTLY NOT IMPLEMENTED).
Definition ope_helmholtz_plus.C:160
Lorene::Tbl
Basic array class.
Definition tbl.h:161
Lorene::Tbl::get_ndim
int get_ndim() const
Gives the number of dimensions (ie dim.ndim).
Definition tbl.h:400
Lorene::Tbl::set_etat_qcq
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
Definition tbl.C:361
Lorene::Tbl::set
double & set(int i)
Read/write of a particular element (index i) (1D case).
Definition tbl.h:281
Lorene
Lorene prototypes.
Definition app_hor.h:64