X-Git-Url: https://hackdaworld.org/gitweb/?a=blobdiff_plain;f=moldyn.c;h=c27fc6215606df4d223290ea2e08e1798d64e95c;hb=4b18a09781c27d786c32a5ada98929c4c2e38f4b;hp=d7d89bc2328506715a1d5d82e439267ce41281fa;hpb=b397d0e80a7191ae1f10a86e9a595dcb8668286e;p=physik%2Fposic.git

diff --git a/moldyn.c b/moldyn.c
index d7d89bc..c27fc62 100644
--- a/moldyn.c
+++ b/moldyn.c
@@ -47,6 +47,12 @@
 #undef PSE_NAME
 #undef PSE_COL
 
+#ifdef PTHREADS
+/* global mutexes */
+pthread_mutex_t *amutex;
+pthread_mutex_t emutex;
+#endif
+
 /*
  * the moldyn functions
  */
@@ -66,13 +72,27 @@ int moldyn_init(t_moldyn *moldyn,int argc,char **argv) {
 	rand_init(&(moldyn->random),NULL,1);
 	moldyn->random.status|=RAND_STAT_VERBOSE;
 
+#ifdef PTHREADS
+	pthread_mutex_init(&emutex,NULL);
+#endif
+
 	return 0;
 }
 
 int moldyn_shutdown(t_moldyn *moldyn) {
 
+#ifdef PTHREADS
+	int i;
+#endif
+
 	printf("[moldyn] shutdown\n");
 
+#ifdef PTHREADS
+	for(i=0;i<moldyn->count;i++)
+		pthread_mutex_destroy(&(amutex[i]));
+	pthread_mutex_destroy(&emutex);
+#endif
+
 	moldyn_log_shutdown(moldyn);
 	link_cell_shutdown(moldyn);
 	rand_close(&(moldyn->random));
@@ -500,6 +520,9 @@ int create_lattice(t_moldyn *moldyn,u8 type,double lc,int element,double mass,
 	void *ptr;
 	t_atom *atom;
 	char name[16];
+#ifdef PTHREADS
+	pthread_mutex_t *mutex;
+#endif
 
 	new=a*b*c;
 	count=moldyn->count;
@@ -522,6 +545,16 @@ int create_lattice(t_moldyn *moldyn,u8 type,double lc,int element,double mass,
 	moldyn->atom=ptr;
 	atom=&(moldyn->atom[count]);
 
+#ifdef PTHREADS
+	ptr=realloc(amutex,(count+new)*sizeof(pthread_mutex_t));
+	if(!ptr) {
+		perror("[moldyn] mutex realloc (add atom)");
+		return -1;
+	}
+	amutex=ptr;
+	mutex=&(amutex[count]);
+#endif
+
 	/* no atoms on the boundaries (only reason: it looks better!) */
 	if(!origin) {
 		orig.x=0.5*lc;
@@ -579,6 +612,9 @@ int create_lattice(t_moldyn *moldyn,u8 type,double lc,int element,double mass,
 		atom[ret].tag=count+ret;
 		check_per_bound(moldyn,&(atom[ret].r));
 		atom[ret].r_0=atom[ret].r;
+#ifdef PTHREADS
+		pthread_mutex_init(&(mutex[ret]),NULL);
+#endif
 	}
 
 	/* update total system mass */
@@ -613,6 +649,16 @@ int add_atom(t_moldyn *moldyn,int element,double mass,u8 brand,u8 attr,
 	moldyn->lc.subcell->list=ptr;
 #endif
 
+#ifdef PTHREADS
+	ptr=realloc(amutex,(count+1)*sizeof(pthread_mutex_t));
+	if(!ptr) {
+		perror("[moldyn] mutex realloc (add atom)");
+		return -1;
+	}
+	amutex=ptr;
+	pthread_mutex_init(&(amutex[count]),NULL);
+#endif
+
 	atom=moldyn->atom;
 
 	/* initialize new atom */