#include "Test.h"

#include "Affine_combination.h"
#include "Atom_impl.h"


namespace MM
{

class Linear_combination_test : public Test
{
    Affine_combination *link;

public:
    Linear_combination_test (const Text& suite_name) : Test (suite_name) { }

    void run()
    {
        test();   count_tests();
    }

    void test()
    {
        /*Atom_impl a_1 (Element ("C")); 
        Atom_impl a_2 (Element ("C"));
        a_1.set_velocity (Vector_3D_impl (3.,4.,-1.));
        a_2.set_velocity (Vector_3D_impl (-1.,0.,2.));

        Arbitrary_atom_group model_1;
        Arbitrary_atom_group model_2;
        model_1.add (a_1);
        model_2.add (a_2);

        Velocities  vector_1 (model_1);
        Velocities  vector_2 (model_2);

        vector_1.linear_combination (2). plus (6., vector_2);
        
        vector_1 -= Vector_3D_impl (0., 8., 10.);
        TEST ("Result is vector (0, 8, 10).", vector_1.length () < 0.000001);//*/
    }
};

Linear_combination_test test_linear_combination ("correctness");

}//MM
