Fix bug: According to the mannual of PETSc, the values are given in row major order by default, which is the opposite of the Fortran convention.
The stiffness matrix of Cahn–Hilliard type equation is unsymmetrical, so we need to transpose the element stiffness matrix before assembling.
This commit is contained in:
parent
e20271b89a
commit
775146030c