21char sol_elliptic_only_zec_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/PDE/sol_elliptic_only_zec.C,v 1.4 2014/10/13 08:53:30 j_novak Exp $" ;
56#include "param_elliptic.h"
66Mtbl_cf elliptic_solver_only_zec (
const Param_elliptic& ope_var,
const Mtbl_cf& source,
double val) {
70 assert (source.get_mg()->get_type_r(nz-1) == UNSURR) ;
82 Mtbl_cf solution_part(source.get_mg(), source.base) ;
83 Mtbl_cf solution_hom_un(source.get_mg(), source.base) ;
84 Mtbl_cf solution_hom_deux(source.get_mg(), source.base) ;
85 Mtbl_cf resultat(source.get_mg(), source.base) ;
87 solution_part.annule_hard() ;
88 solution_hom_un.annule_hard() ;
89 solution_hom_deux.annule_hard() ;
90 resultat.annule_hard() ;
94 for (
int l=0 ; l<nz-1 ; l++) {
95 nt = source.get_mg()->get_nt(l) ;
96 np = source.get_mg()->get_np(l) ;
97 for (
int k=0 ; k<np+1 ; k++)
98 for (
int j=0 ; j<nt ; j++)
101 int conte = conte_start ;
105 nr = source.get_mg()->get_nr(zone) ;
106 nt = source.get_mg()->get_nt(zone) ;
107 np = source.get_mg()->get_np(zone) ;
109 for (
int k=0 ; k<np+1 ; k++)
110 for (
int j=0 ; j<nt ; j++) {
112 if (ope_var.operateurs[conte] != 0x0) {
115 sol_hom =
new Tbl(ope_var.operateurs[conte]->get_solh()) ;
120 for (
int i=0 ; i<nr ; i++)
121 so->set(i) = source(zone, k, j, i) ;
123 sol_part =
new Tbl(ope_var.operateurs[conte]->get_solp(*so)) ;
126 for (
int i=0 ; i<nr ; i++) {
127 solution_part.set(zone, k, j, i) = (*sol_part)(i) ;
128 if (sol_hom->get_ndim()==1)
129 solution_hom_un.set(zone, k, j, i) = (*sol_hom)(i) ;
132 solution_hom_un.set(zone, k, j, i) = (*sol_hom)(0,i) ;
133 solution_hom_deux.set(zone, k, j, i) = (*sol_hom)(1,i) ;
153 int start = conte_start ;
154 for (
int k=0 ; k<1 ; k++)
155 for (
int j=0 ; j<1 ; j++) {
156 if (ope_var.operateurs[start] != 0x0) {
160 double facteur = ((val-ope_var.F_minus(nz-1, k, j))/
161 ope_var.G_minus(nz-1) -
162 ope_var.operateurs[start]->val_sp_minus()) /
163 ope_var.operateurs[start]->val_sh_one_minus() ;
166 nr = source.get_mg()->get_nr(nz-1) ;
167 for (
int i=0 ; i<nr ; i++)
168 resultat.set(nz-1,k,j,i) = solution_part(nz-1,k,j,i) +
169 facteur*solution_hom_un(nz-1,k,j,i) ;
int get_nzone() const
Returns the number of domains.
const Mg3d * get_mg() const
Returns the Mg3d on which the Mtbl_cf is defined.