LSSolver.h

Go to the documentation of this file.

// This file is part of ShapeOp, a lightweight C++ library
// for static and dynamic geometry processing.
//
// Copyright (C) 2014 Sofien Bouaziz <sofien.bouaziz@gmail.com>
// Copyright (C) 2014 LGG EPFL
//
// This Source Code Form is subject to the terms of the Mozilla
// Public License v. 2.0. If a copy of the MPL was not distributed
// with this file, You can obtain one at http://mozilla.org/MPL/2.0/.
#ifndef LSSOLVER_H
#define LSSOLVER_H
#include "Types.h"

#include <iostream>
#include <unsupported/Eigen/src/IterativeSolvers/MINRES.h>

namespace ShapeOp {

class SHAPEOP_API LSSolver {
 public:
  virtual ~LSSolver() {};
  virtual void initialize(const SparseMatrix &A, unsigned int iteration = 1) = 0;
  virtual VectorX solve(const VectorX &b, const VectorX &x0) const = 0;
  virtual Eigen::ComputationInfo info() const = 0;
};

class SHAPEOP_API SimplicialLDLTSolver : public LSSolver {
 public:
  virtual ~SimplicialLDLTSolver() {};
  virtual void initialize(const SparseMatrix &A, unsigned int iteration) override final;
  virtual VectorX solve(const VectorX &b, const VectorX &x0) const override final;
  virtual Eigen::ComputationInfo info() const override final;
 private:
  Eigen::SimplicialLDLT<SparseMatrix> solver_;
};

class SHAPEOP_API CGSolver : public LSSolver {
 public:
  virtual ~CGSolver() {};
  virtual void initialize(const SparseMatrix &A, unsigned int iteration) override final;
  virtual VectorX solve(const VectorX &b, const VectorX &x0) const override final;
  virtual Eigen::ComputationInfo info() const override final;
 private:
  Eigen::ConjugateGradient<SparseMatrix, Eigen::Lower, Eigen::IncompleteLUT<Scalar> > solver_;
};

class SHAPEOP_API MINRESSolver : public LSSolver {
 public:
  virtual ~MINRESSolver() {};
  virtual void initialize(const SparseMatrix &A, unsigned int iteration) override final;
  virtual VectorX solve(const VectorX &b, const VectorX &x0) const override final;
  virtual Eigen::ComputationInfo info() const override final;
 private:
  Eigen::MINRES<SparseMatrix, Eigen::Lower, Eigen::IncompleteLUT<Scalar> > solver_;
};

class SHAPEOP_API SORSolver : public LSSolver {
 public:
  SORSolver(ShapeOp::Scalar relaxation = 1.6);
  virtual ~SORSolver() {};
  virtual void initialize(const SparseMatrix &A, unsigned int iteration) override final;
  virtual VectorX solve(const VectorX &b, const VectorX &x0) const override final;
  virtual Eigen::ComputationInfo info() const override final;
 private:
  SparseMatrixT<Eigen::RowMajor> A_;
  ShapeOp::Scalar relaxation_;
  unsigned int iteration_;
};
} // namespace ShapeOp
#ifdef SHAPEOP_HEADER_ONLY
#include "LSSolver.cpp"
#endif
#endif // LSSOLVER_H

---

Generated on Thu Apr 2 2015 02:37:05 for ShapeOp by   1.8.9.1