"transfer" in PETSc-complex interface

[yuki]

Dear all,

I would like to use “transfer” function in PETSc-complex interface.
In FreeFEM v 4.8, the following code worked well:

//load “PETSc”
load “PETSc-complex”
macro dimension()2// EOM
include “macro_ddm.idp”

meshN Sh1=square(30,30,[x,y]);
meshN Sh2=square(50,50,[x,y]);

buildDmesh(Sh1)
buildDmesh(Sh2)

fespace Ph1(Sh1, P2);
Ph1 u=sin(4x)+cos(5y);

fespace Ph2(Sh2, P2);
Ph2 v;
// v=u;

macro defPlot(u)u//
transfer(Sh1, P2, u, Sh2, P2, v);

This code was borrowed from the previous topic:

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…

However, I could not run the above example with FreeFEM v 4.9.
Would you have any solutions to this issue?

Thank you so much in advance,
Best regards,

Yuki

[prj]

You need this commit Fixes renumv + adds PetscScalar in .edp · FreeFem/FreeFem-sources@723520e · GitHub which is in 4.9, so please consider updating your FreeFEM installation.
Also, your FE function needs to live in Ph1<complex> (resp. Ph2<complex>), not Ph1 (resp. Ph2).

[yuki]

Thank you for your advice.
I could transfer the solution in FreeFEM v 4.9 by changing FE spaces from Phk to Phk (k=1,2).