23char division_xpun_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Non_class_members/Operators/division_xpun.C,v 1.4 2014/10/13 08:53:23 j_novak Exp $" ;
62Cmp division_xpun (
const Cmp& source,
int num_front) {
64 assert (source.get_etat() != ETATNONDEF) ;
65 Cmp resultat (source) ;
67 if (resultat.get_etat() == ETATZERO)
71 resultat.va.set_etat_cf_qcq() ;
72 int base_r = source.va.base.b[num_front+1] &
MSQ_R ;
74 int nr = source.get_mp()->get_mg()->get_nr(num_front+1) ;
75 double* coef =
new double[nr] ;
77 for (
int k=0 ; k<source.get_mp()->get_mg()->get_np(num_front+1)+1 ; k++)
79 for (
int j=0 ; j<source.get_mp()->get_mg()->get_nt(num_front) ; j++) {
80 for (
int i=0 ; i<nr ; i++)
81 coef[i] = (*resultat.va.c_cf)(num_front+1, k, j, i) ;
82 sxpun_1d (nr, &coef, base_r) ;
83 for (
int i=0 ; i<nr ; i++)
84 resultat.va.c_cf->set(num_front+1, k, j, i) = coef[i] ;
#define MSQ_R
Extraction de l'info sur R.