polymer and ellipsoid simulation in 3D

Dear all,

I have a couple of questions. Looks like I am doing something wrong as I get a segmentation fault, and perhaps I am missing something important conceptually.

My plan is to have a brownian dynamics of LJ spheres connected by a bond, and ellipsoids that are only interacting through gay-berne with the bead-spring model but are not bonded to the polymer.

Please find the code attached for your convenience.
If anyone has comments I’d really appreciate it, thank you!

With best wishes,

Anna

ellipsoid.dat (2.35 KB)

ellipsoid_test2.in (803 Bytes)