desystestspring2.h

Go to the documentation of this file.

#ifndef DESYSTESTSPRING2_H
#define DESYSTESTSPRING2_H

#include <print.h>
#include <typedefs.h>

class desystestspring2
{
public:

  doublec * xval;
  doublec * wi;
  double ki[3];
  uint current;

  desystestspring2(doublec * xval_, doublec * wi_)
    : xval(xval_), wi(wi_), current(0) {}
//  ~desystestspring2()
//    { delete[] ki; }

  void kset(doublec *k);

  void equ01(double &  Dz, doublec z, doublec x ) const;
  void equ02(double &  Dz, doublec z, doublec x ) const;

  void operator ++ ()
    { ++current; current %= 2; }
  void operator () (double & Dz, doublec z, doublec x ) const
  { 
    switch (current)
    {
      case 0: equ01(Dz,z,x); break;
      case 1: equ02(Dz,z,x); break;
    }
  }
};

#endif