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
21char 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 :
66namespace Lorene {
68 double alf, double bet, double mas):
69 Ope_elementary(nbr, base, alf, bet), lq (lquant), masse (mas) {
70}
71
72// Constructor by copy :
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}
Time evolution with partial storage (*** under development ***).
Definition evolution.h:371
Matrix handling.
Definition matrice.h:152
Basic class for elementary elliptic operators.
Matrice * ope_mat
Pointer on the matrix representation of the operator.
double ds_two_minus
Value of the derivative of the second homogeneous solution at the inner boundary.
double s_two_plus
Value of the second homogeneous solution at the outer boundary.
double s_one_minus
Value of the first homogeneous solution at the inner boundary.
double beta
Parameter of the associated mapping.
double dsp_plus
Value of the derivative of the particular solution at the outer boundary.
double ds_one_plus
Value of the derivative of the first homogeneous solution at the outer boundary.
double ds_one_minus
Value of the derivative of the first homogeneous solution at the inner boundary.
double alpha
Parameter of the associated mapping.
double s_two_minus
Value of the second homogeneous solution at the inner boundary.
int base_r
Radial basis of decomposition.
double s_one_plus
Value of the first homogeneous solution at the outer boundary.
Matrice * ope_cl
Pointer on the banded-matrix of the operator.
double sp_minus
Value of the particular solution at the inner boundary.
double ds_two_plus
Value of the derivative of the second homogeneous solution at the outer boundary.
Matrice * non_dege
Pointer on the non-degenerated matrix of the operator.
double sp_plus
Value of the particular solution at the outer boundary.
int nr
Number of radial points.
double dsp_minus
Value of the derivative of the particular solution at the inner boundary.
Class for the Helmholtz operator ( ).
virtual void do_non_dege() const
Computes the non-degenerated matrix of the operator.
int lq
The quantum number l.
double masse
The mass parameter m .
Ope_helmholtz_minus(int nbr, int baser, int lq, double alf, double bet, double mas)
Standard constructor.
virtual void do_ope_cl() const
Computes the banded-matrix of the operator.
virtual Tbl get_solh() const
Computes the homogeneous solutions(s).
virtual Tbl get_solp(const Tbl &so) const
Computes the particular solution, given the source so .
virtual ~Ope_helmholtz_minus()
Destructor.
virtual void inc_l_quant()
Increases the quatum number l by one unit (CURRENTLY NOT IMPLEMENTED)
virtual void do_ope_mat() const
Computes the matrix of the operator.
Basic array class.
Definition tbl.h:161
Lorene prototypes.
Definition app_hor.h:64