module solver_report_mod{
  use ddot_mod
  real, private  :: scale
contains
  subroutine solver_report(iter,x,g,rd,rm) {
    integer,            intent(in) :: iter
    real, dimension(:), intent(in) :: x,g
    real, dimension(:), pointer    :: rd,rm
    optional                       ::    rm
    double precision               :: dd,dm,dx,dg

    dd=dot_product(rd,rd)
    dx=dot_product( x, x)
    dg=dot_product( g, g)
    dm=0.; if(present(rm)) dm=dot_product(rm,rm)
    if(iter==1){
       scale=max(dd,dm)/1000.
       write(0,*) "DOT PRODUCT VALUES SCALED BY ",scale
    }
    write (0,*) "iteration ", iter &
    ,           "   res dat", int(dd/scale) &
    ,           "   res mod", int(dm/scale) &
    ,           "       mod", int(dx/scale) &
    ,           "      grad", int(dg/scale)

  }
}