Ignore:
Timestamp:
Aug 18, 2009, 8:48:06 AM (16 years ago)
Author:
Frederik Heber <heber@…>
Children:
66ce7a
Parents:
3de1d2
git-author:
Frederik Heber <heber@…> (08/18/09 08:38:46)
git-committer:
Frederik Heber <heber@…> (08/18/09 08:48:06)
Message:

Incorporation of Unit test on class Vector.

  • new file leastsquaremin.[ch]pp has least square minimisation which is otherwise unclean between classes molecules and Vector

Unit test (later tests rely on good results of earlier ones)

changes to class Vector:

  • Vector::IsNull() -> IsZero()
  • new function Vector::IsOne() similar to IsZero()
  • BUGFIX: Vector::IsNULL() (thx to unit test :)
  • Tesselation::getAngle() changed due to above rename
File:
1 edited

Legend:

Unmodified
Added
Removed
  • molecuilder/src/molecules.hpp

    r3de1d2 rda84d6  
    2828#include "bond.hpp"
    2929#include "element.hpp"
     30#include "leastsquaremin.hpp"
    3031#include "linkedcell.hpp"
    3132#include "parser.hpp"
     
    8384
    8485
    85 // some algebraic matrix stuff
    86 #define RDET3(a) ((a)[0]*(a)[4]*(a)[8] + (a)[3]*(a)[7]*(a)[2] + (a)[6]*(a)[1]*(a)[5] - (a)[2]*(a)[4]*(a)[6] - (a)[5]*(a)[7]*(a)[0] - (a)[8]*(a)[1]*(a)[3])  //!< hard-coded determinant of a 3x3 matrix
    87 #define RDET2(a0,a1,a2,a3) ((a0)*(a3)-(a1)*(a2))                      //!< hard-coded determinant of a 2x2 matrix
    88 
    89 
    90 /** Parameter structure for least square minimsation.
    91  */
    92 struct LSQ_params {
    93   Vector **vectors;
    94   int num;
    95 };
    96 
    97 double LSQ(const gsl_vector * x, void * params);
    98 
    99 /** Parameter structure for least square minimsation.
    100  */
    101 struct lsq_params {
    102   gsl_vector *x;
    103   const molecule *mol;
    104   element *type;
    105 };
    10686
    10787#define MaxThermostats 6      //!< maximum number of thermostat entries in Ions#ThermostatNames and Ions#ThermostatImplemented
Note: See TracChangeset for help on using the changeset viewer.