Hi, So in ex62 of SNES it solves the stokes equation with P2-P1 elements. I want to solve the darcy equation with P2-P1 elements in a similar fashion.
The difference between stokes and darcy is that in the weak form of darcy, the velocity-velocity term is v*u whereas in stokes it's grad[v]*grad[u] In ex62, the function g3_uu (lines 172-183) formulates the grad[v]*grad[u] term but how would I do it for Darcy's equation where its simply v*u (assuming drag coefficient is scalar quantity of 1.0)? Thanks, Justin
