Parallelize Navier Stokes with moving boundary

Dear all,

I already implement the sequential version of NS equation with moving boundary in FreeFem++, the code is
sequential_moving_boundary.edp (2.4 KB)

Now I want to parallelize it using PETSc, just change little, for example changing matrix to Mat. And it shows error like this on my Windows machine

and the parallel version code is
parallel_moving_boundary.edp (2.5 KB)

Can anyone give me any advice? Many thanks.

Your A is uninitialized line 59, so you need to fix that.

Thank you for your reminder, Professor. It’s a stupid mistake. I have fixed the initialization of Mat like this

However, there still exists error as shown below

Do you know what’s wrong? Thanks a lot.
parallel_moving_boundary_2.edp (2.6 KB)

You are trying to initialize with the same macro a Mat with two fields and another one with three fields, hence the error.