ConstrainedLS.H

Go to the documentation of this file.

#ifdef CH_LANG_CC
/*
*      _______              __
*     / ___/ /  ___  __ _  / /  ___
*    / /__/ _ \/ _ \/  V \/ _ \/ _ \
*    \___/_//_/\___/_/_/_/_.__/\___/
*    Please refer to Copyright.txt, in Chombo's root directory.
*/
#endif

#ifndef _CONSTRAINEDLS_H_
#define _CONSTRAINEDLS_H_

#include "Vector.H"
#include "REAL.H"
#include "RealVect.H"
#include "NamespaceHeader.H"

class ConstrainedLS
{
public:
  void LeastSquares(Real**A,
                    Vector<Real>&x,
                    const Vector<Real>&rhs,
                    Real**constraintMat,
                    Real*constraintVec,
                    int nInequality);


  void LeastSquares(Real** A,
                    Vector<Real>&x,
                    const Vector<Real>&rhs);


  void transpose(Real** a_A,
                 Real ** a_Atrans,
                 const int& a_numRows,
                 const int& a_numCols);

  void matMul(Real** a_A,
              Real** a_B,
              Real** a_C,
              const int& a_numRowsA,
              const int& a_numColsA,
              const int& a_numColsB);


  void AtimesX(Real** A,
               const Vector<Real>&x,
               const int& numRowsA,
               Vector<Real>& Ax);

  void allocArray(const int& rows,
                  const int& cols,
                  Real**& A);

  void freeArray(const int& rows,
                 const int& cols,
                 Real**& A);

  void output(const int& rows,
              const int& cols,
              Real**& A,
              char* name);
};

#include "NamespaceFooter.H"
#endif

---