
2D elliptic PDE on unit square
   d^2 u / dx^2 + d^2 u / dy^2 = u^3 - u + 2.0
 + homogeneous Dirichlet boundary conditions

Solution method: Anderson accelerated Picard iteration with band linear solver.
Problem size: 31 x 31 =  961

Computed solution (||F|| = 1.83409e-13):

            0.03125  0.12500  0.21875  0.31250  0.40625  0.50000  0.59375  0.68750  0.78125  0.87500  0.96875  

0.03125     0.00405  0.01165  0.01617  0.01896  0.02051  0.02100  0.02051  0.01896  0.01617  0.01165  0.00405  
0.12500     0.01165  0.03772  0.05461  0.06530  0.07126  0.07318  0.07126  0.06530  0.05461  0.03772  0.01165  
0.21875     0.01617  0.05461  0.08098  0.09813  0.10780  0.11093  0.10780  0.09813  0.08098  0.05461  0.01617  
0.31250     0.01896  0.06530  0.09813  0.11989  0.13229  0.13631  0.13229  0.11989  0.09813  0.06530  0.01896  
0.40625     0.02051  0.07126  0.10780  0.13229  0.14632  0.15089  0.14632  0.13229  0.10780  0.07126  0.02051  
0.50000     0.02100  0.07318  0.11093  0.13631  0.15089  0.15564  0.15089  0.13631  0.11093  0.07318  0.02100  
0.59375     0.02051  0.07126  0.10780  0.13229  0.14632  0.15089  0.14632  0.13229  0.10780  0.07126  0.02051  
0.68750     0.01896  0.06530  0.09813  0.11989  0.13229  0.13631  0.13229  0.11989  0.09813  0.06530  0.01896  
0.78125     0.01617  0.05461  0.08098  0.09813  0.10780  0.11093  0.10780  0.09813  0.08098  0.05461  0.01617  
0.87500     0.01165  0.03772  0.05461  0.06530  0.07126  0.07318  0.07126  0.06530  0.05461  0.03772  0.01165  
0.96875     0.00405  0.01165  0.01617  0.01896  0.02051  0.02100  0.02051  0.01896  0.01617  0.01165  0.00405  

Final Statistics.. 

nni      =      6    nfe     =      7 
nje      =      6    nfeB    =      0 

lenrwB   = 150877    leniwB  =    961 
