void computeForces(SystemInterface &s) {
dds_float box[3];
int per[3];
for (int i=0;i<3;i++)
{
box[i]=s.box()[i];
per[i] = (PERIODIC(i));
}
DipolarDirectSum_kernel_wrapper_force(k,s.npart_gpu(),
s.rGpuBegin(), s.dipGpuBegin(), s.fGpuBegin(),s.torqueGpuBegin(),box,per);
};
virtual void computeForces(SystemInterface &s) {
printf("Trace HarmonicWell computeForces 1");
HarmonicWell_kernel_wrapper(x,y,z,k,s.npart_gpu(),
s.rGpuBegin(), s.fGpuBegin());
};