some ref files updated (again). 60/11260/1
Cedric Delamarre [Thu, 11 Apr 2013 12:13:01 +0000 (14:13 +0200)]
Change-Id: I6776a102351c282aec66ba2260dd1854bcb5c996

scilab/modules/differential_equations/tests/unit_tests/bvode.dia.ref
scilab/modules/differential_equations/tests/unit_tests/feval.dia.ref
scilab/modules/differential_equations/tests/unit_tests/impl.dia.ref
scilab/modules/differential_equations/tests/unit_tests/int2d.dia.ref

index c4a6c4b..2ce01ed 100644 (file)
-
 // =============================================================================
-
 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-
 // Copyright (C) ????-2008 - INRIA
-
 // Copyright (C) 2009 - DIGITEO - Allan CORNET
-
 //
-
 //  This file is distributed under the same license as the Scilab package.
-
 // =============================================================================
-
-
 // <-- CLI SHELL MODE -->
-
 // to check that bvode works
-
 function df=dfsub(x,z)
-
     df=[0,0,-6/x**2,-6/x]
-
 endfunction
-
-
 function f=fsub(x,z)
-
     f=(1 -6*x**2*z(4)-6*x*z(3))/x**3
-
 endfunction
-
-
 function g=gsub(i,z)
-
     g=[z(1),z(3),z(1),z(3)];g=g(i)
-
 endfunction
-
-
 function dg=dgsub(i,z)
-
     dg=[1,0,0,0;0,0,1,0;1,0,0,0;0,0,1,0];
-
     dg=dg(i,:);
-
 endfunction
-
-
 function [z,mpar]=guess(x)// unused here
-
     z=0;
-
     mpar=0;
-
 endfunction
-
-
  //define trusol for testing purposes
-
 function u=trusol(x)
-
     u=0*ones(4,1)
-
     u(1) =  0.25*(10*log(2)-3)*(1-x) + 0.5 *( 1/x   + (3+x)*log(x) - x)
-
     u(2) = -0.25*(10*log(2)-3)       + 0.5 *(-1/x^2 + (3+x)/x      + log(x) - 1)
-
     u(3) = 0.5*( 2/x^3 + 1/x   - 3/x^2)
-
     u(4) = 0.5*(-6/x^4 - 1/x/x + 6/x^3)
-
 endfunction
-
-
 fixpnt=0;
-
 m=4;
-
 ncomp=1;
-
 aleft=1;
-
 aright=2;
-
 zeta=[1,1,2,2];
-
 ipar=zeros(1,11);
-
 ipar(3)=1;
-
 ipar(4)=2;
-
 ipar(5)=2000;
-
 ipar(6)=200;
-
 ipar(7)=1;
-
 ltol=[1,3];
-
 tol=[1.e-11,1.e-11];
-
 res=aleft:0.1:aright;
-
-
 z=bvode(res,ncomp,m,aleft,aright,zeta,ipar,ltol,tol,fixpnt,fsub,dfsub,gsub,dgsub,guess);
-
 z1=[];for x=res,z1=[z1,trusol(x)]; end;
-
 if norm(z-z1) > 1D8*%eps then bugmes();quit;end;
-
-
 clear gsub
-
 clear guess
-
 clear trusol
-
 exec(SCI+'/modules/differential_equations/tests/unit_tests/bvode_tst.sci');
-
-
 [z,zf]=col1();
-
 if max(abs(z-zf))>1.e-5 then bugmes();quit;end
-
-
 [z,zf]=col2(0);
-
 // Fortran Coded version 
-
 if max(abs(z-zf))>1.e-2 then bugmes();quit;end
-
-
 [z,zf]=col2(1);
-
 // Scilab coded version 
-
 if max(abs(z-zf))>1.e-2 then bugmes();quit;end
-
-
index b3cacbb..8410436 100644 (file)
-
 // =============================================================================
-
 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-
 // Copyright (C) ????-2008 - INRIA
-
 //
-
 //  This file is distributed under the same license as the Scilab package.
-
 // =============================================================================
-
-
 // <-- CLI SHELL MODE -->
-
-
 // Test 1
-
-
 function z=f(x,y,t)
-
     z=x^2+y^2+t;
-
 endfunction
-
-
 res  = feval(1:10,1:5,list(f,10));
-
 res1 = [];
-
 for j=1:5
-
   res1=[res1, ((1:10).*(1:10))'+(j**2+10)*ones(10,1)];
-
 end
-
 if res<>res1 then bugmes();quit;end
-
-
 // Test 2
-
-
 function z=g(x,y)
-
     z=x+%i*y;
-
 endfunction
-
-
 res  =feval(1:10,1:5,g);
-
 res1 =[];
-
 for j=1:5
-
   res1=[res1, (1:10)'+%i*j*ones(10,1)];
-
 end
-
 if res<>res1 then bugmes();quit;end
-
-
-
 // An other example from :
-
 // http://numerics.tinabargs.com/2008/07/scilab-functions-as-input-for-other-scilab-functions/
-
 // Thanks to Ma. Cristina R. Bargo for the authorization to include this test 
-
 // into Scilab
-
 function y = plop(x)
-
   y = 2*x - 1;
-
 endfunction
-
-
 function y = plip(x)
-
   y = x^2 - 5*x + 2;
-
 endfunction
-
-
-
 function y = fcninput(fcnname, x)
-
 // fcnname is the name of the function to be evaluated at x
-
   y = feval(x,fcnname)
-
 endfunction
-
-
 fcninput(plip,1:10)
  ans  =
-
   - 2.  - 4.  - 4.  - 2.    2.    8.    16.    26.    38.    52.
-
-
-
 fcninput(plop,1:10)
  ans  =
-
     1.    3.    5.    7.    9.    11.    13.    15.    17.    19.
-
-
index eb7c242..9c8142b 100644 (file)
@@ -1,92 +1,42 @@
-
 // =============================================================================
-
 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-
 // Copyright (C) ????-2008 - INRIA
-
 //
-
 //  This file is distributed under the same license as the Scilab package.
-
 // =============================================================================
-
-
 // <-- CLI SHELL MODE -->
-
-
 y=impl([1;0;0],[-0.04;0.04;0],0,0.4,'resid','aplusp');
-
 if size(y) <> [3 1] then bugmes();quit;end
-
 if norm( y - [0.9851721;0.0000339;0.0147941]) < %eps then bugmes();quit;end
-
-
 //   DAE's
-
 Leps=1.e-6;
-
 // scilab macros
-
 function r=resid(t,y,s)
-
     r(1)=-.04d0*y(1)+1.d4*y(2)*y(3)-s(1);
-
     r(2)=.04d0*y(1)-1.d4*y(2)*y(3)-3.d7*y(2)*y(2)-s(2);
-
     r(3)=y(1)+y(2)+y(3)-1.d0;
-
 endfunction
-
-
 function p=aplusp(t,y,p)
-
     p(1,1)=p(1,1)+1.d0;
-
     p(2,2)=p(2,2)+1.d0;
-
 endfunction
-
-
 function p=dgbydy(t,y,s)
-
     p(1,1)=-.04d0;
-
     p(1,2)=1.d4*y(3);
-
     p(1,3)=1.d4*y(2);
-
     p(2,1)=.04d0;
-
     p(2,2)=-1.d4*y(3)-6.d7*y(2);
-
     p(2,3)=-1.d4*y(2);
-
     p(3,1)=1.d0;
-
     p(3,2)=1.d0;
-
     p(3,3)=1.d0;
-
 endfunction
-
-
 //         calling scilab
-
 //  
-
 yt=impl([1;0;0],[-.04;.04;0],0,0.4,resid,aplusp,dgbydy);
-
 //  
-
 r1=yt-impl([1;0;0],[-0.04;0.04;0],0,0.4,resid,aplusp);
-
 if abs(r1) > 1.e-10 then bugmes();quit;end
-
 //   calling fortran
-
 r2=yt-impl([1;0;0],[-0.04;0.04;0],0,0.4,'resid','aplusp');
-
 if abs(r2) > 1.e-10 then bugmes();quit;end
-
-
index 75b1d0d..b4f52b2 100644 (file)
@@ -1,34 +1,15 @@
-
 // =============================================================================
-
 // Scilab ( http://www.scilab.org/ ) - This file is part of Scilab
-
 // Copyright (C) ????-2008 - INRIA
-
 //
-
 //  This file is distributed under the same license as the Scilab package.
-
 // =============================================================================
-
-
 // <-- CLI SHELL MODE -->
-
-
 X=[0,0;1,1;1,0];
-
 Y=[0,0;0,1;1,1];
-
 function z=f(x,y)
-
     z=cos(x+y)
-
 endfunction
-
-
 [I,e]=int2d(X,Y,f);
-
 if e <%eps then bugmes();quit;end;
-
 if abs(I-0.5) < %eps then bugmes();quit;end;
-