#ifndef LOG_DET_H_
#define LOG_DET_H_

#include <functions/ScalarFunc.h>

/**
 * @short Log determinant 
 * LogDet calculates the log determinant of a square matrix.  The
 * function assumes that the matrix is symmetric positive definite.
 * but currently does not test this.
 * <pre>
 * y <- logdet(x)
 * y <- log|x| for x an n x n symmetric positive definite matrix
 * </pre>
 */
class LogDet : public ScalarFunc
{
public:
  LogDet ();

  double eval (std::vector <SArray const *> const &args) const;
  bool checkParameterDim (std::vector <SArray const *> const &args) const;
};

#endif /* LOG_DET_H_ */