PLUTO
jacobian.c
Go to the documentation of this file.
1 #include "pluto.h"
2 
3 /* ******************************************************** */
4 void Jacobian (double *v, double *rhs, double **dfdy)
5 /*
6  *
7  * Compute the jacobian J(k,l) = dfdy
8  *
9  * k = row index
10  * l = column index
11  *
12  * J(0,0) J(0,1) ... J(0, n-1)
13  * J(1,0) J(1,1) .... J(1, n-1)
14  * . . .
15  * . . .
16  * . . .
17  * J(n-1,0) .... J(n-1, n-1)
18  *
19  *
20  * or,
21  *
22  * +-----------------------+
23  * + | |
24  * + | |
25  * + | |
26  * + dX'/dX | dX'/dp |
27  * + (JXX) | (JXp) |
28  * + | |
29  * + | |
30  * +--------------+--------+
31  * + dp'/dX | dp'/dp |
32  * + (JpX) | Jpp |
33  * +-----------------------+
34  *
35  *
36  *
37  ********************************************************** */
38 {
39  print (" ! Jacobian not defined \n");
40  QUIT_PLUTO(1);
41 }
42 
Jacobian
void Jacobian(real *v, real *rhs, real **dfdy)
Definition: jacobian.c:4
print
void print(const char *fmt,...)
Definition: amrPluto.cpp:497
pluto.h
PLUTO main header file.
QUIT_PLUTO
#define QUIT_PLUTO(e_code)
Definition: macros.h:125