// norm.cc -- compute norms
// $Id$

//
// *** class jtutil::norm::
//

#include <cmath>
#include <algorithm>
#include <assert.h>

#include "util.hh"

//*****************************************************************************

// specify data point
namespace jtutil
	  {
template <typename fp>
  void norm<fp>::data(fp x)
{
++N_;
sum2_ += x*x;
infinity_norm_ = std::max(infinity_norm_, jtutil::abs<fp>(x));
}
	  }	// namespace jtutil::

//******************************************************************************

// get norms
namespace jtutil
	  {
template<typename fp>
  fp norm<fp>::two_norm() const { return std::sqrt(sum2_); }
template<typename fp>
  fp norm<fp>::rms_norm() const
	{ assert(is_nonempty()); return std::sqrt(sum2_/fp(N_)); }
template<typename fp>
  fp norm<fp>::infinity_norm() const { return infinity_norm_; }
	  }	// namespace jtutil::

//******************************************************************************

//
// ***** template instantiations *****
//

template class jtutil::norm<float>;
template class jtutil::norm<double>;