Timings for --r=0 and --r=1/2 for LR and RBF

10_Timings.py

@@ -67,10 +67,10 @@ if args.Architecture is not None:
         raise ValueError('Architecture not recognized')
 
 models = {
-    'LinearRegression': LinearRegressionModel(H_REFINEMENTS=H_REFINEMENTS),
-    'RBF_Linear': RadialBasisRegressionModel(H_REFINEMENTS=H_REFINEMENTS),
-    'RBF_thinplatespline': RadialBasisRegressionModel(H_REFINEMENTS=H_REFINEMENTS),
-    'RBF_cubic': RadialBasisRegressionModel(H_REFINEMENTS=H_REFINEMENTS),
+    # 'LinearRegression': LinearRegressionModel(H_REFINEMENTS=H_REFINEMENTS),
+    # 'RBF_Linear': RadialBasisRegressionModel(H_REFINEMENTS=H_REFINEMENTS),
+    # 'RBF_thinplatespline': RadialBasisRegressionModel(H_REFINEMENTS=H_REFINEMENTS),
+    # 'RBF_cubic': RadialBasisRegressionModel(H_REFINEMENTS=H_REFINEMENTS),
     # 'RBF_quintic': RadialBasisRegressionModel(H_REFINEMENTS=H_REFINEMENTS),
     # 'RBF_multiquadric': RadialBasisRegressionModel(H_REFINEMENTS=H_REFINEMENTS),
     # 'RBF_inversemultiquadric': RadialBasisRegressionModel(H_REFINEMENTS=H_REFINEMENTS),
@@ -115,6 +115,8 @@ else:
 # Loop over the models
 for ModelIndex, model_name in enumerate(models):
     print(f'Model: \t\t\t{model_name}')
+    if not os.path.exists(f'temp_{model_name}'):
+        os.mkdir(f'temp_{model_name}')
 
     model = models[model_name]
 
@@ -181,7 +183,7 @@ for ModelIndex, model_name in enumerate(models):
                 ).communicate()[0]
             end_reconstruct = time.time()
             timings_reconstruct[RunIndex, PODIndex] = end_reconstruct - start_reconstruct
-            
+
             os.remove('pressure_field_patches.xml')
             os.remove('velocity_field_patches.xml')
             os.remove('pressure_field_rec_patches.xml')
@@ -190,7 +192,7 @@ for ModelIndex, model_name in enumerate(models):
     shutil.rmtree(f'temp_{model_name}')
 
     # Save the timings
-    np.savez(f'{TimingsFolder}/{model_name}.npz', timings_train_w_POD=timings_train_w_POD, timings_train_wo_POD=timings_train_wo_POD, timings_evaluate=timings_evaluate)
+    np.savez(f'{TimingsFolder}/{model_name}.npz', timings_train_w_POD=timings_train_w_POD, timings_train_wo_POD=timings_train_wo_POD, timings_evaluate=timings_evaluate, timings_reconstruct=timings_reconstruct)
 
     # Update the counter for the kernel
     if type(model) == RadialBasisRegressionModel: