23char valeur_val_propre_1d_C[] =
"$Header: /cvsroot/Lorene/C++/Source/Valeur/valeur_val_propre_1d.C,v 1.3 2014/10/13 08:53:51 j_novak Exp $" ;
55#include "type_parite.h"
69 so.set_etat_cf_qcq() ;
75 int nt =
so.mg->get_nt(0) ;
89 for (
int i=0 ;
i<nt ;
i++)
90 for (
int j=0 ;
j<nt ;
j++)
94 for (
int j=0 ;
j<nt ;
j++) {
96 for (
int i=1 ;
i<
j ;
i++)
107 inv->set_etat_qcq() ;
110 auxi.set_etat_qcq() ;
112 for (
int i=0 ;
i<nt ;
i++) {
113 for (
int j=0 ;
j<nt ;
j++)
117 for (
int j=0 ;
j<nt ;
j++)
123 for (
int l=0 ;
l<
so.mg->get_nzone() ;
l++)
124 if (
so.c_cf->t[
l]->get_etat() != ETATZERO)
125 for (
int k=0 ;
k<
so.mg->get_np(
l) ;
k++)
126 for (
int i=0 ;
i<
so.mg->get_nr(
l) ;
i++) {
129 coefs.set_etat_qcq() ;
130 for (
int j=0 ;
j<nt ;
j++)
133 prod.set_etat_qcq() ;
134 for (
int j=0 ;
j<nt ;
j++)
139 for (
int j=0 ;
j<nt ;
j++)
140 for (
int jb=0 ;
jb<nt ;
jb++)
145 for (
int j=0 ;
j<nt ;
j++)
146 for (
int jb=0 ;
jb<nt ;
jb++)
150 for (
int j=0 ;
j<nt ;
j++)
168 cout <<
"Problem in rotate_propre_pair" <<
endl ;
182 cout <<
"Problem in rotate_propre_pair" <<
endl ;
198 so.set_etat_cf_qcq() ;
204 int nt =
so.mg->get_nt(0) ;
218 for (
int i=0 ;
i<nt ;
i++)
219 for (
int j=0 ;
j<nt ;
j++)
223 for (
int j=0 ;
j<nt ;
j++) {
224 for (
int i=0 ;
i<
j ;
i++)
235 inv->set_etat_qcq() ;
238 auxi.set_etat_qcq() ;
240 for (
int i=0 ;
i<nt ;
i++) {
241 for (
int j=0 ;
j<nt ;
j++)
245 for (
int j=0 ;
j<nt ;
j++)
252 for (
int l=0 ;
l<
so.mg->get_nzone() ;
l++)
253 if (
so.c_cf->t[
l]->get_etat() != ETATZERO)
254 for (
int k=0 ;
k<
so.mg->get_np(
l) ;
k++)
255 for (
int i=0 ;
i<
so.mg->get_nr(
l) ;
i++) {
258 coefs.set_etat_qcq() ;
259 for (
int j=0 ;
j<nt ;
j++)
262 prod.set_etat_qcq() ;
263 for (
int j=0 ;
j<nt ;
j++)
268 for (
int j=0 ;
j<nt ;
j++)
269 for (
int jb=0 ;
jb<nt ;
jb++)
274 for (
int j=0 ;
j<nt ;
j++)
275 for (
int jb=0 ;
jb<nt ;
jb++)
279 for (
int j=0 ;
j<nt ;
j++)
297 cout <<
"Problem in rotate_propre_impair" <<
endl ;
311 cout <<
"Problem in rotate_propre_impair" <<
endl ;
325 int nz =
mg->get_nzone() ;
328 for (
int i=1 ;
i<nz ;
i++)
353 cout <<
"Unknown basis in Valeur::val_propre_1d" <<
endl ;
361 int nz =
mg->get_nzone() ;
364 for (
int i=1 ;
i<nz ;
i++)
389 cout <<
"Unknown basis in Valeur::val_propre_1d_i" <<
endl ;
Bases of the spectral expansions.
int get_base_t(int l) const
Returns the expansion basis for functions in the domain of index l (e.g.
Time evolution with partial storage (*** under development ***).
Values and coefficients of a (real-value) function.
void val_propre_1d_i()
Inverse transformation of val_propre_1d.
friend void rotate_propre_impair(Valeur &, bool)
Friend fonction.
const Mg3d * mg
Multi-grid Mgd3 on which this is defined.
friend void rotate_propre_pair(Valeur &, bool)
Friend fonction.
Base_val base
Bases on which the spectral expansion is performed.
void val_propre_1d()
Set the basis to the eigenvalues of .
#define T_CL_COS_I
CL of odd cosines.
#define T_CL_SIN_I
CL of odd sines.
#define T_SIN_P
dev. sin seulement, harmoniques paires
#define T_CL_SIN_P
CL of even sines.
#define T_COS_P
dev. cos seulement, harmoniques paires
#define T_CL_COS_P
CL of even cosines.
#define T_SIN_I
dev. sin seulement, harmoniques impaires
#define T_COS_I
dev. cos seulement, harmoniques impaires