Commit 96f2739cd71eb7c1ea5b51837ad1e3ac09761921

fix bug in computeBJacobian
isopatch.h
(8 / 7)
  
150150 for (int k=0; k<ldim; k++)
151151 bJ(k,i) = 0.0;
152152
153 int lk;
153154 for (int i=0; i<dim; i++)
154155 for (int j=0; j<getNumLocalBasis(); j++) {
155 for (int k=0; k<dim; k++)
156 if (k<ndir) {
157 bJ(k,i) += geom->coord(getLocalBasisGlobalIndex(j),i)*rationalgrad(j)[k];
158 } else if (k>ndir) {
159 bJ(k-1,i) += geom->coord(getLocalBasisGlobalIndex(j),i)*rationalgrad(j)[k-1];
160 }
156 for (int k=0; k<dim; k++) {
157 if (k==ndir) continue;
158 if (k<ndir) { lk = k; } else { lk = k-1; }
159 bJ(lk,i) += geom->coord(getLocalBasisGlobalIndex(j),i)*rationalgrad(j)[k];
160 }
161161 }
162
162163}
163164
164165
341341 i /= (ne[j]-2*deg[j]);
342342 }
343343 }
344
344
345345 computeLocalIndices(); // update indices and gindices
346346
347347 computeRationalBasis();