80 if (
aaa.get_etat() == ETATZERO) {
86 int nz =
mgrid.get_nzone();
102 int nt =
mgrid.get_nt(0) ;
103 int np =
mgrid.get_np(0) ;
110 int dzpuis =
source.get_dzpuis();
114 source.set_spectral_va().ylm() ;
121 eta_tilde.set_spectral_va().set_etat_cf_qcq() ;
122 eta_tilde.set_spectral_va().c_cf->annule_hard() ;
124 vr.set_spectral_base(base) ;
125 vr.set_spectral_va().set_etat_cf_qcq() ;
126 vr.set_spectral_va().c_cf->annule_hard() ;
137 for (
int l=0 ;
l<nz ;
l++) {
148 np =
mgrid.get_np(0) ;
149 nt =
mgrid.get_nt(0) ;
156 for (
int k=0 ;
k<np+1 ;
k++)
157 for (
int j=0 ;
j<nt ;
j++)
158 if (nullite_plm(
j, nt,
k, np, base) == 1) {
170 nr =
mgrid.get_nr(0) ;
171 alpha = (*mp_aff).get_alpha()[0] ;
180 if (
source.get_mp().get_mg()->get_type_r(0) == RARE) {
220 else if (l_quant == 1) {
254 if (
source.get_mp().get_mg()->get_type_r(0) == RARE) {
295 else if (l_quant == 1) {
330 if (
source.get_mp().get_mg()->get_type_r(0) == RARE) {
358 if (
source.get_mp().get_mg()->get_type_r(0) == RARE) {
390 for (
int l=1 ;
l<nz-1 ;
l++) {
393 alpha = (*mp_aff).get_alpha()[
l] ;
394 beta = (*mp_aff).get_beta()[
l] ;
444 for (
int i=0 ;
i<nr ;
i++)
488 nr =
mgrid.get_nr(nz-1) ;
489 alpha = (*mp_aff).get_alpha()[nz-1] ;
490 beta = (*mp_aff).get_beta()[nz-1] ;
583 cout <<
"Unknown dzpuis in vector_divfree_A ..." <<
endl ;
593 indic = alpha*alpha ;
621 for (
int l=0 ;
l<nz ;
l++) {
623 for (
int i=0 ;
i<nr ;
i++) {
635 if (
vr.set_spectral_va().c != 0x0)
636 delete vr.set_spectral_va().c ;
637 vr.set_spectral_va().c = 0x0 ;
638 vr.set_spectral_va().ylm_i() ;
void give_quant_numbers(int, int, int, int &, int &, int &) const
Computes the various quantum numbers and 1d radial base.
void sol_Dirac_A_tau(const Scalar &aaa, Scalar &eta, Scalar &vr, const Param *par_bc=0x0) const
Solves via a tau method a system of two-coupled first-order PDEs obtained from the divergence-free co...