SpecSolNID.h

/* SpecSolNID header file
 * Generic class to deal with species involving homogenious speciation
 * must include speciation routine such as speciateNID.cc
 *  $Id: SpecSolNID.h,v 1.17 2006/08/29 21:47:35 kress Exp $
 */

#ifndef SPECSOLNID_H
#define SPECSOLNID_H

#include "SpecSolution.h"
#include "PhaseExceptions.h"

class SpecSolNID: public SpecSolution {
 public:
  SpecSolNID();
  virtual ~SpecSolNID();

  // getting methods //////////////////////////////////////////
  virtual double getGibbs();    
  virtual double getSpecMu(int ispec);
  virtual double getEnthalpy(); 
  double getEntropy();  

 protected:
  double **Wh,**Ws;                                  // nspec x nspec **

  // protected utility methods ///////////////////////////////
  virtual void init();
  virtual int speciate_function(double tol) throw(PhaseError *) ;
  virtual double mufunc(int p,double *x);    //implements speciateNID form
  virtual double dmufunc(int m,int p,double *x);//implements speciateNID form
  virtual double muxs(int ispec,double *x);
  virtual double dmuxs(int m,int p,double *x);
  virtual void getd2Gds2(double **d2gds2);
  virtual void getd3Gds2dT(double **d3gds2dt);
  virtual void getd3Gds3(double ***d3gds3);
 private:
  double *grad,*delxi,**H,*delta,*test,*b,*d,*smalldelta;
};

#endif