mathUtils.h

//
// Created by milinda on 3/31/17.
//

#ifndef SFCSORTBENCH_MATHUTILS_H
#define SFCSORTBENCH_MATHUTILS_H

#include <iostream>
#include <cmath>
#include <math.h>
#include "parUtils.h"

template <typename T>
T normL2(T * vec1,T* vec2, unsigned int n,MPI_Comm comm);

template <typename T>
T normLInfty(T *vec1, T *vec2, unsigned int n, MPI_Comm comm);


template <typename T>
T normLInfty(T *vec1, T *vec2, unsigned int n);



template <typename T>
T normL2(T * vec,unsigned int n, MPI_Comm comm);


template <typename T>
T normL2(T * vec1,T* vec2, unsigned int n);

template <typename T>
T normL2(T * vec,unsigned int n);


template <typename T>
T normLInfty(T * vec,unsigned int n);


template <typename T>
T normLInfty(T * vec,unsigned int n,MPI_Comm comm);

template <typename T>
T vecMin(T * vec,unsigned int n);


template <typename T>
T vecMax(T * vec,unsigned int n);


template <typename T>
T vecMin(T * vec,unsigned int n,MPI_Comm comm);


template <typename T>
T vecMax(T * vec,unsigned int n,MPI_Comm comm);

template <typename T>
T dot(const T* v1, const T*v2,const unsigned int n);

template <typename T>
T dot(const T* v1, const T*v2,const unsigned int n,MPI_Comm comm);


template<typename T>
void mul(const T alpha, const T* v, const unsigned int n, T* out);


template <typename T>
T add(const T* v1, const T*v2,const unsigned int n, T* out);

template <typename T>
T subt(const T* v1, const T*v2,const unsigned int n, T* out);

template<typename T>
void kron(const T* M1, const T* M2, T* out, unsigned int r1,unsigned int c1, unsigned int r2, unsigned int c2);




#include "mathUtils.tcc"


#endif //SFCSORTBENCH_MATHUTILS_H