Changeset 6ccef2a in tol
- Timestamp:
- Apr 2, 2009, 5:28:28 PM (18 years ago)
- Branches:
- 1941-docker-image-to-test, devel, master
- Children:
- 731ece8
- Parents:
- b83ba94d
- File:
-
- 1 edited
Legend:
- Unmodified
- Added
- Removed
-
tol_tests/tol/random_generation/test_0003/test.tol
rb83ba94d r6ccef2a 44 44 Real k1 = 0; 45 45 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); 50 50 Matrix bounds = Tra(a00|b00|a01|b01|a10|b10|a11|b11); 51 51 //Sampling non truncated bidimensional normal to get area for each … … 76 76 77 77 //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) 79 79 { 80 80 Matrix cj = SubCol(c00,[[j]]); 81 81 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) 84 84 { 85 85 Matrix cj = SubCol(c01,[[j]]); 86 86 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) 89 89 { 90 90 Matrix cj = SubCol(c10,[[j]]); 91 91 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) 94 94 { 95 95 Matrix cj = SubCol(c11,[[j]]); 96 96 Real And( MatSum(LE(a11,cj))==n, MatSum(LE(cj,b11))==n) 97 })) ;97 }))); 98 98 //Joining 4 samples 99 99 Matrix c = c00|c01|c10|c11;
Note: See TracChangeset
for help on using the changeset viewer.
