21 char param_elliptic_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Param_elliptic/param_elliptic.C,v 1.21 2014/10/13 08:53:37 j_novak Exp $" ;
110 #include "ope_elementary.h"
111 #include "param_elliptic.h"
112 #include "base_val.h"
128 cout <<
"Unknown mapping in Param_elliptic" << endl ;
134 double Param_elliptic::get_alpha(
int l)
const {
144 cout <<
"Unknown mapping in Param_elliptic" << endl ;
150 double Param_elliptic::get_beta(
int l)
const {
160 cout <<
"Unknown mapping in Param_elliptic" << endl ;
166 int Param_elliptic::get_type(
int l)
const {
176 cout <<
"Unknown mapping in Param_elliptic" << endl ;
184 done_F (so.get_mp().get_mg()->get_nzone(),
185 so.get_mp().get_mg()->get_np(0) + 1,
186 so.get_mp().get_mg()->get_nt(0)),
187 done_G (so.get_mp().get_mg()->get_nzone()),
188 val_F_plus (so.get_mp().get_mg()->get_nzone(),
189 so.get_mp().get_mg()->get_np(0) + 1,
190 so.get_mp().get_mg()->get_nt(0)),
191 val_F_minus (so.get_mp().get_mg()->get_nzone(),
192 so.get_mp().get_mg()->get_np(0) + 1,
193 so.get_mp().get_mg()->get_nt(0)),
194 val_dF_plus (so.get_mp().get_mg()->get_nzone(),
195 so.get_mp().get_mg()->get_np(0) + 1,
196 so.get_mp().get_mg()->get_nt(0)),
197 val_dF_minus (so.get_mp().get_mg()->get_nzone(),
198 so.get_mp().get_mg()->get_np(0) + 1,
199 so.get_mp().get_mg()->get_nt(0)),
200 val_G_plus (so.get_mp().get_mg()->get_nzone()),
201 val_G_minus (so.get_mp().get_mg()->get_nzone()),
202 val_dG_plus (so.get_mp().get_mg()->get_nzone()),
203 val_dG_minus (so.get_mp().get_mg()->get_nzone())
221 if ((map_affine == 0x0) && (map_log == 0x0)) {
222 cout <<
"Param_elliptic not yet defined on this type of mapping" << endl ;
227 if (map_affine != 0x0) {
232 if (map_log != 0x0) {
239 for (
int l=0 ; l<nz ; l++)
240 nbr += (
get_mp().get_mg()->get_np(l)+1)*
245 for (
int l=0 ; l<nbr ; l++)
249 int base_r, m_quant, l_quant ;
252 for (
int l=0 ; l<nz ; l++) {
260 (l, k, j, m_quant, l_quant, base_r) ;
266 get_beta(l), l_quant, dzpuis) ;
272 get_beta(l), l_quant, dzpuis) ;
279 cout <<
"Unknown mapping in Param_elliptic::Param_elliptic"
316 for (
int i=0 ; i<nbr ; i++)
333 for (
int k=0 ; k<np+1 ; k++)
334 for (
int j=0 ; j<nt ; j++) {
351 int m_quant, base_r_1d, l_quant ;
357 for (
int l=0 ; l<nz ; l++) {
367 (l, k, j, m_quant, l_quant, base_r_1d) ;
369 get_beta(l) , masse) ;
379 cout <<
"Operator not define with LOG mapping..." << endl ;
383 for (
int l=0 ; l<nz ; l++) {
393 (l, k, j, m_quant, l_quant, base_r_1d) ;
395 get_alpha(l), get_beta(l), masse) ;
405 cout <<
"Unkown mapping in set_helmhotz_minus" << endl ;
419 int m_quant, base_r_1d, l_quant ;
425 for (
int l=0 ; l<nz ; l++) {
438 (l, k, j, m_quant, l_quant, base_r_1d) ;
440 get_beta(l) , masse) ;
451 cout <<
"Operator not define with LOG mapping..." << endl ;
455 for (
int l=0 ; l<nz ; l++) {
468 (l, k, j, m_quant, l_quant, base_r_1d) ;
470 get_alpha(l), get_beta(l), masse) ;
481 cout <<
"Unkown mapping in set_helmhotz_minus" << endl ;
490 cout <<
"set_poisson_vect_r only defined for an affine mapping..." << endl ;
500 for (
int l=0 ; l<nz ; l++) {
506 if ((
operateurs[conte] != 0x0) && (l==zone)) {
510 assert (op_pois !=0x0) ;
516 if ((!only_l_zero)||(lq_old == 0)) {
518 get_beta(l), lq_old, dzp) ;
538 for (
int l=0 ; l<nz ; l++) {
543 if ((
operateurs[conte] != 0x0) && (l==zone)) {
546 assert (op_pois !=0x0) ;
563 cout <<
"set_poisson_tens_rr only defined for an affine mapping..."
574 for (
int l=0 ; l<nz ; l++) {
580 if ((
operateurs[conte] != 0x0) && (l==zone)) {
584 assert (op_pois !=0x0) ;
591 (nr, old_base,get_alpha(l), get_beta(l), lq_old, dzp) ;
613 for (
int l=0 ; l<nz ; l++) {
619 if ((
operateurs[conte] != 0x0) && (l==zone)) {
626 get_beta(l), a, b, c) ;
636 cout <<
"Operator not define with LOG mapping..." << endl ;
640 for (
int l=0 ; l<nz ; l++) {
646 if ((
operateurs[conte] != 0x0) && (l==zone)) {
653 get_beta(l), a, b, c) ;
663 cout <<
"Unkown mapping in set_sec_order_r2" << endl ;
672 cout <<
"set_sec_order only defined for a log mapping" << endl ;
682 for (
int l=0 ; l<nz ; l++) {
688 if ((
operateurs[conte] != 0x0) && (l==zone)) {
696 get_beta(l), a, b, c) ;
713 int m_quant, base_r_1d, l_quant ;
720 for (
int l=0 ; l<nz ; l++) {
721 int dz = (l==nz-1) ? dzpuis : 0 ;
730 (l, k, j, m_quant, l_quant, base_r_1d) ;
732 get_beta(l), l_quant, dz) ;
742 cout <<
"Operator not define with LOG mapping..." << endl ;
746 for (
int l=0 ; l<nz ; l++) {
747 int dz = (l==nz-1) ? dzpuis : 0 ;
756 (l, k, j, m_quant, l_quant, base_r_1d) ;
758 get_alpha(l), get_beta(l), l_quant, dz) ;
768 cout <<
"Unkown mapping in set_ope_vorton" << endl ;
776 assert (so.
get_etat() != ETATNONDEF) ;
785 assert (so.
get_etat() != ETATNONDEF) ;