Dear all

I learnt from https://github.com/mbarzegary/navier-stokes-solver an efficient @prj @mojtaba.barzegari an efficient PETSC solver for Stokes solution

[u, v, w, p] in FEM space [P2, P2, P2, P1]:

set(A, sparams = "-ksp_monitor -ksp_type fgmres -ksp_converged_reason "

+ "-pc_type fieldsplit -pc_fieldsplit_type schur "

+ "-fieldsplit_velocity_pc_type gamg -fieldsplit_velocity_pc_gamg_sym_graph true -fieldsplit_pressure_ksp_max_it 5 "

+ “-fieldsplit_pressure_pc_type jacobi -fieldsplit_velocity_ksp_type preonly -pc_fieldsplit_schur_fact_type full”,

fields = ux, names = names);

with

string[int] names(2);

names[0] = “velocity”;

names[1] = “pressure”;

Now, I have **two velocity fields** and **two pressure fields**: [u1, v1, w1, u2, v2, w2, p1, p2] in a FEM space [P2,P2,P2,P2,P2,P2, P1,P1], **is the parameter set-up still the same?**

BTW, the second velocity or pressure field is from the adjoint space.

Many thanks in advance!