Changeset 6ccef2a in tol


Ignore:
Timestamp:
Apr 2, 2009, 5:28:28 PM (18 years ago)
Author:
Víctor de Buen <vdebuen@…>
Branches:
1941-docker-image-to-test, devel, master
Children:
731ece8
Parents:
b83ba94d
Message:

Fixed error building bounds

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

File:
1 edited

Legend:

Unmodified
Added
Removed
  • tol_tests/tol/random_generation/test_0003/test.tol

    rb83ba94d r6ccef2a  
    4444  Real   k1  = 0;
    4545  Real   k2  = +5*sigma*sk;
    46   Matrix a00 = Col(k0,k0); b00 = L*Col(k1,k1);
    47   Matrix a01 = Col(k1,k0); b01 = L*Col(k2,k1);
    48   Matrix a10 = Col(k0,k1); b10 = L*Col(k1,k2);
    49   Matrix a11 = Col(k1,k1); b11 = L*Col(k2,k2);
     46  Matrix a00 = Col(k0,k0); b00 = Col(k1,k1);
     47  Matrix a01 = Col(k1,k0); b01 = Col(k2,k1);
     48  Matrix a10 = Col(k0,k1); b10 = Col(k1,k2);
     49  Matrix a11 = Col(k1,k1); b11 = Col(k2,k2);
    5050  Matrix bounds = Tra(a00|b00|a01|b01|a10|b10|a11|b11);
    5151//Sampling non truncated bidimensional normal to get area for each
    … …  
    7676
    7777//Testing coherence of restrictions
    78   Real r00 = SetProd(For(1,Columns(c00),Real(Real j)
     78  Real r00 = If(!Columns(c00),0,SetProd(For(1,Columns(c00),Real(Real j)
    7979  {
    8080    Matrix cj = SubCol(c00,[[j]]);
    8181    Real And( MatSum(LE(a00,cj))==n, MatSum(LE(cj,b00))==n)
    82   }));
    83   Real r01 = SetProd(For(1,Columns(c01),Real(Real j)
     82  })));
     83  Real r01 = If(!Columns(c01),0,SetProd(For(1,Columns(c01),Real(Real j)
    8484  {
    8585    Matrix cj = SubCol(c01,[[j]]);
    8686    Real And( MatSum(LE(a01,cj))==n, MatSum(LE(cj,b01))==n)
    87   }));
    88   Real r10 = SetProd(For(1,Columns(c10),Real(Real j)
     87  })));
     88  Real r10 = If(!Columns(c10),0,SetProd(For(1,Columns(c10),Real(Real j)
    8989  {
    9090    Matrix cj = SubCol(c10,[[j]]);
    9191    Real And( MatSum(LE(a10,cj))==n, MatSum(LE(cj,b10))==n)
    92   }));
    93   Real r11 = SetProd(For(1,Columns(c11),Real(Real j)
     92  })));
     93  Real r11 = If(!Columns(c11),0,SetProd(For(1,Columns(c11),Real(Real j)
    9494  {
    9595    Matrix cj = SubCol(c11,[[j]]);
    9696    Real And( MatSum(LE(a11,cj))==n, MatSum(LE(cj,b11))==n)
    97   }));
     97  })));
    9898//Joining 4 samples
    9999  Matrix c   = c00|c01|c10|c11;
Note: See TracChangeset for help on using the changeset viewer.