partialsolve.h

Go to the documentation of this file.

// ========================================================================= //
// This file is part of MyraMath, copyright (c) 2014-2019 by Ryan A Chilton  //
// and distributed by MyraCore, LLC. See LICENSE.txt for license terms.      //
// ========================================================================= //

#ifndef MYRAMATH_MULTIFRONTAL_LU_PARTIALSOLVE_H
#define MYRAMATH_MULTIFRONTAL_LU_PARTIALSOLVE_H

#include <myramath/utility/Number.h>
#include <myramath/jobgraph/JobGraph.h>
#include <myramath/multifrontal/Options.h>
#include <vector>

namespace myra {

// Forward declarations.
class intCRange;
template<class Number> class Matrix;
template<class Number> class CMatrixRange;
template<class Number> class  MatrixRange;

namespace multifrontal {
namespace detail {
namespace lu {

// Forward declarations.
template<class Kernel> class LUContainer; 

} } } // namespace

namespace multifrontal {
namespace lu {
namespace partialsolve {

// Useful typedef.
typedef ::myra::multifrontal::Options Options;

// Implementation detail of partialsolve()
template<class Kernel> JobGraph partialsolve_jobgraph(
  const ::myra::multifrontal::detail::lu::LUContainer<Kernel>& LU,
  const intCRange& i,
  const intCRange& j,
  const CMatrixRange<typename ReflectNumber<Kernel>::type>& B,
  const  MatrixRange<typename ReflectNumber<Kernel>::type>& X,
  char side, char order, char op_L, char up_U, Options options);

// Implementation detail of partialsolve()
template<class Kernel> void partialsolve_inplace(
  const ::myra::multifrontal::detail::lu::LUContainer<Kernel>& LU,
  const intCRange& i,
  const intCRange& j,
  const CMatrixRange<typename ReflectNumber<Kernel>::type>& B,
  const  MatrixRange<typename ReflectNumber<Kernel>::type>& X,
  char side, char order, char op_L, char up_U, Options options);

} } } } // namespace

#endif