Software Open Access
{ "inLanguage": { "alternateName": "eng", "@type": "Language", "name": "English" }, "description": "<p>Pushing the frontiers of polymer simulation with texture memory<br>\n---------------------------------------------------------------</p>\n\n<p><br>\nWelcome to the Source Code for simulating entangled star polymers.<br>\nThis custom software was used to generate data for the article “Five dimensional entanglement in star polymer dynamics"<br>\npublished in Advanced Theory and Simulations (Wiley)<br>\nA preprint is available here: https://arxiv.org/abs/1805.08508</p>\n\n<p>The demo provided here simulates Brownian dynamics of NS=104 stars with 3 arms, of length N=356 beads per arm.<br>\nIt requires a CUDA-enabled GPU card, which is managed from a Matlab environment using the launchstars.m script.<br>\nThe .cu source files are included and have been pre-compiled for several star sizes using this command:<br>\n mexcuda mystar356.cu</p>\n\n<p>Now simply run:<br>\n launchstars.m</p>\n\n<p>The default settings create the initial polymer configuration in Matlab, which is then transferred to GPU and processed by CUDA for 1000 steps.<br>\nThe resulting configuration is returned to Matlab, and selected data is stored for later analysis.<br>\nThis cycle is repeated 200 times, and the real-space configuration is optionally plotted in 3D (see a snapshot in the Figures folder).<br>\nOnce the demo is complete, you may verify the speed of the computation using this command (it runs an extra numIterations=1000 steps):</p>\n\n<p> tic; [R, randvector, ~] = fh(single(gpuArray(R)), randvector, b, vstrong, rstrong, cosbeta, numIterations, timeStep); toc;<br>\n Elapsed time is 0.398752 seconds.</p>\n\n<p>The above figure is for a box of 104*3*356=111072 beads, 1000 iterations, using a Titan XP card donated by the NVIDIA Corporation.<br>\nThe whole demo lasts about 1 minute 30 seconds.<br>\nThe branch point trajectory may be analysed by running the script autodiff.m provided in the Data Analysis Tools folder.<br>\nA typical result is shown in the MSD.pdf figure. Sub-reptative behaviour (negative slope) is apparent within one minute of simulation.</p>\n\n<p><br>\nAuthor:<br>\nAiridas Korolkovas, PhD<br>\nInstitut Laue-Langevin<br>\nUppsala University<br>\nkorolkovas@ill.fr<br>\nairidas.korolkovas89@gmail.com</p>\n\n<p>July 04, 2018</p>", "license": "https://creativecommons.org/licenses/by/4.0/legalcode", "creator": [ { "affiliation": "Uppsala University", "@id": "https://orcid.org/0000-0001-7904-5639", "@type": "Person", "name": "Airidas Korolkovas" } ], "url": "https://zenodo.org/record/1305227", "datePublished": "2018-07-04", "keywords": [ "GPU, CUDA, textures, polymers, entanglement, reptation" ], "@context": "https://schema.org/", "identifier": "https://doi.org/10.5281/zenodo.1305227", "@id": "https://doi.org/10.5281/zenodo.1305227", "@type": "SoftwareSourceCode", "name": "Pushing the frontiers of polymer simulation with texture memory" }
All versions | This version | |
---|---|---|
Views | 1,570 | 1,572 |
Downloads | 50 | 50 |
Data volume | 1.3 GB | 1.3 GB |
Unique views | 1,470 | 1,472 |
Unique downloads | 44 | 44 |