@@ -201,13 +201,13 @@ def write_control_block(parameters_control: dict[str, Any]) -> str:
201201 _time = default_timestep [parameters_control .get ("units" , "si" )]
202202 control_block = generate_header ("Start of the Control information" )
203203 control_block += "clear\n "
204- control_block += f' units { parameters_control .get (" units" , "si" )} \n '
205- control_block += f' newton { parameters_control .get (" newton" , "on" )} \n '
204+ control_block += f" units { parameters_control .get (' units' , 'si' )} \n "
205+ control_block += f" newton { parameters_control .get (' newton' , 'on' )} \n "
206206 if "processors" in parameters_control :
207207 control_block += (
208- f' processors { join_keywords (parameters_control [" processors" ])} \n '
208+ f" processors { join_keywords (parameters_control [' processors' ])} \n "
209209 )
210- control_block += f' timestep { parameters_control .get (" timestep" , _time )} \n '
210+ control_block += f" timestep { parameters_control .get (' timestep' , _time )} \n "
211211 control_block += generate_header ("End of the Control information" )
212212 return control_block
213213
@@ -245,27 +245,27 @@ def write_potential_block(
245245 potential_block = generate_header ("Start of Potential information" )
246246 potential_block += f"pair_style { potential .pair_style } "
247247 potential_block += (
248- f' { " " .join (parameters_potential .get (" potential_style_options" , ["" ]))} \n '
248+ f" { ' ' .join (parameters_potential .get (' potential_style_options' , ['' ]))} \n "
249249 )
250250
251251 if default_potential [potential .pair_style ].get ("read_from_file" ):
252- potential_block += f' pair_coeff * * { potential_file } { " " .join (kind_symbols )} \n '
252+ potential_block += f" pair_coeff * * { potential_file } { ' ' .join (kind_symbols )} \n "
253253 if not default_potential [potential .pair_style ].get ("read_from_file" ):
254254 data = [
255255 line
256256 for line in potential .get_content ().split ("\n " )
257257 if not line .startswith ("#" ) and line
258258 ]
259- potential_block += f' pair_coeff * * { " " .join (data )} \n '
259+ potential_block += f" pair_coeff * * { ' ' .join (data )} \n "
260260
261261 if "neighbor" in parameters_potential :
262262 potential_block += (
263- f' neighbor { join_keywords (parameters_potential [" neighbor" ])} \n '
263+ f" neighbor { join_keywords (parameters_potential [' neighbor' ])} \n "
264264 )
265265 if "neighbor_modify" in parameters_potential :
266266 potential_block += "neigh_modify"
267267 potential_block += (
268- f' { join_keywords (parameters_potential [" neighbor_modify" ])} \n '
268+ f" { join_keywords (parameters_potential [' neighbor_modify' ])} \n "
269269 )
270270 potential_block += generate_header ("End of Potential information" )
271271 return potential_block
@@ -304,7 +304,7 @@ def write_structure_block(
304304 kind_name_id_map [site .kind_name ] = len (kind_name_id_map ) + 1
305305
306306 structure_block = generate_header ("Start of the Structure information" )
307- structure_block += f' box tilt { parameters_structure .get (" box_tilt" , " small" )} \n '
307+ structure_block += f" box tilt { parameters_structure .get (' box_tilt' , ' small' )} \n "
308308
309309 # Set the dimensions of the structure
310310 if "dimension" in parameters_structure :
@@ -319,12 +319,12 @@ def write_structure_block(
319319 structure_block += f"boundary { ' ' .join (['p' if entry else 'f' for entry in structure .pbc ])} \n "
320320
321321 # Set the atom style for the structure
322- structure_block += f' atom_style { parameters_structure [" atom_style" ]} \n '
322+ structure_block += f" atom_style { parameters_structure [' atom_style' ]} \n "
323323
324324 # Set the atom modify for the structure
325325 if "atom_modify" in parameters_structure :
326326 structure_block += "atom_modify"
327- structure_block += f' { parameters_structure [" atom_modify" ]} \n '
327+ structure_block += f" { parameters_structure [' atom_modify' ]} \n "
328328 # Write the command to read the structure from a file
329329 structure_block += f"read_data { structure_filename } \n "
330330 # Set the groups which will be used for the calculations
@@ -338,7 +338,7 @@ def write_structure_block(
338338 raise ValueError ("atom type not defined" )
339339 # Set the current group
340340 structure_block += (
341- f' group { _group [" name" ]} { join_keywords (_group [" args" ])} \n '
341+ f" group { _group [' name' ]} { join_keywords (_group [' args' ])} \n "
342342 )
343343 # Store the name of the group for later usage
344344 group_names .append (_group ["name" ])
@@ -362,11 +362,11 @@ def write_minimize_block(parameters_minimize: dict[str, Union[str, float, int]])
362362 """
363363
364364 minimize_block = generate_header ("Start of the Minimization information" )
365- minimize_block += f' min_style { parameters_minimize .get (" style" , "cg" )} \n '
366- minimize_block += f' minimize { parameters_minimize .get (" energy_tolerance" , 1e-4 )} '
367- minimize_block += f' { parameters_minimize .get (" force_tolerance" , 1e-4 )} '
368- minimize_block += f' { parameters_minimize .get (" max_iterations" , 1000 )} '
369- minimize_block += f' { parameters_minimize .get (" max_evaluations" , 1000 )} \n '
365+ minimize_block += f" min_style { parameters_minimize .get (' style' , 'cg' )} \n "
366+ minimize_block += f" minimize { parameters_minimize .get (' energy_tolerance' , 1e-4 )} "
367+ minimize_block += f" { parameters_minimize .get (' force_tolerance' , 1e-4 )} "
368+ minimize_block += f" { parameters_minimize .get (' max_iterations' , 1000 )} "
369+ minimize_block += f" { parameters_minimize .get (' max_evaluations' , 1000 )} \n "
370370 minimize_block += generate_header ("End of the Minimization information" )
371371
372372 return minimize_block
@@ -403,18 +403,18 @@ def write_md_block(parameters_md: dict[str, Any]) -> str:
403403
404404 md_block = generate_header ("Start of the MD information" )
405405 _key = parameters_md ["integration" ].get ("style" , "nve" )
406- md_block += f' fix { generate_id_tag (_key , " all" )} all { _key } { integration_options } \n '
406+ md_block += f" fix { generate_id_tag (_key , ' all' )} all { _key } { integration_options } \n "
407407 if "velocity" in parameters_md :
408408 md_block += (
409- f' { generate_velocity_string (parameters_velocity = parameters_md [" velocity" ])} '
409+ f" { generate_velocity_string (parameters_velocity = parameters_md [' velocity' ])} "
410410 )
411411 md_block += f"reset_timestep { time_step } \n "
412412 if parameters_md .get ("run_style" , "verlet" ) == "respa" :
413- md_block += f' run_style { parameters_md .get (" run_style" , " verlet" )} '
414- md_block += f' { join_keywords (parameters_md [" respa_options" ])} \n '
413+ md_block += f" run_style { parameters_md .get (' run_style' , ' verlet' )} "
414+ md_block += f" { join_keywords (parameters_md [' respa_options' ])} \n "
415415 else :
416- md_block += f' run_style { parameters_md .get (" run_style" , " verlet" )} \n '
417- md_block += f' run { parameters_md .get (" max_number_steps" , 100 )} \n '
416+ md_block += f" run_style { parameters_md .get (' run_style' , ' verlet' )} \n "
417+ md_block += f" run { parameters_md .get (' max_number_steps' , 100 )} \n "
418418 md_block += generate_header ("End of the MD information" )
419419
420420 return md_block
@@ -476,29 +476,29 @@ def generate_velocity_string(parameters_velocity: list[dict[str, Any]]) -> str:
476476 for entry in parameters_velocity :
477477 _options = generate_velocity_options (entry )
478478 if "create" in entry :
479- options += f' velocity { entry .get (" group" , " all" )} create'
480- options += f' { entry [" create" ].get (" temp" ) } '
479+ options += f" velocity { entry .get (' group' , ' all' )} create"
480+ options += f" { entry [' create' ].get (' temp' ) } "
481481 options += (
482- f' { entry [" create" ].get (" seed" , np .random .randint (10000 ))} { _options } \n '
482+ f" { entry [' create' ].get (' seed' , np .random .randint (10000 ))} { _options } \n "
483483 )
484484 if "set" in entry :
485- options += f' velocity { entry .get (" group" , " all" )} set'
486- options += f' { entry [" set" ].get ("vx" , " NULL" ) } '
487- options += f' { entry [" set" ].get ("vy" , " NULL" ) } '
488- options += f' { entry [" set" ].get ("vz" , " NULL" )} { _options } \n '
485+ options += f" velocity { entry .get (' group' , ' all' )} set"
486+ options += f" { entry [' set' ].get ('vx' , ' NULL' ) } "
487+ options += f" { entry [' set' ].get ('vy' , ' NULL' ) } "
488+ options += f" { entry [' set' ].get ('vz' , ' NULL' )} { _options } \n "
489489 if "scale" in entry :
490- options += f' velocity { entry .get (" group" , " all" )} scale'
491- options += f' { entry [" scale" ]} { _options } \n '
490+ options += f" velocity { entry .get (' group' , ' all' )} scale"
491+ options += f" { entry [' scale' ]} { _options } \n "
492492 if "ramp" in entry :
493- options += f' velocity { entry .get (" group" , " all" )} ramp'
494- options += f' { entry [" ramp" ].get (" vdim" )} { entry [" ramp" ].get (" vlo" ) } '
495- options += f' { entry [" ramp" ].get (" vhi" )} { entry [" ramp" ].get (" dim" ) } '
493+ options += f" velocity { entry .get (' group' , ' all' )} ramp"
494+ options += f" { entry [' ramp' ].get (' vdim' )} { entry [' ramp' ].get (' vlo' ) } "
495+ options += f" { entry [' ramp' ].get (' vhi' )} { entry [' ramp' ].get (' dim' ) } "
496496 options += (
497- f' { entry [" ramp" ].get (" clo" )} { entry [" ramp" ].get (" chi" )} { _options } \n '
497+ f" { entry [' ramp' ].get (' clo' )} { entry [' ramp' ].get (' chi' )} { _options } \n "
498498 )
499499 if "zero" in entry :
500- options += f' velocity { entry .get (" group" , " all" )} zero'
501- options += f' { entry [" zero" ]} { _options } \n '
500+ options += f" velocity { entry .get (' group' , ' all' )} zero"
501+ options += f" { entry [' zero' ]} { _options } \n "
502502 return options
503503
504504
@@ -618,31 +618,31 @@ def generate_integration_options(
618618 _value = integration_parameters .get (_option )
619619 if _value :
620620 _value = [str (val ) for val in _value ]
621- options += f' { _option } { " " .join (_value ) if isinstance (_value , list ) else _value } '
621+ options += f" { _option } { ' ' .join (_value ) if isinstance (_value , list ) else _value } "
622622 # Set the options that depend on the pressure
623623 if style in pressure_dependent :
624624 for _option in pressure_options :
625625 if _option in integration_parameters :
626626 _value = integration_parameters .get (_option )
627627 if _value :
628628 _value = [str (val ) for val in _value ]
629- options += f' { _option } { " " .join (_value ) if isinstance (_value , list ) else _value } '
629+ options += f" { _option } { ' ' .join (_value ) if isinstance (_value , list ) else _value } "
630630 # Set the options that depend on the 'uef' parameters
631631 if style in uef_dependent :
632632 for _option in uef_options :
633633 if _option in integration_parameters :
634634 _value = integration_parameters .get (_option )
635635 if _value :
636636 _value = [str (val ) for val in _value ]
637- options += f' { _option } { " " .join (_value ) if isinstance (_value , list ) else _value } '
637+ options += f" { _option } { ' ' .join (_value ) if isinstance (_value , list ) else _value } "
638638 # Set the options that depend on the 'nve/limit' parameters
639639 if style in ["nve/limit" ]:
640- options += f' { integration_parameters .get (" xmax" , 0.1 )} '
640+ options += f" { integration_parameters .get (' xmax' , 0.1 )} "
641641 # Set the options that depend on the 'langevin' parameters
642642 if style in ["nve/dotc/langevin" ]:
643- options += f' { integration_parameters .get (" temp" ) } '
644- options += f' { integration_parameters .get (" seed" ) } '
645- options += f' angmom { integration_parameters .get (" angmom" ) } '
643+ options += f" { integration_parameters .get (' temp' ) } "
644+ options += f" { integration_parameters .get (' seed' ) } "
645+ options += f" angmom { integration_parameters .get (' angmom' ) } "
646646 return options
647647
648648
@@ -685,7 +685,7 @@ def write_fix_block(
685685 f'group name "{ _group } " is not the defined groups { [* group_names , "all" ]} '
686686 )
687687 fix_block += f"fix { generate_id_tag (key , _group )} { _group } { key } "
688- fix_block += f' { join_keywords (entry [" type" ])} \n '
688+ fix_block += f" { join_keywords (entry [' type' ])} \n "
689689 fix_block += generate_header ("End of the Fix information" )
690690 return fix_block
691691
@@ -721,7 +721,7 @@ def write_compute_block(
721721 if _group not in [* group_names , "all" ]:
722722 raise ValueError (f'group name "{ _group } " is not the defined groups' )
723723 compute_block += f"compute { generate_id_tag (key , _group )} { _group } { key } "
724- compute_block += f' { join_keywords (entry [" type" ])} \n '
724+ compute_block += f" { join_keywords (entry [' type' ])} \n "
725725 compute_block += generate_header ("End of the Compute information" )
726726 return compute_block
727727
@@ -780,12 +780,12 @@ def write_dump_block(
780780 if atom_style == "charge" :
781781 num_double += 1
782782 dump_block = generate_header ("Start of the Dump information" )
783- dump_block += f' dump aiida all custom { parameters_dump .get (" dump_rate" , 10 )} '
783+ dump_block += f" dump aiida all custom { parameters_dump .get (' dump_rate' , 10 )} "
784784 dump_block += f"{ trajectory_filename } id type element x y z "
785- dump_block += f' { "q " if atom_style == " charge" else "" } '
786- dump_block += f' { " " .join (computes_list )} \n '
785+ dump_block += f" { 'q ' if atom_style == ' charge' else '' } "
786+ dump_block += f" { ' ' .join (computes_list )} \n "
787787 dump_block += "dump_modify aiida sort id\n "
788- dump_block += f' dump_modify aiida element { " " .join (kind_symbols )} \n '
788+ dump_block += f" dump_modify aiida element { ' ' .join (kind_symbols )} \n "
789789 dump_block += "dump_modify aiida format int ' %d ' \n "
790790 dump_block += "dump_modify aiida format float ' %16.10e ' \n "
791791 dump_block += generate_header ("End of the Dump information" )
@@ -855,9 +855,9 @@ def write_thermo_block(
855855
856856 thermo_block = generate_header ("Start of the Thermo information" )
857857 thermo_block += (
858- f' thermo_style custom { " " .join (fixed_thermo )} { " " .join (computes_list )} \n '
858+ f" thermo_style custom { ' ' .join (fixed_thermo )} { ' ' .join (computes_list )} \n "
859859 )
860- thermo_block += f' thermo { parameters_thermo .get (" printing_rate" , 1000 )} \n '
860+ thermo_block += f" thermo { parameters_thermo .get (' printing_rate' , 1000 )} \n "
861861 thermo_block += generate_header ("End of the Thermo information" )
862862
863863 printing_variables = fixed_thermo + list (
@@ -898,7 +898,7 @@ def write_restart_block(
898898 "Start of the intermediate write restart information"
899899 )
900900 restart_block ["intermediate" ] += (
901- f"restart { parameters_restart .get ('num_steps' , int (max_number_steps / 10 ))} { restart_filename } \n "
901+ f"restart { parameters_restart .get ('num_steps' , int (max_number_steps / 10 ))} { restart_filename } \n "
902902 )
903903 restart_block ["intermediate" ] += generate_header (
904904 "End of the intermediate write restart information"
@@ -1000,7 +1000,7 @@ def generate_id_tag(name: str, group: str) -> str:
10001000 :rtype: str
10011001 """
10021002
1003- return f"{ name .replace ('/' ,'_' )} _{ group } _aiida"
1003+ return f"{ name .replace ('/' , '_' )} _{ group } _aiida"
10041004
10051005
10061006def join_keywords (value : list [Any ]) -> str :
0 commit comments