28char sol_Dirac_A_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Tensor/vector_divfree_A.C,v 1.6 2014/10/13 08:53:44 j_novak Exp $" ;
85 if (
aaa.get_etat() == ETATZERO) {
91 int nz =
mgrid.get_nzone() ;
106 int nt =
mgrid.get_nt(0) ;
107 int np =
mgrid.get_np(0) ;
114 source.set_spectral_va().ylm() ;
121 tilde_vr.set_spectral_va().set_etat_cf_qcq() ;
122 tilde_vr.set_spectral_va().c_cf->annule_hard() ;
125 tilde_eta.set_spectral_va().set_etat_cf_qcq() ;
126 tilde_eta.set_spectral_va().c_cf->annule_hard() ;
146 for (
int k=0 ;
k<np+1 ;
k++) {
147 for (
int j=0 ;
j<nt ;
j++) {
150 if ( (nullite_plm(
j, nt,
k, np, base) == 1) && (
l_q > 0)) {
178 ope.set(2*nr-1, nr) = 1;
185 ope.set(2*nr-1,
col) = 0 ;
196 ope.set(nr-1, nr-1) = 1 ;
197 ope.set(2*nr-1, 2*nr-1) = 1 ;
212 sec.set(2*nr-1) = 0 ;
218 for (
int i=0;
i<nr;
i++) {
227 for (
int i=0;
i<nr;
i++) {
249 for (
int k=0 ;
k<np+1 ;
k++) {
250 for (
int j=0 ;
j<nt ;
j++) {
253 if ( (nullite_plm(
j, nt,
k, np, base) == 1) && (
l_q > 0)) {
293 for (
int i=0;
i<nr;
i++) {
304 for (
int i=0;
i<nr;
i++) {
314 for (
int i=0;
i<nr;
i++) {
334 for (
int k=0 ;
k<np+1 ;
k++) {
335 for (
int j=0 ;
j<nt ;
j++) {
338 if ( (nullite_plm(
j, nt,
k, np, base) == 1) && (
l_q > 0)) {
384 for (
int i=0;
i<nr;
i++) {
391 for (
int i=0;
i<nr;
i++) {
400 int taille = 2*
nz_bc + 1;
401 if (
cedbc) taille-- ;
435 for (
int k=0 ;
k<np+1 ;
k++)
436 for (
int j=0 ;
j<nt ;
j++) {
438 if ((nullite_plm(
j, nt,
k, np, base) == 1) && (
l_q > 0)) {
445 if ((nz==1)&&(
mgrid.get_type_r(0) == RARE)) {
564 nr =
mgrid.get_nr(0) ;
565 for (
int i=0 ;
i<nr ;
i++) {
574 for (
int i=0 ;
i<nr ;
i++) {
587 for (
int i=0 ;
i<nr ;
i++) {
600 if (
tilde_vr.set_spectral_va().c != 0x0)
601 delete tilde_vr.set_spectral_va().c ;
602 tilde_vr.set_spectral_va().c = 0x0 ;
603 tilde_vr.set_spectral_va().ylm_i() ;
605 if (
tilde_eta.set_spectral_va().c != 0x0)
Bases of the spectral expansions.
void mult_x()
The basis is transformed as with a multiplication by .
void give_quant_numbers(int, int, int, int &, int &, int &) const
Computes the various quantum numbers and 1d radial base.
Class for the elementary differential operator (see the base class Diff ).
Class for the elementary differential operator Identity (see the base class Diff ).
Class for the elementary differential operator division by (see the base class Diff ).
Class for the elementary differential operator (see the base class Diff ).
Time evolution with partial storage (*** under development ***).
const double * get_beta() const
Returns the pointer on the array beta.
const double * get_alpha() const
Returns the pointer on the array alpha.
const Mg3d * get_mg() const
Gives the Mg3d on which the mapping is defined.
Coefficients storage for the multi-domain spectral method.
Tensor field of valence 0 (or component of a tensorial field).
void sol_Dirac_A(const Scalar &aaa, Scalar &eta, Scalar &vr, const Param *par_bc=0x0) const
Solves a system of two-coupled first-order PDEs obtained from the divergence-free condition and the r...
#define R_CHEBP
base de Cheb. paire (rare) seulement
const Map *const mp
Mapping on which the numerical values at the grid points are defined.