Changeset 26131 for issm/trunk-jpl/src/c/analyses/AdjointBalancethicknessAnalysis.cpp

Timestamp:

03/23/21 11:52:52 (4 years ago)

Author:

Mathieu Morlighem

Message:

CHG: big simplification here: we now carry the vectors of fset and sset using the same size as gset, with -1 when they are not in the right set. Example: old_fset = [1 2 3], new_fset=[1 -1 2 3 -1 -1]; enjoy

File:

• issm/trunk-jpl/src/c/analyses/AdjointBalancethicknessAnalysis.cpp (modified) (1 diff)

issm/trunk-jpl/src/c/analyses/AdjointBalancethicknessAnalysis.cpp

@@ -227 +227 @@
         /*Transpose and return Ke*/
         Ke->Transpose();
-        _assert_(Ke->nrows == numvertices && Ke->ncols == numvertices);
+        _assert_(Ke->nrows == numvertices);
 
         for(int i=0;i<numvertices;i++){
                 for(int j=0;j<numvertices;j++){
-                        ge[i] += Ke->values[i*Ke->ncols + j] * lambda[j];
+                        ge[i] += Ke->values[i*Ke->nrows + j] * lambda[j];
                 }
                 //ge[i]=-lambda[i];