Changeset 1112 in Sophya


Ignore:
Timestamp:
Jul 28, 2000, 6:34:00 PM (25 years ago)
Author:
ansari
Message:

adaptation aux nouvelles Sequences et pb de int/uint cmv 28/7/00

Location:
trunk/SophyaProg/Tests
Files:
2 edited

Legend:

Unmodified
Added
Removed
  • trunk/SophyaProg/Tests/tsttmat.cc

    r857 r1112  
    112112TMatrix<r_8> a(A,false);
    113113cout<<endl<<"Matrices AAA*=C"<<endl;
    114 {for(int i=0;i<AAA.NRows();i++) for(int j=0;j<C.NCols();j++)
    115   {a(i,j) = 0; for(int k=0;k<AAA.NCols();k++) a(i,j) += AAA(i,k)*C(k,j);}}
     114{for(uint_4 i=0;i<AAA.NRows();i++) for(uint_4 j=0;j<C.NCols();j++)
     115  {a(i,j) = 0; for(uint_4 k=0;k<AAA.NCols();k++) a(i,j) += AAA(i,k)*C(k,j);}}
    116116AAA *= C;
    117117cout<<"AAA:"<<endl<<AAA;
     
    119119TMatrix<r_8> b(B,false);
    120120cout<<"Matrices BBB*=D"<<endl;
    121 {for(int i=0;i<BBB.NRows();i++) for(int j=0;j<D.NCols();j++)
    122   {b(i,j) = 0; for(int k=0;k<BBB.NCols();k++) b(i,j) += BBB(i,k)*D(k,j);}}
     121{for(uint_4 i=0;i<BBB.NRows();i++) for(uint_4 j=0;j<D.NCols();j++)
     122  {b(i,j) = 0; for(uint_4 k=0;k<BBB.NCols();k++) b(i,j) += BBB(i,k)*D(k,j);}}
    123123BBB *= D;
    124124cout<<"B:"<<endl<<BBB;
     
    158158//-----------------------------------------------
    159159TMatrix<r_8> RR(2,4);
    160 {for(int i=0;i<A.NRows();i++) for(int j=0;j<E.NCols();j++)
    161   {RR(i,j) = 0; for(int k=0;k<A.NCols();k++) RR(i,j) += A(i,k)*E(k,j);}}
     160{for(uint_4 i=0;i<A.NRows();i++) for(uint_4 j=0;j<E.NCols();j++)
     161  {RR(i,j) = 0; for(uint_4 k=0;k<A.NCols();k++) RR(i,j) += A(i,k)*E(k,j);}}
    162162cout<<endl<<"Matrices R = A * E"<<endl;
    163163R = A * E;
  • trunk/SophyaProg/Tests/tsttminv.cc

    r1008 r1112  
    138138Vind = 0;
    139139if(nalea>0) for(i=0;i<nalea;i++) drand01();
    140 A = Sequence(RandomSequence(RandomSequence::Flat,0.,1.));
     140A = RandomSequence(RandomSequence::Flat,0.,1.);
    141141if(nbig>0) for(i=0;i<nbig;i++) {
    142142  k=(uint_8)(drand01()*N*N); if(k>=N*N) k--;
     
    147147#if defined(COMPLEX)
    148148Vind = 0;
    149 B = Sequence(RandomSequence(RandomSequence::Flat,0.,1.));
     149B = RandomSequence(RandomSequence::Flat,0.,1.);
    150150if(nbig>0) for(i=0;i<nbig;i++) {
    151151  k=(uint_8)(drand01()*N*N); if(k>=N*N) k--;
Note: See TracChangeset for help on using the changeset viewer.