Publication details
Implementing efficient concerted rotations using Mathematica and C code
| Authors | |
|---|---|
| Year of publication | 2018 |
| Type | Article in Periodical |
| Magazine / Source | EUROPEAN PHYSICAL JOURNAL E |
| MU Faculty or unit | |
| Citation | |
| Doi | http://dx.doi.org/10.1140/epje/i2018-11694-7 |
| Keywords | INTRINSICALLY UNSTRUCTURED PROTEINS; MONTE-CARLO SIMULATIONS; COARSE-GRAINED MODELS; MOLECULAR-DYNAMICS; SYSTEMS; ALGORITHM |
| Attached files | |
| Description | In this article we demonstrate a general and efficient metaprogramming implementation of concerted rotations using Mathematica. Concerted rotations allow the movement of a fixed portion of a polymer backbone with fixed bending angles, like a protein, while maintaining the correct geometry of the backbone and the initial and final points of the portion fixed. Our implementation uses Mathematica to generate a C code which is then wrapped in a library by a Python script. The user can modify the Mathematica notebook to generate a set of concerted rotations suited for a particular backbone geometry, without having to write the C code himself. The resulting code is highly optimized, performing on the order of thousands of operations per second. |
| Related projects: |