#include <boost/numeric/ublas/matrix.hpp>
#include <boost/numeric/ublas/io.hpp>
// 
// Written by:
// -- 
// John L. Weatherwax                2010-02-25
// 
// email: wax@alum.mit.edu
// 
// Please send comments and especially bug reports to the
// above email address.
// 
//-----

using namespace boost::numeric::ublas;

void scale_regularization( matrix<double> & M, double gamma ){
  assert( (0.0 <= gamma) && (gamma <= 1.0) ); 
  
  size_t N = M.size1();
  for( size_t ii=0; ii<N; ++ii ){
    for( size_t jj=0; jj<N; ++jj ){
      if( ii==jj ) continue;  
      M(ii,jj) *= (1.0-gamma);         // modify the non-diagonal elements 
    }
  }
}