desystestspring2.cpp

Go to the documentation of this file.

#include <desystestspring2.h>

/*
desystestspring2::desystestspring2()
{
  wi[0] = wi[1] = wi[2] = 1.0;
  // Some default k.
  ki[0] = ki[1] = ki[2] = 0.1;
}
*/


void desystestspring2::kset(doublec * k)
{
  ki[0] = k[0];
  ki[1] = k[1];
  ki[2] = k[2];
}

void desystestspring2::equ01
(
  double & Dz, 
  doublec z, 
  doublec x 
) const
{
  // (wi[0],wi[1],wi[2],wi[3]) corresponds with (w0,w0',w1,w1')
  // w0'' = k1*w1 - k1*w0 - k0*w0
  Dz = ki[1]*wi[2]-ki[1]*wi[0]-ki[0]*wi[0];
}

void desystestspring2::equ02
(
  double &  Dz, 
  doublec z, 
  doublec x 
) const
{
  // (wi[0],wi[1],wi[2],wi[3]) corresponds with (w0,w0',w1,w1')
  // w1'' = -k2*w1 - k1*w1 + k1*w0
  Dz = -ki[2]*wi[2]-ki[1]*wi[2]+ki[1]*wi[0];
}