PetscVector.hpp

Go to the documentation of this file.

/*
 * ------------------------------------------------------------------------------------------------------------
 * SPDX-License-Identifier: LGPL-2.1-only
 *
 * Copyright (c) 2016-2024 Lawrence Livermore National Security LLC
 * Copyright (c) 2018-2024 Total, S.A
 * Copyright (c) 2018-2024 The Board of Trustees of the Leland Stanford Junior University
 * Copyright (c) 2023-2024 Chevron
 * Copyright (c) 2019-     GEOS/GEOSX Contributors
 * All rights reserved
 *
 * See top level LICENSE, COPYRIGHT, CONTRIBUTORS, NOTICE, and ACKNOWLEDGEMENTS files for details.
 * ------------------------------------------------------------------------------------------------------------
 */
 
#ifndef GEOS_LINEARALGEBRA_INTERFACES_PETSCVECTOR_HPP_
#define GEOS_LINEARALGEBRA_INTERFACES_PETSCVECTOR_HPP_
 
#include "linearAlgebra/interfaces/VectorBase.hpp"
 
 
extern "C" struct _p_Vec;
 
 
namespace geos
{
 
class PetscVector final : private VectorBase< PetscVector >
{
public:
 
 
  PetscVector();
 
  PetscVector( PetscVector const & src );
 
  PetscVector( PetscVector && src ) noexcept;
 
  PetscVector & operator=( PetscVector const & src );
 
  PetscVector & operator=( PetscVector && src ) noexcept;
 
  ~PetscVector();
 
 
 
  using VectorBase::setName;
  using VectorBase::closed;
  using VectorBase::ready;
  using VectorBase::open;
  using VectorBase::zero;
  using VectorBase::values;
 
  virtual bool created() const override;
 
  virtual void create( localIndex const localSize,
                       MPI_Comm const & comm ) override;
 
  virtual void close() override;
 
  virtual void touch() override;
 
  virtual void reset() override;
 
  virtual void set( real64 const value ) override;
 
  virtual void rand( unsigned const seed ) override;
 
  virtual void scale( real64 const scalingFactor ) override;
 
  virtual void reciprocal() override;
 
  virtual real64 dot( PetscVector const & vec ) const override;
 
  virtual void copy( PetscVector const & x ) override;
 
  virtual void axpy( real64 const alpha,
                     PetscVector const & x ) override;
 
  virtual void axpby( real64 const alpha,
                      PetscVector const & x,
                      real64 const beta ) override;
 
  virtual void pointwiseProduct( PetscVector const & x,
                                 PetscVector & y ) const override;
 
  virtual real64 norm1() const override;
 
  virtual real64 norm2() const override;
 
  virtual real64 normInf() const override;
 
  virtual globalIndex globalSize() const override;
 
  virtual localIndex localSize() const override;
 
  virtual globalIndex ilower() const override;
 
  virtual globalIndex iupper() const override;
 
  virtual void print( std::ostream & os = std::cout ) const override;
 
  virtual void write( string const & filename,
                      LAIOutputFormat const format = LAIOutputFormat::MATRIX_MARKET ) const override;
 
  virtual MPI_Comm comm() const override;
 
 
private:
 
  using Vec = struct _p_Vec *;
 
public:
 
  const Vec & unwrapped() const;
 
  Vec & unwrapped();
 
private:
 
  Vec m_vec;
};
 
} // end geos namespace
 
#endif /*GEOS_LINEARALGEBRA_INTERFACES_PETSCVECTOR_HPP_*/