See https://www.tbeardsley.com/projects/lfts/fts_gpu for a detailed discussion of this project.
Lamellar | Cylindrical | Spherical | Gyroid | Fddd |
---|---|---|---|---|
GSL - GNU Scientific Library (https://www.gnu.org/software/gsl/)
CUDA Toolkit (https://developer.nvidia.com/cuda-toolkit/)
Two methods of compiling the program are available:
- comp.sh
A simple bash script to create a 'build' directory containing the compiled program code: lfts-gpu.
On a Linux system, run the bash script from the top directory via:
sh comp.sh
- CMake
CMakeLists.txt specifies the required commands for CMake to create (and run) Makefiles, which create a 'build' directory and compile the program code as: lfts-gpu.
From the top directory, run:
cmake -B build
cmake --build build
After compilation the executable file, lfts-gpu, resides in the 'build' directory. An input file must be supplied to the executable at the command line, examples of which are contained in the 'input_files' folder.
For example, from the top level of the directory tree, the program could be run via:
./build/lfts-gpu ./input_files/input
The input_files directory contains example input files that can be supplied to the program from the command line.
Line 1: N NA XN C Ndt isXeN
Line 2: mx my mz Lx Ly Lz
Line 3: n_eq n_st n_smpl save_freq loadType
Lines 4->(M+3): W-(r)
Lines (M+4)->(2M+3): w+(r)
Note: A real-space position r = (x,y,z) corresponds to a mesh point position r_m = (i,j,k), where i=0->mx-1, j=0->my-1 and k=0->mz-1 are integers. The elements of the fields, W-(r) and w+(r), are then written in ascending order of the row-major index: p = mx*(i*my+j)+k.
Parameter | Type | Description |
---|---|---|
N | Integer | Number of monomers in a single polymer chain |
NA | Integer | Number of monomers in the A-block of a polymer chain |
XN | Double | Interaction strength between A and B-type monomers |
C | Double | Square root of the invariant polymerisation index, Nbar |
Ndt | Double | Size of the time step in the Langevin update of W-(r) |
isXeN | Integer | Whether the parameter XN is in terms of bare (isXeN=0) or effective (isXeN=1) chi |
mx, my, mz | Integers | Number of mesh points in the x, y, and z dimensions of the simulation box |
Lx, Ly, Lz | Doubles | Dimensions of the simulation box (in units of the polymer end-to-end length, R0) in the x, y, and z dimensions |
n_eq | Integer | Number of langevin steps performed to equilibrate the system |
n_st | Integer | Number of langevin steps performed after equilibration has ended, during which statistics are sampled |
n_smpl | Integer | Number of steps between samples being taken in the statistics period |
save_freq | Integer | Number of steps between saving outputs to file |
loadType | Integer | Whether to load the W-(r) and w+(r) fields from the proceeding file lines (loadType=1), start from a disordered state (loadType=0) or start from a (300) lamellar phase (loadType=2) |
M | Integer | Total number of mesh points (M = mx*my*mz), such that the proceeding 2*M lines of the file can hold the W-(r) and w+(r) fields to be loaded |
The state of the W-(r) and w+(r) fields at simulation step number <step_number> during the equilibration period. First three lines are simulation parameters so it can be used as an input file.
The state of the W-(r) and w+(r) fields at simulation step number <step_number> during the statistics gathering period. First three lines are simulation parameters so it can be used as an input file.
The state of the phi-(r) and phi+(r) fields at simulation step number <step_number> during the equilibration period.
The state of the phi-(r) and phi+(r) fields at simulation step number <step_number> during the statistics gathering period.