VectorSpaceBasis.h

// Copyright (C) 2013-2019 Patrick E. Farrell and Garth N. Wells
//
// This file is part of DOLFINx (https://www.fenicsproject.org)
//
// SPDX-License-Identifier:    LGPL-3.0-or-later
 
#pragma once
 
#include <memory>
#include <petscmat.h>
#include <vector>
 
namespace dolfinx::la
{
class PETScVector;
 
 
class VectorSpaceBasis
{
public:
  VectorSpaceBasis(const std::vector<std::shared_ptr<PETScVector>>& basis);
 
  VectorSpaceBasis(const VectorSpaceBasis& basis) = delete;
 
  VectorSpaceBasis(VectorSpaceBasis&& basis) = default;
 
  ~VectorSpaceBasis() = default;
 
  void orthonormalize(double tol = 1.0e-10);
 
  bool is_orthonormal(double tol = 1.0e-10) const;
 
  bool is_orthogonal(double tol = 1.0e-10) const;
 
  bool in_nullspace(const Mat A, double tol = 1.0e-10) const;
 
  void orthogonalize(PETScVector& x) const;
 
  int dim() const;
 
  std::shared_ptr<const PETScVector> operator[](int i) const;
 
private:
  const std::vector<std::shared_ptr<PETScVector>> _basis;
};
} // namespace dolfinx::la