>         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.


>         }
>

Reply via email to