Parallel Krylov adjoint sensitivity analysis example
2D Advection diffusion PDE with homogeneous Neumann B.C.
Computes gradient of G = int_t_Omega ( c_i^2 ) dt dOmega
with respect to the source values at each grid point.
 
Domain:
   0.000000 < x < 20.000000   mx = 20  npe_x = 2
   0.000000 < y < 20.000000   my = 40  npe_y = 2
 
Begin forward integration... done.   G = 3.723818e+03
 
Final Statistics..
 
lenrw   =   8746     leniw =    212
llrw    =   8656     lliw  =     80
nst     =    104
nfe     =    108     nfel  =    126
nni     =    105     nli   =    126
nsetups =     16     netf  =      0
npe     =      2     nps   =    215
ncfn    =      0     ncfl  =      0
 
Begin backward integration... done.
 
Final Statistics..
 
lenrw   =  17316     leniw =    212
llrw    =   8656     lliw  =     80
nst     =     78
nfe     =     90     nfel  =    138
nni     =     87     nli   =    138
nsetups =     17     netf  =      0
npe     =      2     nps   =    217
ncfn    =      0     ncfl  =      0