Skip to content

Commit 5a751a7

Browse files
committed
pass 2e-5 converge test
1 parent 24c7b97 commit 5a751a7

File tree

5 files changed

+160
-1498
lines changed

5 files changed

+160
-1498
lines changed

src/minimize/minimize.cu

Lines changed: 24 additions & 51 deletions
Original file line numberDiff line numberDiff line change
@@ -20,7 +20,6 @@ The driver class for minimizers.
2020
#include "force/force.cuh"
2121
#include "minimize.cuh"
2222
#include "minimizer_fire.cuh"
23-
#include "minimizer_fire2.cuh"
2423
#include "minimizer_fire_box_change.cuh"
2524
#include "minimizer_sd.cuh"
2625
#include "utilities/error.cuh"
@@ -139,30 +138,18 @@ void Minimize::parse_minimize(
139138
printf(" with fixed box.\n");
140139
printf(" with a force tolerance of %g eV/A.\n", force_tolerance);
141140
printf(" for maximally %d steps.\n", number_of_steps);
142-
if (number_of_atoms > 999) {
143-
minimizer.reset(new Minimizer_FIRE(number_of_atoms, number_of_steps, force_tolerance));
144-
145-
minimizer->compute(
146-
force,
147-
box,
148-
position_per_atom,
149-
type,
150-
group,
151-
potential_per_atom,
152-
force_per_atom,
153-
virial_per_atom);
154-
} else {
155-
minimizer.reset(new Minimizer_FIRE2(number_of_atoms, number_of_steps, force_tolerance));
156-
minimizer->compute(
157-
force,
158-
box,
159-
position_per_atom,
160-
type,
161-
group,
162-
potential_per_atom,
163-
force_per_atom,
164-
virial_per_atom);
165-
}
141+
142+
minimizer.reset(new Minimizer_FIRE(number_of_atoms, number_of_steps, force_tolerance));
143+
144+
minimizer->compute(
145+
force,
146+
box,
147+
position_per_atom,
148+
type,
149+
group,
150+
potential_per_atom,
151+
force_per_atom,
152+
virial_per_atom);
166153

167154
break;
168155
case 2:
@@ -175,32 +162,18 @@ void Minimize::parse_minimize(
175162
printf(" with a force tolerance of %g eV/A.\n", force_tolerance);
176163
printf(" for maximally %d steps.\n", number_of_steps);
177164

178-
if (number_of_atoms > 999) {
179-
minimizer.reset(new Minimizer_FIRE_Box_Change(
180-
number_of_atoms, number_of_steps, force_tolerance, hydrostatic_strain));
181-
182-
minimizer->compute(
183-
force,
184-
box,
185-
position_per_atom,
186-
type,
187-
group,
188-
potential_per_atom,
189-
force_per_atom,
190-
virial_per_atom);
191-
} else {
192-
minimizer.reset(new Minimizer_FIRE2(
193-
number_of_atoms, number_of_steps, force_tolerance, true, hydrostatic_strain));
194-
minimizer->compute(
195-
force,
196-
box,
197-
position_per_atom,
198-
type,
199-
group,
200-
potential_per_atom,
201-
force_per_atom,
202-
virial_per_atom);
203-
}
165+
minimizer.reset(new Minimizer_FIRE_Box_Change(
166+
number_of_atoms, number_of_steps, force_tolerance, hydrostatic_strain));
167+
168+
minimizer->compute(
169+
force,
170+
box,
171+
position_per_atom,
172+
type,
173+
group,
174+
potential_per_atom,
175+
force_per_atom,
176+
virial_per_atom);
204177

205178
break;
206179
default:

0 commit comments

Comments
 (0)