Changeset 747dbd9 in tol for tol_tests

Timestamp:

May 6, 2014, 10:21:34 AM (12 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

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

Children:

2741c88

Parents:

a915a03

Message:

Refs #1770

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

Location:

tol_tests/tol/LinearAlgebra/PolMatrix

Files:

• test_0001/check.tol (modified) (1 diff)
• test_0002/check.tol (modified) (1 diff)

tol_tests/tol/LinearAlgebra/PolMatrix/test_0001/check.tol

@@ -42 +42 @@
 }));
 
+Matrix x = GetNumeric([[X]]);
+Matrix xP = GetNumeric([[XP]]);
 
+Real quality.dif.xP = 1-MatMax(Abs(x-xP));
+
+Set X2 = DeepCopy([[X]],x*2);
+Matrix x2 = GetNumeric([[X2]]);
+Real quality.dif.x2 = 1-MatMax(Abs(x2-x*2));
+
+

tol_tests/tol/LinearAlgebra/PolMatrix/test_0002/check.tol

@@ -1 +1 @@
 #Require VariableChange;
+#Require AlgLib;
+#Require MPM;
 
 Real dim = 2;
-Real deg = 3;
+Real deg = 2;
+Real T = 100;
 
-Matrix x = VariableChange::StableMatrixPolyn::
+Matrix x = VariableChange::@StableMatrixPolyn::
   SquareSetGen.Rand(dim, deg, -1, 1, 0.001);
 
-Set PHI = VariableChange::StableMatrixPolyn::
+Set PHI = VariableChange::@StableMatrixPolyn::
   DuffinSchurBinFact(dim, deg, (x));
 Polyn determinant = VecSer::PolynMatrixDeterminant(PHI);
-
 
 PolMatrix P = TabPol2PolMat(PHI);
 
 Polyn determinant_ = Determinant(P);
+Set PHI.det.roots = MPM::@ModuleVAR::PolynRootsReport(determinant_);
+Real root.abs.min = PHI.det.roots[1][3];
+Real P.isStationary = PolMatIsStationary(P);
+
+Real inequality_constraint_ = -Log(StationaryValue(determinant));
+
+Matrix X = Rand(T,dim,-1,1);
+
+Matrix Y = P:X;
+
+
+VMatrix C = CompanionMatrix(P);
+Matrix eigen = AlgLib::CppTools::rmatrixevd.eigen(VMat2Mat(C));
+Matrix eigen.module = Sqrt(eigen^2 * Col(1,1));
+Real eigen.module.max = MatMax(eigen.module);
+Real eigen.module.max.inv = 1/eigen.module.max;
+
+
+/*
+Real m = VRows(C);
+Real S=10000*m;
+VMatrix U = Gaussian(m,S,0,1);
+VMatrix V = C*U;
+VMatrix u2 = (U^2)*Constant(S,1,1/S);
+VMatrix v2 = (V^2)*Constant(S,1,1/S);
+VMatrix r2 = v2 $/ u2;
+//VMatrix R = MtMSqr(V) $/ MtMSqr(U);
+//Real VMatMax(R)
+/* */