X-Git-Url: https://www.hackdaworld.org/gitweb/?a=blobdiff_plain;f=moldyn.h;h=30c4bf89527b6d023534ca02f2df1845711f70ba;hb=eceebe3ee412aa8cea3e6a7f0038883707f78460;hp=65f0d8da10cfa0753b72ffdeac86dadf7d407917;hpb=3ffe2a08e25fc091b6241885055450009267e2d8;p=physik%2Fposic.git diff --git a/moldyn.h b/moldyn.h index 65f0d8d..30c4bf8 100644 --- a/moldyn.h +++ b/moldyn.h @@ -22,9 +22,33 @@ typedef struct s_atom { } t_atom; -/* defines */ +typedef struct s_moldyn { + int count; + t_atom *atom; + double (*potential)(void *ptr); + int (*force)(struct s_moldyn *moldyn,void *ptr); + unsigned char status; +} t_moldyn; -#define K_BOLTZMANN 1.3807E-23 +typedef struct s_lj_params { + double sigma6; + double sigma12; + double epsilon; +} t_lj_params; + +/* + * defines + */ + +/* general defines */ + +#define MOLDYN_STAT_POTENTIAL 0x01 +#define MOLDYN_STAT_FORCE 0x02 + +/* phsical values */ + +//#define K_BOLTZMANN 1.3807E-23 +#define K_BOLTZMANN 1.0 #define FCC 0x01 #define DIAMOND 0x02 @@ -44,5 +68,8 @@ int destroy_lattice(t_atom *atom); int thermal_init(t_atom *atom,t_random *random,int count,double t); int scale_velocity(t_atom *atom,int count,double t); double get_e_kin(t_atom *atom,int count); +t_3dvec get_total_p(t_atom *atom,int count); + +double potential_lennard_jones_2(t_moldyn *moldyn,void *ptr); #endif