C++ et éléments finis Note de cours de DEA (version provisoire)

38 CHAPITRE 5. MODÉLISATION DES VECTEURS, MATRICES, TENSEURS
} ;
void MatriceLaplacien1D::addMatMul(const KN_ & x, KN_ & Ax) const {
int n= x.N(),n_1=n-1 ;
double h=1./(n_1), h2= h*h, d = 2/h, d1 = -1/h ;
R Ax0=Ax[0], Axn_1=Ax[n_1] ;
Ax=0 ;
for (int i=1 ;i< n_1 ; i++)
Ax[i] = (x[i-1] +x[i+1]) * d1 + x[i]*d ;
Ax[0]=x[0] ;
Ax[n_1]=x[n_1] ;
}
int main(int argc,char ** argv)
{
typedef KN Rn ;
typedef KN_ Rn_ ;
typedef KNM Rnm ;
typedef KNM_ Rnm_ ;
}
// CL
{ int n=10 ;
Rnm A(n,n),C(n,n),Id(n,n);
A=-1 ;
C=0 ;
Id=0 ;
Rn_ Aii(A,SubArray(n,0,n+1)); // la diagonal de la matrice A sans copy
Rn_ Cii(C,SubArray(n,0,n+1)); // la diagonal de la matrice C sans copy
Rn_ Idii(Id,SubArray(n,0,n+1)); // la diagonal de la matrice Id sans copy
for (int i=0 ;i