Dear All,
Could someone help to answer a question about snes/ex22.c?
I runed the program twice with two different option sets.
In the first run, I used:
mpirun -np 1 ./ex22 -da_grid_x 10 -snes_monitor
this will invoke all the matrix_free_options at line 96. It is shown that
on all grid levels
Ryan,
This is because the explicit Jacobian computed in ex22.c is
not correct. Run with -snes_type test -dmmg_nlevels 1
This is because the code cannot determine the nonzero structure
of the Jacobian and (also) hence a coloring of the Jacobian.
If you run with