Changeset 7dd69db in tol for tol_pkg/FunRn2R/BaseC2.LinComb.tol

Timestamp:

Jan 5, 2010, 7:27:53 PM (17 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

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

Children:

c62037e

Parents:

a693055

Message:

Renombrado de métodos

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

File:

• tol_pkg/FunRn2R/BaseC2.LinComb.tol (modified) (2 diffs)

tol_pkg/FunRn2R/BaseC2.LinComb.tol

@@ -34 +34 @@
 
   ////////////////////////////////////////////////////////////////////////////
-  Static Text _.autodoc.member.Compact.Ratio = "Creates an instance of "
-  "@BaseC2.LinComb that is close to be orthonormal and has at least "
-  "the specified ratio of information of original base";
-  Static @BaseC2.LinComb Compact.Ratio(
+  Static Text _.autodoc.member.Orthogonal = "Creates an instance of "
+  "@BaseC2.LinComb that is orthonormal for given points and has exactly K "
+  "elements";
+  Static @BaseC2.LinComb Orthogonal(
     @BaseC2 base, 
-    Real informationRatio)
+    Real K, 
+    VMatrix u)
   ////////////////////////////////////////////////////////////////////////////
   {
-    Real K1 = base::_.K;
-
-    VMatrix u = Rand(K1*10, K1, -1, 1);
-
-    VMatrix Pu = Group("ConcatColumns", For(1, K1, VMatrix(Real k)
-    {
-      base::function.v(k,u)
-    }));
-
-    Set Pu.svd = SVD(VMat2Mat(Pu));
-
-    Matrix Pu.eigenvalue = Tra(SubDiag(Pu.svd[2],0));
-    Matrix Pu.eigenvalueCum = 
-    {
-       Matrix aux = DifEq(1/(1-B),Pu.eigenvalue);
-       aux*(1/MatSum(Pu.eigenvalue))
-    };
-    Matrix ratio = Rand(K1,1,informationRatio,informationRatio);
-    Real K = MatSum(LE(Pu.eigenvalueCum,ratio));
-    Set range.K = Range(1,K,1);
-    VMatrix V = Mat2VMat(SubCol(Pu.svd[3],range.K));
-    VMatrix D = Eye(K,K,0,RPow(Mat2VMat(SubRow(Pu.eigenvalue,range.K)),+1));
-    VMatrix Di = Eye(K,K,0,RPow(Mat2VMat(SubRow(Pu.eigenvalue,range.K)),-1));
-    VMatrix W = V*Di;
-    VMatrix Wi = D*Tra(V);
-    New(base, W, Wi)
-  };
-
-  ////////////////////////////////////////////////////////////////////////////
-  Static Text _.autodoc.member.Compact.K = "Creates an instance of "
-  "@BaseC2.LinComb that is close to be orthonormal and has exactly K "
-  "elements";
-  Static @BaseC2.LinComb Compact.K(
-    @BaseC2 base, 
-    Real K)
-  ////////////////////////////////////////////////////////////////////////////
-  {
-    Set aux = [[
     Set range.K = Range(1,K,1);
     Real n = base::_.n;
     Real K1 = base::_.K;
-
-    VMatrix u = Gaussian(K1, n, 0, 1);
-  //VMatrix u = Eye(K1,n);
-  //WriteLn("TRACE [@BaseC2.LinComb Compact.K] u:"<<u);
     Real base::upgrade.cache(u,0,0);
     VMatrix Pu = Group("ConcatColumns", For(1, K1, VMatrix(Real k)
@@ -100 +59 @@
     VMatrix Di = Eye(K,K,0,RPow(Mat2VMat(SubRow(Pu.eigenvalue,range.K)),-1));
     VMatrix W = V*Di;
-    VMatrix Wi = D*Tra(V) ]];
+    VMatrix Wi = D*Tra(V);
     @BaseC2.LinComb new = New(base, W, Wi);
-    Set new::_.aux := aux;
+    Set new::_.aux :=
+    [[
+      K,n,K1,u,Pu,Pu.svd,Pu.eigenvalue,V,D,Di,W,Wi 
+    ]];
     new
+  };
+
+  ////////////////////////////////////////////////////////////////////////////
+  Static Text _.autodoc.member.Orthogonal.Rand = "Creates an instance of "
+  "@BaseC2.LinComb that is orthonormal for random points and has exactly K "
+  "elements";
+  Static @BaseC2.LinComb Orthogonal.Rand(
+    @BaseC2 base, 
+    Real K)
+  ////////////////////////////////////////////////////////////////////////////
+  {
+    Orthogonal(base,K,Gaussian(10*base::_.K, base::_.n, 0, 1))
   };