testcollapse.cpp

Go to the documentation of this file.

#include <iostream>
#include "CollapseModel.h"
#include "UniformRNG.h"
#include <boost/numeric/ublas/io.hpp>

using namespace std;
int main()
{
        const unsigned int size = 4;
        const unsigned int tests = 10;
        
        UniformRNG Random;
        cout << "One identical layer"; 
        for (int t = 0; t < tests; ++t)
        {
                ttranscribed thick(size), value(size);
                
                for (int i = 0; i < size; ++i)
                {
                        thick(i) = Random.GetNumber();
                        value(i) = Random.GetNumber();
                }
                value(1) = value(2);
                cout << "Before: " << endl;
                cout << "  Thick: " << thick << endl;
                cout << "  Value: " << value << endl;
                CollapseModel(thick,value);
                cout << "After: " << endl;
                 cout << "  Thick: " << thick << endl;
                cout << "  Value: " << value << endl << endl;
        }
        cout << "One nearly identical layer"; 
        double eps = 1e-10;
        
        
                
        for (int t = 0; t < tests; ++t)
        {
                ttranscribed thick(size), value(size);
                for (int i = 0; i < size; ++i)
                {
                        thick(i) = Random.GetNumber();
                        value(i) = Random.GetNumber();
                }
                
                
                value(1) = value(2)+eps*t;
                cout << "Eps: " << eps*t << endl;
                cout << "Before: " << endl;
                cout << "  Thick: " << thick << endl;
                cout << "  Value: " << value << endl;
                CollapseModel(thick,value);
                cout << "After: " << endl;
                 cout << "  Thick: " << thick << endl;
                cout << "  Value: " << value << endl << endl;
        }
}

---