fix: missing bracket around multiple variables

GPUSimulators/Common.py

@@ -35,8 +35,6 @@ import gc
 import netCDF4
 import json
 
-from tqdm import trange
-
 #import pycuda.compiler as cuda_compiler
 #import pycuda.gpuarray
 #import pycuda.driver as cuda
@@ -180,8 +178,8 @@ def runSimulation(simulator, simulator_args, outfile, save_times, save_var_names
         profiling_data_sim_runner["end"]["t_sim_init"] = time.time()
 
         #Start simulation loop
-        # progress_printer = ProgressPrinter(save_times[-1], print_every=10)
-        for k in trange(len(save_times)):
+        progress_printer = ProgressPrinter(save_times[-1], print_every=10)
+        for k in range(len(save_times)):
             #Get target time and step size there
             t_step = t_steps[k]
             t_end = save_times[k]
@@ -213,9 +211,9 @@ def runSimulation(simulator, simulator_args, outfile, save_times, save_var_names
             profiling_data_sim_runner["end"]["t_nc_write"] += time.time()
 
             #Write progress to screen
-            # print_string = progress_printer.getPrintString(t_end)
-            # if (print_string):
-            #     logger.debug(print_string)
+            print_string = progress_printer.getPrintString(t_end)
+            if (print_string):
+                logger.debug(print_string)
                 
         logger.debug("Simulated to t={:f} in {:d} timesteps (average dt={:f})".format(t_end, sim.simSteps(), sim.simTime() / sim.simSteps()))
         
@@ -435,58 +433,58 @@ class DataDumper(object):
         
 
 
-# class ProgressPrinter(object):
-#     """
-#     Small helper class for 
-#     """
-#     def __init__(self, total_steps, print_every=5):
-#         self.logger = logging.getLogger(__name__)
-#         self.start = time.time()
-#         self.total_steps = total_steps
-#         self.print_every = print_every
-#         self.next_print_time = self.print_every
-#         self.last_step = 0
-#         self.secs_per_iter = None
+class ProgressPrinter(object):
+    """
+    Small helper class for 
+    """
+    def __init__(self, total_steps, print_every=5):
+        self.logger = logging.getLogger(__name__)
+        self.start = time.time()
+        self.total_steps = total_steps
+        self.print_every = print_every
+        self.next_print_time = self.print_every
+        self.last_step = 0
+        self.secs_per_iter = None
         
-#     def getPrintString(self, step):
-#         elapsed =  time.time() - self.start
-#         if (elapsed > self.next_print_time):            
-#             dt = elapsed - (self.next_print_time - self.print_every)
-#             dsteps = step - self.last_step
-#             steps_remaining = self.total_steps - step
+    def getPrintString(self, step):
+        elapsed =  time.time() - self.start
+        if (elapsed > self.next_print_time):            
+            dt = elapsed - (self.next_print_time - self.print_every)
+            dsteps = step - self.last_step
+            steps_remaining = self.total_steps - step
                         
-#             if (dsteps == 0):
-#                 return
+            if (dsteps == 0):
+                return
                 
-#             self.last_step = step
-#             self.next_print_time = elapsed + self.print_every
+            self.last_step = step
+            self.next_print_time = elapsed + self.print_every
             
-#             if not self.secs_per_iter:
-#                 self.secs_per_iter = dt / dsteps
-#             self.secs_per_iter = 0.2*self.secs_per_iter + 0.8*(dt / dsteps)
+            if not self.secs_per_iter:
+                self.secs_per_iter = dt / dsteps
+            self.secs_per_iter = 0.2*self.secs_per_iter + 0.8*(dt / dsteps)
             
-#             remaining_time = steps_remaining * self.secs_per_iter
+            remaining_time = steps_remaining * self.secs_per_iter
 
-#             return "{:s}. Total: {:s}, elapsed: {:s}, remaining: {:s}".format(
-#                 ProgressPrinter.progressBar(step, self.total_steps), 
-#                 ProgressPrinter.timeString(elapsed + remaining_time), 
-#                 ProgressPrinter.timeString(elapsed), 
-#                 ProgressPrinter.timeString(remaining_time))
+            return "{:s}. Total: {:s}, elapsed: {:s}, remaining: {:s}".format(
+                ProgressPrinter.progressBar(step, self.total_steps), 
+                ProgressPrinter.timeString(elapsed + remaining_time), 
+                ProgressPrinter.timeString(elapsed), 
+                ProgressPrinter.timeString(remaining_time))
 
-#     def timeString(seconds):
-#         seconds = int(max(seconds, 1))
-#         minutes, seconds = divmod(seconds, 60)
-#         hours, minutes = divmod(minutes, 60)
-#         periods = [('h', hours), ('m', minutes), ('s', seconds)]
-#         time_string = ' '.join('{}{}'.format(value, name)
-#                                 for name, value in periods
-#                                 if value)
-#         return time_string
+    def timeString(seconds):
+        seconds = int(max(seconds, 1))
+        minutes, seconds = divmod(seconds, 60)
+        hours, minutes = divmod(minutes, 60)
+        periods = [('h', hours), ('m', minutes), ('s', seconds)]
+        time_string = ' '.join('{}{}'.format(value, name)
+                                for name, value in periods
+                                if value)
+        return time_string
 
-#     def progressBar(step, total_steps, width=30):
-#         progress = np.round(width * step / total_steps).astype(np.int32)
-#         progressbar = "0% [" + "#"*(progress) + "="*(width-progress) + "] 100%"
-#         return progressbar
+    def progressBar(step, total_steps, width=30):
+        progress = np.round(width * step / total_steps).astype(np.int32)
+        progressbar = "0% [" + "#"*(progress) + "="*(width-progress) + "] 100%"
+        return progressbar
 
 
 """

GPUSimulators/Simulator.py

@@ -25,7 +25,6 @@ import numpy as np
 import math
 import logging
 from enum import IntEnum
-from tqdm import tqdm
 
 #import pycuda.compiler as cuda_compiler
 #import pycuda.gpuarray
@@ -157,7 +156,7 @@ class BaseSimulator(object):
         self.num_substeps = num_substeps
         
         #Handle autotuning block size
-        if self.context.autotuner:
+        if (self.context.autotuner):
             peak_configuration = self.context.autotuner.get_peak_performance(self.__class__)
             block_width = int(peak_configuration["block_width"])
             block_height = int(peak_configuration["block_height"])
@@ -196,45 +195,42 @@ class BaseSimulator(object):
         Requires that the step() function is implemented in the subclasses
         """
 
-        # printer = Common.ProgressPrinter(t)
+        printer = Common.ProgressPrinter(t)
         
         t_start = self.simTime()
         t_end = t_start + t
         
         update_dt = True
-        if dt is not None:
+        if (dt is not None):
             update_dt = False
             self.dt = dt
         
-        for _ in tqdm(range(math.ceil(t_end / self.dt))):
+        while(self.simTime() < t_end):
             # Update dt every 100 timesteps and cross your fingers it works
             # for the next 100
-            # TODO this is probably broken now after fixing the "infinite" loop
-            if update_dt and (self.simSteps() % 100 == 0):
+            if (update_dt and (self.simSteps() % 100 == 0)):
                 self.dt = self.computeDt()*self.cfl_scale
-
+        
             # Compute timestep for "this" iteration (i.e., shorten last timestep)
             current_dt = np.float32(min(self.dt, t_end-self.simTime()))
 
             # Stop if end reached (should not happen)
-            if current_dt <= 0.0:
+            if (current_dt <= 0.0):
                 self.logger.warning("Timestep size {:d} is less than or equal to zero!".format(self.simSteps()))
                 break
-
+        
             # Step forward in time
             self.step(current_dt)
 
             #Print info
-            # print_string = printer.getPrintString(self.simTime() - t_start)
-            # if (print_string):
-            #     self.logger.info("%s: %s", self, print_string)
-            #     try:
-            #         self.check()
-            #     except AssertionError as e:
-            #         e.args += ("Step={:d}, time={:f}".format(self.simSteps(), self.simTime()),)
-            #         raise
-
-        print("Done")
+            print_string = printer.getPrintString(self.simTime() - t_start)
+            if (print_string):
+                self.logger.info("%s: %s", self, print_string)
+                try:
+                    self.check()
+                except AssertionError as e:
+                    e.args += ("Step={:d}, time={:f}".format(self.simSteps(), self.simTime()),)
+                    raise
 
 
     def step(self, dt):

Jobs/job_lumi.slrum

@@ -22,7 +22,7 @@ Myapplication=${Mydir}/FiniteVolumeGPU_hip/mpiTesting.py
 ml LUMI/24.03 partition/G
 ml lumi-container-wrapper
 ml cray-python/3.11.7
-ml rocm/6.2.2
+ml rocm/5.4.6
 
 ml craype-accel-amd-gfx90a
 ml cray-mpich/8.1.29

conda_environment_lumi.yml

@@ -13,9 +13,8 @@ dependencies:
 - pytools
 - netcdf4
 - scipy
-- tqdm
 - pip:
-  - hip-python==6.2.0.499.16
+  - hip-python==5.4.3.470.16
   - -i https://test.pypi.org/simple/