Changeset 1e064ff in tol

Timestamp:

Dec 30, 2009, 6:22:11 PM (17 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

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

Children:

4a13114

Parents:

7874ccc

Message:

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

File:

• tol_pkg/MulVarFunAprox/test/test_0001/test.tol (modified) (7 diffs)

tol_pkg/MulVarFunAprox/test/test_0001/test.tol

@@ -29 +29 @@
   {
     Real min = Rand(2,4);
-    Real len = Rand(5,10)*((i/n)^(1/n));
+    Real len = Rand(5,15)*((i/n)^(1/n));
     Rand(J,1,min,min+len)
   }))
 };
-VMatrix X = generate.nodal.points(15000,2);
+VMatrix X = generate.nodal.points(5000,2);
 /* * /
 MulVarFunAprox::@DesignC2 design = {
@@ -47 +47 @@
 /* */
 Real design::max.deg := 6;
-Real design::use.compact.base := True;
+Real design::use.compact.base := False;
 Real design::min.contrast.surface := 2.0;
 
@@ -122 +122 @@
 };
 
-Real use_extra = False;
+Real use_extra = True;
 
 Code fun.eval   = If(use_extra,  fun.eval.extra.base,  fun.eval.intra.base);
@@ -133 +133 @@
 VMatrix F_X= fun.eval(design::_.X);
 
-Real t0 := Copy(Time);
-WriteLn("Running regression model for design matrix "<<design::_.Q);
+Real box.cox.exp = 1;
 MulVarFunAprox::@ModelC2 model = MulVarFunAprox::@ModelC2::New(
    design,
-   F_X);
+   RPow(F_X,box.cox.exp));
+WriteLn("Running regression model for design matrix "<<design::_.Q);
+Real t0 := Copy(Time);
+Real model::estimate(0);
 Real t.estimating = Copy(Time)-t0;
 WriteLn("Estimated model in "<<t.estimating+" seconds");
-
-
+/*
 VMatrix dF_X1 = der.1.eval(design::_.X);
 VMatrix dF_X2 = der.2.eval(design::_.X);
@@ -159 +160 @@
 Real dF_X2_err.ste = VMatStDs(dF_X2_err);
 
-
-VMatrix X2 = { design::generate.internal.points(
-  Real m = design::_.J/2,
-  True)};
+*/
+VMatrix X2 = { design::generate.internal.points(Real m = design::_.J/2)};
 
 Real b::upgrade.cache(X2,0,0);
@@ -168 +167 @@
 Real t0 := Copy(Time);
 WriteLn("Evaluating forecast for "<<VRows(X2)+" points...");
-VMatrix F2.frc = model::evaluate(X2);
+VMatrix F2.frc = model::evaluate(X2,0);
 Real t.forecasting = Copy(Time)-t0;
 WriteLn("Evaluated forecast in "<<t.forecasting+" seconds");
@@ -205 +204 @@
 
 VMatrix Y.n.m = model::_.design::_.Yn | model::_.design::_.Ym;
+
+VMatrix xyz = X|F_X;
+
+Matrix MatWriteFile("real.xyz.bbm",VMat2Mat(xyz));
+
+VMatrix _.X = X;
+Real _.n = VColumns(_.X);
+    Set _.X.col.info = For(1,_.n,Set(Real i)
+    {[[
+      VMatrix col = SubCol(_.X,[[i]]);
+      Real min = VMatMin(col);
+      Real max = VMatMax(col);
+      Real length = max-min;
+      Set Range(min,max,length/99)
+    ]]});
+VMatrix grid.X = Mat2VMat(SetMat(Group("CartProd",Traspose(Extract(_.X.col.info,5))[1])));
+VMatrix grid.F = RPow(model::evaluate(grid.X,?),1/box.cox.exp);
+
+//Rango de latitudes
+Matrix x = SetCol(_.X.col.info[1][5]);
+//Rango de latitudes
+Matrix y = SetCol(_.X.col.info[2][5]);
+Matrix z = 
+{
+  Matrix aux= Group("ConcatColumns",For(1,Rows(x),Matrix(Real j)
+  {
+    VMat2Mat(SubRow(grid.F, Range(1+(j-1)*Rows(y),(j)*Rows(y),1)))
+  }));
+  Tra(aux)
+};
+
+Matrix MatWriteFile("aprox.x.bbm",x);
+Matrix MatWriteFile("aprox.y.bbm",y);
+Matrix MatWriteFile("aprox.z.bbm",z);
+
+Set BMTFile([[z]],"aprox.grid.csv");
+
+Text currentDir = GetAbsolutePath(GetFilePath(GetSourcePath(z)));
+Text WriteFile("show.R", 
+"
+library(rgl)
+
+################################################################################
+read.bbm <-
+################################################################################
+function(fbbm) {
+  fd = file(fbbm, open='rb')
+  dim <- readBin(fd,integer(), n=2)
+  data<-readBin(fd,double(),n=dim[1]*dim[2])
+  close(fd)
+  matrix(data, nrow=dim[1], ncol=dim[2], byrow=TRUE)
+}
+
+x<-t(read.bbm(\""+currentDir+"aprox.x.bbm\")) 
+ y<-t(read.bbm(\""+currentDir+"aprox.y.bbm\")) 
+ z<-read.bbm(\""+currentDir+"aprox.z.bbm\")
+xyz.real <- read.bbm(\""+currentDir+"real.xyz.bbm\")
+open3d()
+bg3d(\"white\")
+
+persp3d(x[,],y[,],z, xlim = range(x, na.rm = TRUE), ylim = range(y, na.rm = TRUE), col=\"green\")
+points3d(xyz.real)
+");
+
+Real OSCmdWait("Rgui CMD BATCH --vanilla show.R show.R.log");
+/*
 Real numErr1 = Copy(NError);
 Real numWar1 = Copy(NWarning);