|
24 | 24 |
|
25 | 25 | jl_parameters = LennardJonesParameters(ϵOO, σOO, R) |
26 | 26 | e_parameters = ElectrostaticParameters(k_el, Rel) |
27 | | - spc_paramters = SPCFwParameters(rOH, ∠HOH, k_bond, k_angle) |
| 27 | + spc_parameters = SPCFwParameters(rOH, ∠HOH, k_bond, k_angle) |
28 | 28 | pbc = CubicPeriodicBoundaryConditions(L) |
29 | 29 |
|
30 | 30 | @testset "Analyzing simulation result" begin |
|
43 | 43 |
|
44 | 44 | bodies = [p1, p2, p3] |
45 | 45 | water = WaterSPCFw(bodies, mH, mO, qH, qO, jl_parameters, e_parameters, |
46 | | - spc_paramters) |
| 46 | + spc_parameters) |
47 | 47 | simulation = NBodySimulation(water, (t1, t2), pbc, kb) |
48 | 48 |
|
49 | 49 | result = run_simulation(simulation, VelocityVerlet(), dt = τ) |
|
96 | 96 |
|
97 | 97 | bodies = [p1, p2] |
98 | 98 | water = WaterSPCFw(bodies, mH, mO, qH, qO, jl_parameters, e_parameters, |
99 | | - spc_paramters) |
| 99 | + spc_parameters) |
100 | 100 |
|
101 | 101 | thermostat = BerendsenThermostat(T0, 200τ) |
102 | 102 | simulation = NBodySimulation(water, (t1, t2), pbc, kb) |
|
123 | 123 | T0 = 275 |
124 | 124 | bodies = generate_bodies_in_cell_nodes(N, mH2O, v_dev, L) |
125 | 125 | water = WaterSPCFw(bodies, mH, mO, qH, qO, jl_parameters, e_parameters, |
126 | | - spc_paramters) |
| 126 | + spc_parameters) |
127 | 127 | thermostat = LangevinThermostat(T0, 100) |
128 | 128 | simulation = NBodySimulation(water, (t1, t2), pbc, thermostat, kb) |
129 | 129 | result = run_simulation(simulation, EM(), dt = τ) |
|
146 | 146 |
|
147 | 147 | bodies = [p1, p2] |
148 | 148 | water = WaterSPCFw(bodies, mH, mO, qH, qO, jl_parameters, e_parameters, |
149 | | - spc_paramters) |
| 149 | + spc_parameters) |
150 | 150 | simulation = NBodySimulation(water, (t1, t2), pbc, kb) |
151 | 151 | result = run_simulation(simulation, VelocityVerlet(), dt = τ) |
152 | 152 |
|
153 | 153 | io = IOBuffer() |
154 | 154 | pdb_data = sprint(io -> NBodySimulator.write_pdb_data(io, result)) |
155 | | - splitted_data = split(pdb_data, '\n') |
| 155 | + split_data = split(pdb_data, '\n') |
156 | 156 |
|
157 | 157 | hetatm_count = 0 |
158 | 158 | timestep_count = 0 |
159 | 159 |
|
160 | | - for s in splitted_data |
| 160 | + for s in split_data |
161 | 161 | if length(s) >= 10 && s[1:10] == "REMARK 250" |
162 | 162 | timestep_count += 1 |
163 | 163 | elseif length(s) >= 6 && s[1:6] == "HETATM" |
@@ -199,7 +199,7 @@ ENDMDL") |
199 | 199 |
|
200 | 200 | bodies = NBodySimulator.extract_from_pdb(pdb_data) |
201 | 201 | water = WaterSPCFw(bodies, mH, mO, qH, qO, jl_parameters, e_parameters, |
202 | | - spc_paramters) |
| 202 | + spc_parameters) |
203 | 203 | t1 = 0τ |
204 | 204 | t2 = 2τ |
205 | 205 | simulation = NBodySimulation(water, (t1, t2), pbc, kb) |
|
0 commit comments