> for (PetscInt i = 0; i < nConstants; ++i) { > > VecRestoreArray(constants[i], &rawConstants[i]); > > VecAssemblyBegin(constants[i]); > > VecAssemblyEnd(constants[i]); > > VecNormalize(constants[i], nullptr); >
This does not look correct. > } >
> for (PetscInt i = 0; i < nConstants; ++i) { > > VecRestoreArray(constants[i], &rawConstants[i]); > > VecAssemblyBegin(constants[i]); > > VecAssemblyEnd(constants[i]); > > VecNormalize(constants[i], nullptr); >
This does not look correct. > } >