-- test routine for interface to chapter d01
indexing
	description: " Test routines for d01"
	legal: "See notice at end of class."
	status: "See notice at end of class.";
	date: "$Date$";
	revision: "$Revision$"

class

	DRIVE
inherit
	EIFFELMATH_TESTING_FRAMEWORK

create
	make

feature -- The main program

	make (args: ARRAY [STRING]) is
			-- test of various routines from d01
		do
			test_root;
			test_sys;
		end;


feature -- Test of c05adc, root finder

	testf: TEST_CLASS;
			-- A class that contains feature f1 such that f1(1.5)=0.

	test_root is
			-- test zero in an interval
		local
			gr: FUNCTION_ROOT_FINDER;
			-- An zero finder for finite intervals.
		do
	 		-- create an instance of the target class for passing to the root finder.
			create testf

			create gr.make;
			gr.make_target (testf, "f1");
			gr.set_interval(1., 2.);

			print("Solving for zero of f1. Answer should be 1.5%N");
			gr.solve;

			print ("Root: ") ;
			print (gr.last_root);
			print_nl ("");
		end;

	test_sys is
			-- test zero of non-linear system with and without Jacobian
		local
			gr: NONLINEAR_SYSTEM_SOLVER;
			t: TEST_SYSTEM;
			x: ARRAY [DOUBLE];
		do
	 		-- create an instance of the target class for passing to the root finder.
			create t;

			create gr.make (t);

			print("Solving for zero of system. Answer should be %N");
			x := <<-.5707, -.6816, -.7017, -.7042, -.7014,
				-.6919, -.6658, -.5960, -.4164>>;
			print("Checking function and Jacobian from Eiffel%N")
			print(x.tagged_out);
			print(t.value(x).tagged_out);
			print(t.jacobian(x).tagged_out);

			x := << -1., -1., -1., -1., -1., -1., -1., -1., -1.>>;
			gr.do_not_use_jacobian;
			gr.solve (x);
			
			print ("Solution without Jacobian, last_root: ") ;
			print (gr.last_root.tagged_out);
			print ("last_value: ") ;
			print (gr.last_value.tagged_out);
			print_nl ("");

			gr.use_jacobian;
			gr.check_jacobian (x);
			print("Jacobian checked against system, ok.%N");
			gr.solve (x);

			print ("Solution with Jacobian, last_root: ") ;
			print (gr.last_root.tagged_out);
			print ("last_value: ") ;
			print (gr.last_value.tagged_out);
			print_nl ("");
		end;

indexing
	copyright:	"Copyright (c) 1984-2006, Eiffel Software and others"
	license:	"Eiffel Forum License v2 (see http://www.eiffel.com/licensing/forum.txt)"
	source: "[
			 Eiffel Software
			 356 Storke Road, Goleta, CA 93117 USA
			 Telephone 805-685-1006, Fax 805-685-6869
			 Website http://www.eiffel.com
			 Customer support http://support.eiffel.com
		]"


end -- class DRIVE