Mercurial > octave-nkf
view examples/oregonator.cc @ 18496:80b8873c90ba
Fix pinv bug with 0 values introduced in cset 16b0cd465ecd.
* pinv.cc (Fpinv): Add %!tests for diagonal matrices
* CDiagMatrix.cc (pseudo_inverse), dDiagMatrix.cc (pseudo_inverse),
fCDiagMatrix.cc (pseudo_inverse), fDiagMatrix.cc (pseudo_inverse):
Don't skip assigning 0 to retval element when matrix element is 0 since
retval memory is not initialized to 0.
author | Rik <rik@octave.org> |
---|---|
date | Thu, 20 Feb 2014 09:26:45 -0800 |
parents | 224e76250443 |
children |
line wrap: on
line source
#include <octave/oct.h> DEFUN_DLD (oregonator, args, , "The `oregonator'.\n\ \n\ Reference:\n\ \n\ Oscillations in chemical systems. IV. Limit cycle behavior in a\n\ model of a real chemical reaction. Richard J. Field and Richard\n\ M. Noyes, The Journal of Chemical Physics, Volume 60 Number 5,\n\ March 1974.") { ColumnVector dx (3); ColumnVector x (args(0).vector_value ()); dx(0) = 77.27 * (x(1) - x(0)*x(1) + x(0) - 8.375e-06*pow (x(0), 2.0)); dx(1) = (x(2) - x(0)*x(1) - x(1)) / 77.27; dx(2) = 0.161*(x(0) - x(2)); return octave_value (dx); }