Changeset 59b329d in tol

Timestamp:

Oct 17, 2012, 4:18:28 PM (14 years ago)

Author:

Víctor de Buen <vdebuen@…>

Branches:

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

Children:

c2d5bcee

Parents:

6ca24fa

Message:

Refs #1136
Fixed bug in RandStationary

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

File:

• tol/btol/polynomial_type/polgra.cpp (modified) (2 diffs)

tol/btol/polynomial_type/polgra.cpp

@@ -510 +510 @@
   BPol p1 = BPol::One();
   BPol p2 = BPol::One();
-  BUniformDist u(-1,1);
   BDat a, b, d, m1, m2;
   BPol B0 = BPol::One();
   BPol BP = BPol::X() ^ period;
   BPol factor;
+  BUniformDist u(-1,1);
+
   if(maxInvRootModule>1-2*DEpsilon()) { maxInvRootModule = 1-2*DEpsilon(); }
   if(maxInvRootModule<  2*DEpsilon()) { maxInvRootModule =   2*DEpsilon(); }
@@ -525 +526 @@
   if(n1)
   {
+    BUniformDist v(0,1);
+    BDat vr = v.Random();
     BDat ur = u.Random();
-    a = maxInvRootModule*ur + minInvRootModule*(1-ur);
+    a = maxInvRootModule*vr + minInvRootModule*(1-vr);
+    if(ur<0) { a*=-1.0; }
     factor = B0-a*BP;
   //Std(BText("\nRandStationary(")+degree+","+period+","+maxInvRootModule+") monome a="+a+"; factor="+factor.Name());