#ifndef  IntfLapack_H_SEEN
#define  IntfLapack_H_SEEN

#include "machdefs.h"
#include "tarray.h"

namespace SOPHYA {

template <class T>
class LapackServer { 
public:
  TArray<T>& LinSolve(TArray<T>& a, TArray<T> & b); 
};

template <class T>
inline TArray<T>& LapackLinSolve(TArray<T>& a, TArray<T> & b)
  { LapackServer<T> lps; return( lps.LinSolve(a, b) );  }


} // Fin du namespace

void rztest_lapack(TArray<r_4>& a, TArray<r_4>& b);

#endif