[Octave-cvsupdate] SF.net SVN: octave:[10612] trunk/octave-forge/extra/secs1d/inst/ secs1d_dd_newton.m

[<cd...@us...>]

Revision: 10612
          http://octave.svn.sourceforge.net/octave/?rev=10612&view=rev
Author:   cdf
Date:     2012-06-11 07:29:39 +0000 (Mon, 11 Jun 2012)
Log Message:
-----------
bug fix in newton solver

Modified Paths:
--------------
    trunk/octave-forge/extra/secs1d/inst/secs1d_dd_newton.m

Modified: trunk/octave-forge/extra/secs1d/inst/secs1d_dd_newton.m
===================================================================
--- trunk/octave-forge/extra/secs1d/inst/secs1d_dd_newton.m	2012-06-09 15:38:32 UTC (rev 10611)
+++ trunk/octave-forge/extra/secs1d/inst/secs1d_dd_newton.m	2012-06-11 07:29:39 UTC (rev 10612)
@@ -128,12 +128,12 @@
   A21 = - bim1a_laplacian (x, un .* nm, 1);
   A22 = bim1a_advection_diffusion (x, un, 1, 1, V) + bim1a_reaction (x, 1, p .* fact);
   A23 = bim1a_reaction (x, 1, n .* fact);
-  R2  = A22 * n + bim1a_rhs (x, 1, (p .* n - theta .^ 2) .* fact);
+  R2  = A22 * n + bim1a_rhs (x, 1, (- theta .^ 2) .* fact);
 
   A31 = bim1a_laplacian (x, up .* pm, 1);
   A32 = bim1a_reaction (x, 1, p .* fact);
   A33 = bim1a_advection_diffusion (x, up, 1, 1, -V) + bim1a_reaction (x, 1, n .* fact);
-  R3  = A33 * p + bim1a_rhs (x, 1, (p .* n - theta .^ 2) .* fact);
+  R3  = A33 * p + bim1a_rhs (x, 1, (- theta .^ 2) .* fact);
 
   res = [R1(2:end-1); R2(2:end-1); R3(2:end-1)];
 

This was sent by the SourceForge.net collaborative development platform, the world's largest Open Source development site.