Changeset aa7db24 in tol for tol_tests

Timestamp:

Apr 9, 2010, 8:20:35 AM (17 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

1941-docker-image-to-test, devel, master

Children:

4db909c

Parents:

8a6d7f1

Message:

Introducción del profiler para depurar la eficiencia del algoritmo

git-svn-id: ​https://www.tol-project.org/svn/tolp/trunk@2085 be4846ae-f051-483a-831f-b64e7d08527f

File:

• tol_tests/tol/ARIMA/mcmc/mtrymh.tol (modified) (2 diffs)

tol_tests/tol/ARIMA/mcmc/mtrymh.tol

@@ -28 +28 @@
     Set  PeriodsTable     = [[ 1,12 ]]
   ]];
-  Set Include("./SimulArima/SimulArima.tol")
+  Set Include("../SimulArima/SimulArima.tol")
 };
   
@@ -634 +634 @@
 Real numSim = 0;
 
+Real append_draw(VMatrix draw)
+{
+  If(!(numSim%10), WriteLn("
+##############################################################################
+Simulation "<<numSim+ " _.acceptRatio="<<mtm_arima::_.acceptRatio));
+  Matrix mcmc := mcmc << VMat2Mat(Tra(draw));
+  Real numSim := numSim+1
+};
+
 /* */ 
 Real mtm_arima::_.sFactor := 1.01;
 Real numSim := 0;
-While(numSim<1000,
-{
-  If(!(numSim%10), WriteLn("
-##############################################################################
-Simulation MH "<<numSim+ " _.acceptRatio="<<mtm_arima::_.acceptRatio));
-  Matrix mcmc := mcmc << VMat2Mat(Tra(mtm_arima::draw_single_try(0)));
-  Real numSim := numSim+1
-});
+Real mh.t0 = Copy(Time);
+While(numSim<1000, append_draw(mtm_arima::draw_single_try(0)));
+Real mh.time = (Copy(Time) - mh.t0)/numSim;
+WriteLn("Time by MH simulation "<<mh.time);
+
 /* */
 Real mtm_arima::_.sFactor := 1.0;
 Real numSim := 0;
-While(numSim<2500,
-{
-  If(!(numSim%10), WriteLn("
-##############################################################################
-Simulation MTM "<<numSim+ " _.acceptRatio="<<mtm_arima::_.acceptRatio));
-  Matrix mcmc := mcmc << VMat2Mat(Tra(mtm_arima::draw_multiple_try(0)));
-  Real numSim := numSim+1
-});
+Real mtm.t0 = Copy(Time);
+Real TolOprProfiler.Enabled := True;
+While(numSim<2500, append_draw(mtm_arima::draw_multiple_try(0)));
+Real TolOprProfiler.Dump("draw_multiple_try.csv");
+Real mtm.time = (Copy(Time) - mtm.t0)/numSim;
+WriteLn("Time by MTM simulation "<<mtm.time);
+/* */
 
 Matrix real_param = Group("ConcatColumns", For(1,Columns(PARAM),Matrix(Real j)