Parent Directory | Revision Log

Revision **468** -
(**show annotations**)
(**download**)
(**as text**)

*Mon Apr 17 03:20:06 2006 UTC*
(18 years, 2 months ago)
by *ben.allan*

File MIME type: text/x-ascend

File size: 2026 byte(s)

File MIME type: text/x-ascend

File size: 2026 byte(s)

new file.

1 | (*********************************************************************\ |

2 | sensitivity_test.asc |

3 | by Ben Allan |

4 | Part of the Ascend Library |

5 | |

6 | This file is part of the Ascend modeling library. |

7 | |

8 | Copyright (C) 1996 |

9 | |

10 | The Ascend modeling library is free software; you can redistribute |

11 | it and/or modify it under the terms of the GNU General Public License as |

12 | published by the Free Software Foundation; either version 2 of the |

13 | License, or (at your option) any later version. |

14 | |

15 | The Ascend Language Interpreter is distributed in hope that it will be |

16 | useful, but WITHOUT ANY WARRANTY; without even the implied warranty of |

17 | MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |

18 | General Public License for more details. |

19 | |

20 | You should have received a copy of the GNU General Public License along with |

21 | the program; if not, write to the Free Software Foundation, Inc., 675 |

22 | Mass Ave, Cambridge, MA 02139 USA. Check the file named COPYING. |

23 | |

24 | This module is used to calculate the sensitity matrix dX/dU where |

25 | X is a vector OF vapor compositions and U is a vector OF liquid |

26 | compositions. This sensitivity matrix is THEN used to calculate |

27 | J, the Jacobian, as described IN Fidkowski et al. (1995). |

28 | |

29 | \*********************************************************************) |

30 | |

31 | REQUIRE "system.a4l"; |

32 | |

33 | IMPORT PackSolve_Init FROM "libsensitivity"; |

34 | |

35 | MODEL sensitivity_test; |

36 | nc IS_A integer_constant; |

37 | X[1..nc] IS_A solver_var; |

38 | U[1..nc], Unew[1..nc] IS_A solver_var; |

39 | dx_du[1..nc][1..nc] IS_A solver_var; |

40 | (* system *) |

41 | k,x,y,u IS_A solver_var; |

42 | k*x - y = 7; |

43 | y = u^2; |

44 | (* dx/du = 2*u/k *) |

45 | nc:==1; |

46 | u,U[1] ARE_THE_SAME; |

47 | x,X[1] ARE_THE_SAME; |

48 | METHODS |

49 | METHOD values; |

50 | END values; |

51 | METHOD clear; |

52 | END clear; |

53 | METHOD specify; |

54 | u.fixed := TRUE; |

55 | k.fixed := TRUE; |

56 | END specify; |

57 | METHOD reset; |

58 | RUN clear; |

59 | RUN specify; |

60 | END reset; |

61 | METHOD analyze; |

62 | EXTERNAL do_sensitivity(SELF,U[1..nc],X[1..nc],dx_du[1..nc][1..nc]); |

63 | END analyze; |

64 | END sensitivity_test; |

john.pye@anu.edu.au | ViewVC Help |

Powered by ViewVC 1.1.22 |