ElmNorm.h

Go to the documentation of this file.

// $Id$
//==============================================================================
//==============================================================================
 
#ifndef _ELM_NORM_H
#define _ELM_NORM_H
 
#include "LocalIntegral.h"
#include <cstddef>
 
 
class ElmNorm : public LocalIntegral
{
public:
  ElmNorm(double* p, size_t n) : ptr(p), nnv(n) {}
  explicit ElmNorm(size_t n) : buf(n,0.0), nnv(n) { ptr = &buf.front(); }
  virtual ~ElmNorm() {}
 
  double& operator[](size_t i) { return ptr[i]; }
  const double& operator[](size_t i) const { return ptr[i]; }
 
  size_t size() const { return nnv; }
 
  bool externalStorage() const { return buf.empty(); }
 
  virtual void destruct()
  {
    if (buf.empty())
    {
      vec.clear();
      psol.clear();
    }
    else
      delete this; // The internal buffer has been used, delete the whole thing
  }
 
private:
  RealArray buf; 
  double*   ptr; 
  size_t    nnv; 
 
public:
  Vectors  psol; 
};
 
#endif