CharacteristicFlux.h

#ifndef CHARACTERISTICFLUX_H
#define CHARACTERISTICFLUX_H
#include "TNT2.h"
#include "globals.h"
#include "ConfigParams.h"

using namespace TNT;
using namespace JAMA;

//du/dt + A du/dx +Bu=0

//dw/dt + Lamb dw/dx + ?=0 characteristic form

class CharacteristicFlux
{
 public:
  vector<double> AV; //see above
  vector<double> AtrimmedV; //part of A that has nonzero derivative dependence
  vector<double> SmatrixV; //matrix with eigenvectors of A in columns
  vector<double> SinvV; //inverse of S matrix
  vector<double> LambdaV; //eigenvalue of A matrix (characteristic form)
  //  Array2D<double> trimA(); 
  int Adimension; //dimension of A
  int Ddimension; //dimension of non-zero derivative variables

 public:
  CharacteristicFlux(vector<double> Amatrix, vector<double> Atrimmed);
  vector<double> getA();
  vector<double> getS();
  vector<double> getSinv();
  vector<double> getLambda();
  vector<double> getAtrimmed();
  int getAdim();//untrimmed A dimension
  int getDdim();//derivative dimension (Atrimmed.dim())
};

#endif