cpu/VecOp.h

#ifndef RPC_VEC_OP_H
#define RPC_VEC_OP_H
 
#include <util/containers/Array.h>
 
/*
* PSCF - Polymer Self-Consistent Field 
*
* Copyright 2015 - 2025, The Regents of the University of Minnesota
* Distributed under the terms of the GNU General Public License.
*/
 
// Forward declaration
namespace Util {
   template <typename T> class Array;
}
 
namespace Pscf {
namespace Prdc {
namespace Cpu {
 
   using namespace Util;

   namespace VecOp {
   
      // Assignment
      
      void eqV(Array<double>& a, Array<double> const & b);
      
      void eqS(Array<double>& a, double b);
      
      // Addition
      
      void addVV(Array<double>& a, Array<double> const & b, 
                 Array<double> const & c);
      
      void addVS(Array<double>& a, Array<double> const & b, double c);
      
      // Subtraction
      
      void subVV(Array<double>& a, 
                 Array<double> const & b, Array<double> const & c);
      
      void subVS(Array<double>& a, Array<double> const & b, double c);
      
      // Multiplication
      
      void mulVV(Array<double>& a, 
                 Array<double> const & b, Array<double> const & c);
      
      void mulVS(Array<double>& a, Array<double> const & b, double c);
      
      // Division
      
      void divVV(Array<double>& a, 
                 Array<double> const & b, Array<double> const & c);
      
      void divVS(Array<double>& a, Array<double> const & b, double c);
      
      void divSV(Array<double>& a, double b, Array<double> const & c);
      
      // Exponentiation
      
      void expV(Array<double>& a, Array<double> const & b);
      
      
      // Compound addition
      
      void addEqV(Array<double>& a, Array<double> const & b);
      
      void addEqS(Array<double>& a, double b);
      
      // Compound subtraction
      
      void subEqV(Array<double>& a, Array<double> const & b);
      
      void subEqS(Array<double>& a, double b);
      
      // Compound multiplication
      
      void mulEqV(Array<double>& a, Array<double> const & b);
      
      void mulEqS(Array<double>& a, double b);
      
      // Compound division
      
      void divEqV(Array<double>& a, Array<double> const & b);
      
      void divEqS(Array<double>& a, double b);
   
 
   } // namespace VecOp
 
} // namespace Cpu
} // namespace Prdc
} // namespace Pscf
 
#endif