Call PETSc matrix operator MatDiagonalSet (…)

[prj]

Copy variables between different meshes with PETSc

Dear Sir or Madam, I would like to use different 2 meshes with PETSc. Additionally, variables want to be copied from each other. However, as shown following codes, the direct copy, such as “v=u”, does not work. Would you please let me know how to do it? Thank you very much in advance. Best regards, Tak // $ ff-mpirun -np 4 script.edp -wg -ne load “PETSc” macro dimension()2// EOM include “macro_ddm.idp” meshN Sh1=square(30,30,[x,y]); meshN Sh2=square(50,50,[x,y]); buildDmesh(Sh1) b…