[dune-pdelab] [dune-pdelab-commit] dune-pdelab r1988 - in trunk: . dune/pdelab/common dune/pdelab/finiteelementmap dune/pdelab/localoperator dune/pdelab/newton dune/pdelab/stationary

[Christian Engwer]

Hi Markus,

> Modified: trunk/dune/pdelab/newton/newton.hh
> ==============================================================================
> --- trunk/dune/pdelab/newton/newton.hh	Sun Apr 22 13:40:26 2012	(r1987)
> +++ trunk/dune/pdelab/newton/newton.hh	Sun Apr 22 15:43:51 2012	(r1988)
> @@ -13,6 +13,7 @@
>  #include <dune/common/ios_state.hh>
>  #include <dune/common/timer.hh>
>  
> +#include<dune/istl/matrixmarket.hh>
>  namespace Dune
>  {
>      namespace PDELab
> @@ -164,6 +165,8 @@
>                  if (this->verbosity_level >= 4)
>                      std::cout << "      Solving linear system..." << std::endl;
>                  z = 0.0;                                        // TODO: vector interface
> +                storeMatrixMarket(A.base(), std::string("A"));
> +                storeMatrixMarket(r.base(), std::string("z"));
>                  this->solver.apply(A, z, r, this->linear_reduction);        // TODO: solver interface
>  
>                  ios_base_all_saver restorer(std::cout); // store old ios flags
> 

are you sure this is what you wanted to commit?

I guess it is debug code, or you should add a switch to dis-/en-able
the matrix market output.

Christian