jacobian.h

Go to the documentation of this file.

#ifndef JACOBIANH
#define JACOBIANH
 
#include "variables.h"
#include "equations.h"
#include "csmapping.h"
 
#include <stdio.h>
 
typedef struct {
  Tvariables *vars;   
  Tequations *eqs;    
  unsigned int neqs;  
  unsigned int nvars; 
  Tequations *J;      
} TJacobian;
 
void InitJacobian(Tvariables *vs,Tequations *eqs,TJacobian *j);
 
void InitPositionJacobian(Tvariables *vs,Tequations *eqs,TJacobian *j);
 
void InitScalarPositionJacobian(Tvariables *vs,Tequations *eqs,TJacobian *j);
 
void CopyJacobian(TJacobian *j_dst,TJacobian *j_src);
 
void CopyNonDynamicsJacobian(TJacobian *j_dstP,TJacobian *j_dstO,
                             TJacobian *j_src);
 
void CopyPositionJacobian(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 RewriteJacobian2Orig(Tmapping *map,TJacobian *j);
 
void JacobianLinearCombination(boolean *selectedVars,unsigned int *nv,
                               Tequations *eqs,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