Software Open Access

Pushing the frontiers of polymer simulation with texture memory

Airidas Korolkovas


DataCite XML Export

<?xml version='1.0' encoding='utf-8'?>
<resource xmlns:xsi="http://www.w3.org/2001/XMLSchema-instance" xmlns="http://datacite.org/schema/kernel-4" xsi:schemaLocation="http://datacite.org/schema/kernel-4 http://schema.datacite.org/meta/kernel-4.1/metadata.xsd">
  <identifier identifierType="DOI">10.5281/zenodo.1305227</identifier>
  <creators>
    <creator>
      <creatorName>Airidas Korolkovas</creatorName>
      <nameIdentifier nameIdentifierScheme="ORCID" schemeURI="http://orcid.org/">0000-0001-7904-5639</nameIdentifier>
      <affiliation>Uppsala University</affiliation>
    </creator>
  </creators>
  <titles>
    <title>Pushing the frontiers of polymer simulation with texture memory</title>
  </titles>
  <publisher>Zenodo</publisher>
  <publicationYear>2018</publicationYear>
  <subjects>
    <subject>GPU, CUDA, textures, polymers, entanglement, reptation</subject>
  </subjects>
  <dates>
    <date dateType="Issued">2018-07-04</date>
  </dates>
  <language>en</language>
  <resourceType resourceTypeGeneral="Software"/>
  <alternateIdentifiers>
    <alternateIdentifier alternateIdentifierType="url">https://zenodo.org/record/1305227</alternateIdentifier>
  </alternateIdentifiers>
  <relatedIdentifiers>
    <relatedIdentifier relatedIdentifierType="URL" relationType="IsCitedBy">https://arxiv.org/abs/1805.08508</relatedIdentifier>
    <relatedIdentifier relatedIdentifierType="DOI" relationType="IsCitedBy">10.1002/adts.201800078</relatedIdentifier>
    <relatedIdentifier relatedIdentifierType="DOI" relationType="IsVersionOf">10.5281/zenodo.1305226</relatedIdentifier>
  </relatedIdentifiers>
  <rightsList>
    <rights rightsURI="https://creativecommons.org/licenses/by/4.0/legalcode">Creative Commons Attribution 4.0 International</rights>
    <rights rightsURI="info:eu-repo/semantics/openAccess">Open Access</rights>
  </rightsList>
  <descriptions>
    <description descriptionType="Abstract">&lt;p&gt;Pushing the frontiers of polymer simulation with texture memory&lt;br&gt;
---------------------------------------------------------------&lt;/p&gt;

&lt;p&gt;&lt;br&gt;
Welcome to the Source Code for simulating entangled star polymers.&lt;br&gt;
This custom software was used to generate data for the article &amp;ldquo;Five dimensional entanglement in star polymer dynamics&amp;quot;&lt;br&gt;
published in Advanced Theory and Simulations (Wiley)&lt;br&gt;
A preprint is available here: https://arxiv.org/abs/1805.08508&lt;/p&gt;

&lt;p&gt;The demo provided here simulates Brownian dynamics of NS=104 stars with 3 arms, of length N=356 beads per arm.&lt;br&gt;
It requires a CUDA-enabled GPU card, which is managed from a Matlab environment using the launchstars.m script.&lt;br&gt;
The .cu source files are included and have been pre-compiled for several star sizes using this command:&lt;br&gt;
&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp; mexcuda mystar356.cu&lt;/p&gt;

&lt;p&gt;Now simply run:&lt;br&gt;
&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp;&amp;nbsp; launchstars.m&lt;/p&gt;

&lt;p&gt;The default settings create the initial polymer configuration in Matlab, which is then transferred to GPU and processed by CUDA for 1000 steps.&lt;br&gt;
The resulting configuration is returned to Matlab, and selected data is stored for later analysis.&lt;br&gt;
This cycle is repeated 200 times, and the real-space configuration is optionally plotted in 3D (see a snapshot in the Figures folder).&lt;br&gt;
Once the demo is complete, you may verify the speed of the computation using this command (it runs an extra numIterations=1000 steps):&lt;/p&gt;

&lt;p&gt;&amp;nbsp;&amp;nbsp; &amp;nbsp;tic; [R, randvector, ~] = fh(single(gpuArray(R)), randvector, b, vstrong, rstrong, cosbeta, numIterations, timeStep); toc;&lt;br&gt;
&amp;nbsp;&amp;nbsp; &amp;nbsp;Elapsed time is 0.398752 seconds.&lt;/p&gt;

&lt;p&gt;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.&lt;br&gt;
The whole demo lasts about 1 minute 30 seconds.&lt;br&gt;
The branch point trajectory may be analysed by running the script autodiff.m provided in the Data Analysis Tools folder.&lt;br&gt;
A typical result is shown in the MSD.pdf figure. Sub-reptative behaviour (negative slope) is apparent within one minute of simulation.&lt;/p&gt;

&lt;p&gt;&lt;br&gt;
Author:&lt;br&gt;
Airidas Korolkovas, PhD&lt;br&gt;
Institut Laue-Langevin&lt;br&gt;
Uppsala University&lt;br&gt;
korolkovas@ill.fr&lt;br&gt;
airidas.korolkovas89@gmail.com&lt;/p&gt;

&lt;p&gt;July 04, 2018&lt;/p&gt;</description>
  </descriptions>
</resource>
1,576
50
views
downloads
All versions This version
Views 1,5761,578
Downloads 5050
Data volume 1.3 GB1.3 GB
Unique views 1,4761,478
Unique downloads 4444

Share

Cite as