//-*- mode: c++; indent-tabs-mode: t; coding: utf-8; show-trailing-whitespace: t -*-

// file algebra.hpp

#ifndef _ALGEBRA_HPP
#define _ALGEBRA_HPP

// System includes.

#include <iostream>
#include <iomanip>
#include <cmath>

// Local includes.

#include "blkcom.hpp"
#include "utilities.hpp"
#include "movecopy.hpp"

namespace algebra {
  
  bool cdivz(std::vector<double> &, \
	     std::vector<double> &, \
	     const double &,	    \
	     const double &,	    \
	     const double &,	    \
	     const double &,	    \
	     const long int &);
  bool cmultz(std::vector<double> &,		\
	      std::vector<double> &,		\
	      const double &,			\
	      const double &,			\
	      const double &,			\
	      const double &,			\
	      const long int &);
  bool trgwnd(const double &, double *);
  bool addmxd(std::vector<double> &,		\
	      const double &,			\
	      std::vector<double> &,		\
	      const size_t &);
  bool multmx(std::vector<double> &, \
	      std::vector<double> &, \
	      std::vector<double> &, \
	      std::vector<double> &, \
	      const size_t &);
  bool mult(std::vector<double> &, \
	    std::vector<double> &, \
	    std::vector<double> &, \
	    const size_t &,	   \
	    long int &);
  void dgelg(std::vector<double> &, \
	     std::vector<double> &, \
	     const size_t &,	    \
	     const size_t &,	    \
	     const float &,	    \
	     long int &);
  void matmul(float [ 3 ][ 3 ], float [ 3 ][ 3 ]);
  void matvec(float [ 3 ][ 3 ], float [ 15 ]);
  
}

#endif // _ALGEBRA_HPP

// end of file algebra.hpp