vnl_trace.cxx

InsightToolkit-1.4.0(有大量的优化算法程序)

// This is vxl/vnl/vnl_trace.cxx
#ifdef VCL_NEEDS_PRAGMA_INTERFACE
#pragma implementation
#endif
//:
// \file
// \author fsm

#include "vnl_trace.h"

#include <vnl/vnl_matrix.h>

template <class T>
T vnl_trace(vnl_matrix<T> const &M) {
  T sum(0);
  unsigned N = (M.rows() < M.cols() ? M.rows() : M.cols());
  for (unsigned i=0; i<N; ++i)
    sum += M(i, i);
  return sum;
}

#define VNL_TRACE_INSTANTIATE(T) \
template T vnl_trace(vnl_matrix<T > const &)

VNL_TRACE_INSTANTIATE(float);
VNL_TRACE_INSTANTIATE(double);