diff --git a/bfps/FluidConvert.py b/bfps/FluidConvert.py index cbb9d9c3f478463426afbb63c18bb26abba49285..9a06f2e9dfdb80a21bb86a0f08db508567862a06 100644 --- a/bfps/FluidConvert.py +++ b/bfps/FluidConvert.py @@ -37,14 +37,14 @@ class FluidConvert(_fluid_particle_base): """ def __init__( self, - name = 'FluidConvert', + name = 'FluidConvert-v' + bfps.__version__, work_dir = './', simname = 'test', fluid_precision = 'single', use_fftw_wisdom = True): _fluid_particle_base.__init__( self, - name = name, + name = name + '-' + fluid_precision, work_dir = work_dir, simname = simname, dtype = fluid_precision, diff --git a/bfps/FluidResize.py b/bfps/FluidResize.py index 8e03e1657e2f1551a4e27e8442dd0fc1032bc876..057b40f6876ff3d79931c67ea45a99a70371c87b 100644 --- a/bfps/FluidResize.py +++ b/bfps/FluidResize.py @@ -38,14 +38,14 @@ class FluidResize(_fluid_particle_base): """ def __init__( self, - name = 'FluidResize', + name = 'FluidResize-v' + bfps.__version__, work_dir = './', simname = 'test', dtype = np.float32, use_fftw_wisdom = False): _fluid_particle_base.__init__( self, - name = name, + name = name + '-' + fluid_precision, work_dir = work_dir, simname = simname, dtype = dtype, diff --git a/bfps/NavierStokes.py b/bfps/NavierStokes.py index 5f4cbc90e40964d2ba06b331c5e33d1a8f170ab4..e32f53123354cc5f850cb99f02b10d129b8bebf5 100644 --- a/bfps/NavierStokes.py +++ b/bfps/NavierStokes.py @@ -41,7 +41,7 @@ class NavierStokes(_fluid_particle_base): """ def __init__( self, - name = 'NavierStokes', + name = 'NavierStokes-v' + bfps.__version__, work_dir = './', simname = 'test', fluid_precision = 'single', @@ -54,7 +54,7 @@ class NavierStokes(_fluid_particle_base): self.fftw_plan_rigor = fftw_plan_rigor _fluid_particle_base.__init__( self, - name = name, + name = name + '-' + fluid_precision, work_dir = work_dir, simname = simname, dtype = fluid_precision, @@ -62,8 +62,8 @@ class NavierStokes(_fluid_particle_base): self.parameters['nu'] = 0.1 self.parameters['fmode'] = 1 self.parameters['famplitude'] = 0.5 - self.parameters['fk0'] = 1.5 - self.parameters['fk1'] = 3.0 + self.parameters['fk0'] = 2.0 + self.parameters['fk1'] = 4.0 self.parameters['forcing_type'] = 'linear' self.parameters['histogram_bins'] = 256 self.parameters['max_velocity_estimate'] = 1.0 @@ -984,6 +984,11 @@ class NavierStokes(_fluid_particle_base): self.QR_stats_on = opt.QR_stats self.parameters['nu'] = (opt.kMeta * 2 / opt.n)**(4./3) self.parameters['dt'] = (opt.dtfactor / opt.n) + # custom famplitude for 288 and 576 + if opt.n == 288: + self.parameters['famplitude'] = 0.45 + elif opt.n == 576: + self.parameters['famplitude'] = 0.47 if ((self.parameters['niter_todo'] % self.parameters['niter_out']) != 0): self.parameters['niter_out'] = self.parameters['niter_todo'] if self.QR_stats_on: @@ -993,6 +998,7 @@ class NavierStokes(_fluid_particle_base): meantrS2 = (opt.n//2 / opt.kMeta)**4 * self.parameters['nu']**2 self.parameters['max_Q_estimate'] = meantrS2 self.parameters['max_R_estimate'] = .4*meantrS2**1.5 + self.name += 'QR', if len(opt.src_work_dir) == 0: opt.src_work_dir = opt.work_dir diff --git a/bfps/_code.py b/bfps/_code.py index 755c51d6d948b4d2845a8b19ca007cfbf98f4d12..3589689cebd20104dd5bc034ecd6bcba63a2df97 100644 --- a/bfps/_code.py +++ b/bfps/_code.py @@ -198,7 +198,7 @@ class _code(_base): self.write_sge_file( file_name = os.path.join(self.work_dir, qsub_script_name), nprocesses = ncpu, - name_of_run = self.name + '_' + suffix, + name_of_run = suffix, command_atoms = command_atoms[3:], hours = hours, minutes = minutes, @@ -210,7 +210,7 @@ class _code(_base): qsub_atoms += ['-hold_jid', job_name_list[-1]] subprocess.call(qsub_atoms + [qsub_script_name]) os.chdir(current_dir) - job_name_list.append(self.name + '_' + suffix) + job_name_list.append(suffix) elif self.host_info['type'] == 'pc': os.chdir(self.work_dir) os.environ['LD_LIBRARY_PATH'] += ':{0}'.format(bfps.lib_dir) diff --git a/bfps/_fluid_base.py b/bfps/_fluid_base.py index d82f83c3a71b3fa8dbfbd77b2d855c0b941de5ee..82e073d2c04959f9f418ebc98ed02c168a1aeda5 100644 --- a/bfps/_fluid_base.py +++ b/bfps/_fluid_base.py @@ -48,7 +48,7 @@ class _fluid_particle_base(_code): work_dir = work_dir, simname = simname) self.use_fftw_wisdom = use_fftw_wisdom - self.name = name + '_' + simname + self.name = name self.particle_species = 0 if dtype in [np.float32, np.float64]: self.dtype = dtype