LORENE
ope_helmholtz_minus.C
1 /*
2  * Copyright (c) 2003 Philippe Grandclement
3  *
4  * This file is part of LORENE.
5  *
6  * LORENE is free software; you can redistribute it and/or modify
7  * it under the terms of the GNU General Public License version 2
8  * as published by the Free Software Foundation.
9  *
10  * LORENE is distributed in the hope that it will be useful,
11  * but WITHOUT ANY WARRANTY; without even the implied warranty of
12  * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
13  * GNU General Public License for more details.
14  *
15  * You should have received a copy of the GNU General Public License
16  * along with LORENE; if not, write to the Free Software
17  * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA 02111-1307 USA
18  *
19  */
20 
21 char ope_helmholtz_minus_C[] = "$Header: /cvsroot/Lorene/C++/Source/Ope_elementary/ope_helmholtz_minus.C,v 1.8 2014/10/13 08:53:32 j_novak Exp $" ;
22 
23 /*
24  * $Id: ope_helmholtz_minus.C,v 1.8 2014/10/13 08:53:32 j_novak Exp $
25  * $Log: ope_helmholtz_minus.C,v $
26  * Revision 1.8 2014/10/13 08:53:32 j_novak
27  * Lorene classes and functions now belong to the namespace Lorene.
28  *
29  * Revision 1.7 2014/10/06 15:13:15 j_novak
30  * Modified #include directives to use c++ syntax.
31  *
32  * Revision 1.6 2008/07/10 10:59:17 p_grandclement
33  * forgot another ones
34  *
35  * Revision 1.5 2004/08/24 09:14:45 p_grandclement
36  * Addition of some new operators, like Poisson in 2d... It now requieres the
37  * GSL library to work.
38  *
39  * Also, the way a variable change is stored by a Param_elliptic is changed and
40  * no longer uses Change_var but rather 2 Scalars. The codes using that feature
41  * will requiere some modification. (It should concern only the ones about monopoles)
42  *
43  * Revision 1.4 2004/01/15 09:15:38 p_grandclement
44  * Modification and addition of the Helmholtz operators
45  *
46  * Revision 1.3 2003/12/11 16:12:10 e_gourgoulhon
47  * Changed sqrt(2) to sqrt(double(2)).
48  *
49  * Revision 1.2 2003/12/11 15:57:27 p_grandclement
50  * include stdlib.h encore ...
51  *
52  * Revision 1.1 2003/12/11 14:48:50 p_grandclement
53  * Addition of ALL (and that is a lot !) the files needed for the general elliptic solver ... UNDER DEVELOPEMENT...
54  *
55  *
56  * $Header: /cvsroot/Lorene/C++/Source/Ope_elementary/ope_helmholtz_minus.C,v 1.8 2014/10/13 08:53:32 j_novak Exp $
57  *
58  */
59 #include <cmath>
60 #include <cstdlib>
61 
62 #include "proto.h"
63 #include "ope_elementary.h"
64 
65 // Standard constructor :
66 namespace Lorene {
67 Ope_helmholtz_minus::Ope_helmholtz_minus (int nbr, int base, int lquant,
68  double alf, double bet, double mas):
69  Ope_elementary(nbr, base, alf, bet), lq (lquant), masse (mas) {
70 }
71 
72 // Constructor by copy :
74  Ope_elementary(so), lq(so.lq), masse (so.masse) {
75 }
76 
77 // Destructor :
79 
80 // True functions :
82  if (ope_mat != 0x0)
83  delete ope_mat ;
84 
85  ope_mat = new Matrice
86  (helmholtz_minus_mat(nr, lq, alpha, beta, masse, base_r)) ;
87 }
88 
90  if (ope_mat == 0x0)
91  do_ope_mat() ;
92 
93  if (ope_cl != 0x0)
94  delete ope_cl ;
95 
96  ope_cl = new Matrice
97  (cl_helmholtz_minus(*ope_mat, base_r)) ;
98 }
99 
101  if (ope_cl == 0x0)
102  do_ope_cl() ;
103 
104  if (non_dege != 0x0)
105  delete non_dege ;
106 
107  non_dege = new Matrice
108  (prepa_helmholtz_minus_nondege(*ope_cl, base_r)) ;
109 }
110 
112 
113  if (non_dege == 0x0)
114  do_non_dege() ;
115 
116  Tbl res(solp_helmholtz_minus (*ope_mat, *non_dege, so, alpha, beta, lq, base_r));
117 
118  Tbl valeurs (val_solp (res, alpha, base_r)) ;
119  sp_plus = valeurs(0) ;
120  sp_minus = valeurs(1) ;
121  dsp_plus = valeurs(2) ;
122  dsp_minus = valeurs(3) ;
123 
124  return res ;
125 }
126 
128 
129  Tbl res (solh_helmholtz_minus (nr, lq, alpha, beta, masse, base_r)) ;
130 
131  // Un peu tricky...
132  if (res.get_ndim() == 1) {
133  Tbl val_lim (val_solp (res, alpha, base_r)) ;
134 
135  s_one_plus = val_lim(0) ;
136  s_one_minus = val_lim(1) ;
137  ds_one_plus = val_lim(2) ;
138  ds_one_minus = val_lim(3) ;
139 
140  }
141  else {
142  Tbl auxi (nr) ;
143  auxi.set_etat_qcq() ;
144  for (int i=0 ; i<nr ; i++)
145  auxi.set(i) = res(0,i) ;
146 
147  Tbl val_one (val_solp (auxi, alpha, base_r)) ;
148 
149  s_one_plus = val_one(0) ;
150  s_one_minus = val_one(1) ;
151  ds_one_plus = val_one(2) ;
152  ds_one_minus = val_one(3) ;
153 
154  for (int i=0 ; i<nr ; i++)
155  auxi.set(i) = res(1,i) ;
156 
157  Tbl val_two (val_solp (auxi, alpha, base_r)) ;
158 
159  s_two_plus = val_two(0) ;
160  s_two_minus = val_two(1) ;
161  ds_two_plus = val_two(2) ;
162  ds_two_minus = val_two(3) ;
163  }
164  return res ;
165 }
166 
167 
168 
170 
171  cout << "inc_l_quant not implemented for Helmholtz operator." << endl ;
172  abort() ;
173 }
174 }
Lorene::Ope_elementary::ope_cl
Matrice * ope_cl
Pointer on the banded-matrix of the operator.
Definition: ope_elementary.h:118
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::Tbl::set
double & set(int i)
Read/write of a particular element (index i) (1D case)
Definition: tbl.h:281
Lorene::Ope_helmholtz_minus::masse
double masse
The mass parameter m .
Definition: ope_elementary.h:394
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_helmholtz_minus::get_solh
virtual Tbl get_solh() const
Computes the homogeneous solutions(s).
Definition: ope_helmholtz_minus.C:127
Lorene
Lorene prototypes.
Definition: app_hor.h:64
Lorene::Ope_elementary::base_r
int base_r
Radial basis of decomposition.
Definition: ope_elementary.h:107
Lorene::Ope_elementary::nr
int nr
Number of radial points.
Definition: ope_elementary.h:106
Lorene::Tbl::get_ndim
int get_ndim() const
Gives the number of dimensions (ie dim.ndim)
Definition: tbl.h:400
Lorene::Ope_helmholtz_minus::inc_l_quant
virtual void inc_l_quant()
Increases the quatum number l by one unit (CURRENTLY NOT IMPLEMENTED)
Definition: ope_helmholtz_minus.C:169
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_helmholtz_minus
Class for the Helmholtz operator ( ).
Definition: ope_elementary.h:390
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::beta
double beta
Parameter of the associated mapping.
Definition: ope_elementary.h:109
Lorene::Tbl::set_etat_qcq
void set_etat_qcq()
Sets the logical state to ETATQCQ (ordinary state).
Definition: tbl.C:361
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::sp_minus
double sp_minus
Value of the particular solution at the inner boundary.
Definition: ope_elementary.h:165
Lorene::Ope_helmholtz_minus::do_ope_mat
virtual void do_ope_mat() const
Computes the matrix of the operator.
Definition: ope_helmholtz_minus.C:81
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_helmholtz_minus::Ope_helmholtz_minus
Ope_helmholtz_minus(int nbr, int baser, int lq, double alf, double bet, double mas)
Standard constructor.
Definition: ope_helmholtz_minus.C:67
Lorene::Tbl
Basic array class.
Definition: tbl.h:161
Lorene::Ope_helmholtz_minus::get_solp
virtual Tbl get_solp(const Tbl &so) const
Computes the particular solution, given the source so .
Definition: ope_helmholtz_minus.C:111
Lorene::Ope_helmholtz_minus::lq
int lq
The quantum number l.
Definition: ope_elementary.h:393
Lorene::Ope_helmholtz_minus::do_ope_cl
virtual void do_ope_cl() const
Computes the banded-matrix of the operator.
Definition: ope_helmholtz_minus.C:89
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::sp_plus
double sp_plus
Value of the particular solution at the outer boundary.
Definition: ope_elementary.h:169
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_helmholtz_minus::~Ope_helmholtz_minus
virtual ~Ope_helmholtz_minus()
Destructor.
Definition: ope_helmholtz_minus.C:78
Lorene::Matrice
Matrix handling.
Definition: matrice.h:152
Lorene::Ope_elementary
Basic class for elementary elliptic operators.
Definition: ope_elementary.h:102
Lorene::Ope_helmholtz_minus::do_non_dege
virtual void do_non_dege() const
Computes the non-degenerated matrix of the operator.
Definition: ope_helmholtz_minus.C:100
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::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::dsp_minus
double dsp_minus
Value of the derivative of the particular solution at the inner boundary.
Definition: ope_elementary.h:173
Lorene::Ope_elementary::ope_mat
Matrice * ope_mat
Pointer on the matrix representation of the operator.
Definition: ope_elementary.h:114
Lorene::Ope_elementary::alpha
double alpha
Parameter of the associated mapping.
Definition: ope_elementary.h:108