+#ifdef PTHREADS
+ kdata[k].moldyn=moldyn;
+ kdata[k].ai=&(itom[i]);
+ kdata[k].aj=jtom;
+ kdata[k].ak=ktom;
+ kdata[k].bc=bc_ik;
+ ret=pthread_create(&(kthread[k]),NULL,moldyn->func3b_k1,&(kdata[k]));
+ if(ret) {
+ perror("[moldyn] create k1 thread");
+ return ret;
+ }
+#else