75 if (
etat == ETATZERO) {
82 int nz =
mgrid.get_nzone() ;
93 if (
par_bc.get_n_int() > 2) {
102 int nt =
mgrid.get_nt(0) ;
103 int np =
mgrid.get_np(0) ;
132 if (
par_mat->get_n_matrice_mod() == 4)
145 int index = 0 ;
int index01 = 0 ;
146 int nr =
mgrid.get_nr(0);
151 system_odd.set(index, index) = -(2.*nr - 5.)/alpha ;
152 system_odd.set(index, index+1) = (2.*nr - 3.)/alpha ;
174 system_even.set(index+1, index-1) = 4*(nr-2)*(nr-2)/alpha ;
175 system_even.set(index+1, index) = 4*(nr-1)*(nr-1)/alpha ;
176 system_odd.set(index+1, index-1) = (2*nr-5)*(2*nr-5)/alpha ;
177 system_odd.set(index+1, index) = (2*nr-3)*(2*nr-3)/alpha ;
182 nr =
mgrid.get_nr(1) ;
185 alpha = -0.25/alpha ;
202 system_even.set(index, index) = -(nr-2)*(nr-2)/alpha ;
203 system_even.set(index, index+1) = (nr-1)*(nr-1)/alpha ;
204 system_odd.set(index, index) = -(nr-2)*(nr-2)/alpha ;
205 system_odd.set(index, index+1) = (nr-1)*(nr-1)/alpha ;
235 system_even.set(index+1, index-1) = (nr-2)*(nr-2)/alpha ;
236 system_even.set(index+1, index) = (nr-1)*(nr-1)/alpha ;
237 system_odd.set(index+1, index-1) = (nr-2)*(nr-2)/alpha ;
238 system_odd.set(index+1, index) = (nr-1)*(nr-1)/alpha ;
244 nr =
mgrid.get_nr(1) ;
247 alpha = -0.25/alpha ;
264 system_even.set(index+1, index) = (nr-1)*(nr-1)/alpha ;
265 system_odd.set(index+1, index) = (nr-1)*(nr-1)/alpha ;
274 alpha = -0.25/alpha ;
291 system_even.set(index, index) = -(nr-2)*(nr-2)/alpha ;
292 system_even.set(index, index+1) = (nr-1)*(nr-1)/alpha ;
293 system_odd.set(index, index) = -(nr-2)*(nr-2)/alpha ;
294 system_odd.set(index, index+1) = (nr-1)*(nr-1)/alpha ;
324 system_even.set(index+1, index-1) = (nr-2)*(nr-2)/alpha ;
325 system_even.set(index+1, index) = (nr-1)*(nr-1)/alpha ;
326 system_odd.set(index+1, index-1) = (nr-2)*(nr-2)/alpha ;
327 system_odd.set(index+1, index) = (nr-1)*(nr-1)/alpha ;
353 par_mat->get_matrice_mod(2) ) ;
355 par_mat->get_matrice_mod(3) ) ;
366 for (
int k=0;
k<np+2;
k++) {
367 for (
int j=0;
j<nt;
j++) {
369 if ((nullite_plm(
j, nt,
k, np, base) == 1)&&(
l_q >=
l_min)) {
372 int nl = (
l_q < 2 ? 1 : 2) ;
377 int nr =
mgrid.get_nr(0) ;
384 for (
int i=0;
i<nr-
nl;
i++) {
393 for (
int i=0;
i<nr-
nl-1;
i++) {
402 for (
int i=0;
i<nr-
nl;
i++) {
409 for (
int i=0;
i<nr-
nl-1;
i++) {
426 alpha = -0.25/alpha ;
427 nr =
mgrid.get_nr(1) ;
465 alpha = -0.25/alpha ;
466 nr =
mgrid.get_nr(1) ;
490 alpha = -0.25/alpha ;
540 nr =
mgrid.get_nr(0) ;
549 coef.
set(0,
k,
j, nr-1) = 0 ;
552 nr =
mgrid.get_nr(1) ;