jacobian.h

Go to the documentation of this file.

#ifndef JACOBIANH
#define JACOBIANH

#include "variables.h"
#include "equations.h"

#include <stdio.h>

typedef struct {
  unsigned int neqs;  
  unsigned int nvars; 
  Tequations *J;      
} TJacobian;

void InitJacobian(Tvariables *vs,Tequations *eqs,TJacobian *j);

void CopyJacobian(TJacobian *j_dst,TJacobian *j_src);

void GetJacobianSize(unsigned int *nr,unsigned int *nc,TJacobian *j);

Tequations *GetJacobianColumn(unsigned int nv,TJacobian *J);

Tequation *GetJacobianEquation(unsigned int r,unsigned int c,TJacobian *j);

void AllocateJacobianEvaluation(double ***m,TJacobian *j);

void EvaluateJacobian(double *v,double **m,TJacobian *j);

void EvaluateJacobianInVector(double *v,
                              unsigned int nr,unsigned int nc,double *m,TJacobian *j);

void EvaluateJacobianSubSetInVector(double *v,boolean *sr,
                                    unsigned int nr,unsigned int nc,double *m,TJacobian *j);

void EvaluateTransposedJacobianInVector(double *v,
                                        unsigned int nr,unsigned int nc,double *m,TJacobian *j);

void EvaluateTransposedJacobianSubSetInVector(double *v,boolean *sr,
                                              unsigned int nr,unsigned int nc,double *m,TJacobian *j);

void PrintJacobianEvaluation(FILE *f,double **m,TJacobian *j);

void FreeJacobianEvaluation(double **m,TJacobian *j);

void EvaluateJacobianXVectors(double *p,
                              unsigned int n,unsigned int ng,unsigned int *g,double *v,
                              unsigned int *nr,unsigned int *nc,double ***m,
                              TJacobian *j);

void PrintJacobian(FILE *f,char **varNames,TJacobian *j);

void DeleteJacobian(TJacobian *j);

#endif