FXMat4d.h

/********************************************************************************
*                                                                               *
*            D o u b l e - P r e c i s i o n   4 x 4   M a t r i x              *
*                                                                               *
*********************************************************************************
* Copyright (C) 1994,2024 by Jeroen van der Zijp.   All Rights Reserved.        *
*********************************************************************************
* This library is free software; you can redistribute it and/or modify          *
* it under the terms of the GNU Lesser General Public License as published by   *
* the Free Software Foundation; either version 3 of the License, or             *
* (at your option) any later version.                                           *
*                                                                               *
* This library is distributed in the hope that it will be useful,               *
* but WITHOUT ANY WARRANTY; without even the implied warranty of                *
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the                 *
* GNU Lesser General Public License for more details.                           *
*                                                                               *
* You should have received a copy of the GNU Lesser General Public License      *
* along with this program.  If not, see <http://www.gnu.org/licenses/>          *
********************************************************************************/
#ifndef FXMAT4D_H
#define FXMAT4D_H

namespace FX {


class FXQuatd;
class FXMat3d;


class FXAPI FXMat4d {
protected:
  FXVec4d m[4];
public:

  FXMat4d(){}

  FXMat4d(FXdouble s);

  FXMat4d(const FXMat3d& s);

  FXMat4d(const FXMat4d& s);

  FXMat4d(const FXdouble s[]);

  FXMat4d(FXdouble a,FXdouble b,FXdouble c,FXdouble d);

  FXMat4d(FXdouble a00,FXdouble a01,FXdouble a02,FXdouble a03,
          FXdouble a10,FXdouble a11,FXdouble a12,FXdouble a13,
          FXdouble a20,FXdouble a21,FXdouble a22,FXdouble a23,
          FXdouble a30,FXdouble a31,FXdouble a32,FXdouble a33);

  FXMat4d(const FXVec4d& a,const FXVec4d& b,const FXVec4d& c,const FXVec4d& d);

  FXMat4d(const FXQuatd& quat);

  FXMat4d& operator=(FXdouble s);

  FXMat4d& operator=(const FXMat3d& s);
  FXMat4d& operator=(const FXMat4d& s);

  FXMat4d& operator=(const FXQuatd& quat);

  FXMat4d& operator=(const FXdouble s[]);

  FXMat4d& set(FXdouble s);

  FXMat4d& set(const FXMat3d& s);

  FXMat4d& set(const FXMat4d& s);

  FXMat4d& set(const FXdouble s[]);

  FXMat4d& set(FXdouble a,FXdouble b,FXdouble c,FXdouble d);

  FXMat4d& set(FXdouble a00,FXdouble a01,FXdouble a02,FXdouble a03,
               FXdouble a10,FXdouble a11,FXdouble a12,FXdouble a13,
               FXdouble a20,FXdouble a21,FXdouble a22,FXdouble a23,
               FXdouble a30,FXdouble a31,FXdouble a32,FXdouble a33);

  FXMat4d& set(const FXVec4d& a,const FXVec4d& b,const FXVec4d& c,const FXVec4d& d);

  FXMat4d& set(const FXQuatd& quat);

  FXMat4d& operator+=(const FXMat4d& w);
  FXMat4d& operator-=(const FXMat4d& w);
  FXMat4d& operator*=(const FXMat4d& w);
  FXMat4d& operator*=(FXdouble w);
  FXMat4d& operator/=(FXdouble w);

  FXVec4d& operator[](FXint i){return m[i];}
  const FXVec4d& operator[](FXint i) const {return m[i];}

  operator FXdouble*(){return m[0];}
  operator const FXdouble*() const {return m[0];}

  FXMat4d operator-() const;

  FXMat4d& identity();

  FXbool isIdentity() const;

  FXMat4d& setOrtho(FXdouble xlo,FXdouble xhi,FXdouble ylo,FXdouble yhi,FXdouble zlo,FXdouble zhi);

  void getOrtho(FXdouble& xlo,FXdouble& xhi,FXdouble& ylo,FXdouble& yhi,FXdouble& zlo,FXdouble& zhi) const;

  FXMat4d& setInverseOrtho(FXdouble xlo,FXdouble xhi,FXdouble ylo,FXdouble yhi,FXdouble zlo,FXdouble zhi);

  FXMat4d& setFrustum(FXdouble xlo,FXdouble xhi,FXdouble ylo,FXdouble yhi,FXdouble zlo,FXdouble zhi);

  void getFrustum(FXdouble& xlo,FXdouble& xhi,FXdouble& ylo,FXdouble& yhi,FXdouble& zlo,FXdouble& zhi) const;

  FXMat4d& setInverseFrustum(FXdouble xlo,FXdouble xhi,FXdouble ylo,FXdouble yhi,FXdouble zlo,FXdouble zhi);

  FXMat4d& left();

  FXMat4d& rot(const FXMat3d& r);

  FXMat4d& rot(const FXQuatd& q);

  FXMat4d& rot(const FXVec3d& v,FXdouble c,FXdouble s);

  FXMat4d& rot(const FXVec3d& v,FXdouble phi);

  FXMat4d& xrot(FXdouble c,FXdouble s);
  FXMat4d& xrot(FXdouble phi);

  FXMat4d& yrot(FXdouble c,FXdouble s);
  FXMat4d& yrot(FXdouble phi);

  FXMat4d& zrot(FXdouble c,FXdouble s);
  FXMat4d& zrot(FXdouble phi);

  FXMat4d& look(const FXVec3d& from,const FXVec3d& to,const FXVec3d& up);

  FXMat4d& trans(FXdouble tx,FXdouble ty,FXdouble tz);
  FXMat4d& trans(const FXVec3d& v);

  FXMat4d& scale(FXdouble sx,FXdouble sy,FXdouble sz);
  FXMat4d& scale(const FXVec3d& v);
  FXMat4d& scale(FXdouble s);

  FXMat4d& mirror(const FXVec3d& n);

  FXdouble det() const;

  FXMat4d transpose() const;

  FXMat4d invert() const;

  FXMat4d affineInvert() const;

  FXMat4d rigidInvert() const;

  FXMat3d normalMatrix() const;

 ~FXMat4d(){}
  };


extern FXAPI FXVec3d operator*(const FXMat4d& m,const FXVec3d& v);
extern FXAPI FXVec4d operator*(const FXMat4d& m,const FXVec4d& v);

extern FXAPI FXVec3d operator*(const FXVec3d& v,const FXMat4d& m);
extern FXAPI FXVec4d operator*(const FXVec4d& v,const FXMat4d& m);

extern FXAPI FXMat4d operator+(const FXMat4d& a,const FXMat4d& b);
extern FXAPI FXMat4d operator-(const FXMat4d& a,const FXMat4d& b);

extern FXAPI FXMat4d operator*(const FXMat4d& a,const FXMat4d& b);

extern FXAPI FXMat4d operator*(FXdouble x,const FXMat4d& a);
extern FXAPI FXMat4d operator*(const FXMat4d& a,FXdouble x);
extern FXAPI FXMat4d operator/(const FXMat4d& a,FXdouble x);
extern FXAPI FXMat4d operator/(FXdouble x,const FXMat4d& a);

extern FXAPI FXbool operator==(const FXMat4d& a,const FXMat4d& b);
extern FXAPI FXbool operator!=(const FXMat4d& a,const FXMat4d& b);
extern FXAPI FXbool operator==(const FXMat4d& a,FXdouble n);
extern FXAPI FXbool operator!=(const FXMat4d& a,FXdouble n);
extern FXAPI FXbool operator==(FXdouble n,const FXMat4d& a);
extern FXAPI FXbool operator!=(FXdouble n,const FXMat4d& a);

extern FXAPI FXMat4d orthogonalize(const FXMat4d& m);

extern FXAPI FXStream& operator<<(FXStream& store,const FXMat4d& m);

extern FXAPI FXStream& operator>>(FXStream& store,FXMat4d& m);

}

#endif