Location Matrix with vectorial space

General Discussion
1

Hello,

I’m trying to build a location matrix for a FE vectorial space. Am I doing it right like this?

fespace LambdaH(calP, [Pm, Pm]);
LambdaH [l10, l20];
LambdaH [l1X, l2X] = [x,x], [l1Y, l2Y] = [y,y];

int nDoFl0 = LambdaH.ndof;
int nDoFl0K = LambdaH.ndofK;

// cout << nDoFl0 << "\n";
// cout << nDoFl0K << "\n";

real[int,int] coordbL(nDoFl0,2);
for(int i=0; i<nDoFl0; i++){
    coordbL(i,0) = l1X[][i];
    coordbL(i,1) = l1Y[][i];
}
// cout << coordbL << "\n";

//-------------------------------------LM matrix
int[int,int] locationMatrix(NbTriangles,nDoFl0K);
for(int k=0; k < NbTriangles; k++){
    for(int i=0;i<nDoFl0K;i++){
        locationMatrix(k,i) = LambdaH(k,i);
    }
}
// cout << "LM array Lambda0" << "\n"; 
// cout << locationMatrix << "\n";