Changeset 9 for tests

Show
Ignore:
Timestamp:
03/14/07 10:27:52 (2 years ago)
Author:
pernet
Message:

v1.1.2
Change the detection of field representation using double (previous one was not working)

Files:
1 modified

Legend:

Unmodified
Added
Removed
  • tests/test-invert.C

    r1 r9  
    2727        int n; 
    2828        int nbit=atoi(argv[3]); // number of times the product is performed 
    29         long alphalong=1, betalong=0; 
    3029        cerr<<setprecision(10); 
    31         Field::Element alpha,beta; 
    3230        Field::Element zero, one; 
    3331 
     
    4543        A = read_field(F,argv[2],&n,&n); 
    4644        Ab = new Field::Element[n*n]; 
    47         for (size_t i=0; i<n*n;++i) 
     45        for (int i=0; i<n*n;++i) 
    4846                F.assign(*(Ab+i),*(A+i)); 
    4947        Field::Element * X = new Field::Element[n*n]; 
    5048 
    5149        Timer tim,t; t.clear();tim.clear();  
    52         int nullity; 
     50        int nullity=0; 
     51         
    5352        for(int i = 0;i<nbit;++i){ 
    5453                t.clear(); 
     
    5857                tim+=t; 
    5958                if (i+1<nbit) 
    60                         for (size_t i=0; i<n*n;++i) 
     59                        for (int i=0; i<n*n;++i) 
    6160                                F.assign(*(A+i),*(Ab+i)); 
    6261        } 
     
    6867        bool wrong = false; 
    6968 
    70         for (size_t i=0;i<n;++i) 
    71                 for (size_t j=0;j<n;++j) 
     69        for (int i=0;i<n;++i) 
     70                for (int j=0;j<n;++j) 
    7271                        if ( ((i!=j) && !F.areEqual(*(I+i*n+j),zero)) 
    7372                             ||((i==j) &&!F.areEqual(*(I+i*n+j),one)))