############################################################################### ######################## Creating the TrackGenerator ###################### ############################################################################### log.py_printf('NORMAL', 'Initializing the track generator...') track_generator = TrackGenerator(geometry, num_azim, track_spacing) track_generator.generateTracks() ############################################################################### ########################### Running a Simulation ########################## ############################################################################### solver = CPUSolver(geometry, track_generator) solver.setSourceConvergenceThreshold(tolerance) solver.setNumThreads(num_threads) solver.convergeSource(max_iters) solver.printTimerReport() ############################################################################### ############################ Generating Plots ############################# ############################################################################### log.py_printf('NORMAL', 'Plotting data...') plotter.plot_cells(geometry) plotter.plot_materials(geometry) plotter.plot_flat_source_regions(geometry) plotter.plot_fluxes(geometry, solver, energy_groups=[1,2,3,4,5,6,7,8])
############################################################################### ######################## Creating the TrackGenerator ###################### ############################################################################### log.py_printf('NORMAL', 'Initializing the track generator...') track_generator = ModularTrackGenerator(geometry, num_azim, track_spacing) track_generator.setLatticeStructure(25,25) track_generator.generateTracks() plotter.plot_flat_source_regions(geometry, gridsize=1000) ############################################################################### ########################### Running a Simulation ########################## ############################################################################### solver = ModularCPUSolver(geometry, track_generator) solver.setSourceConvergenceThreshold(tolerance) solver.setNumThreads(num_threads) solver.convergeSource(max_iters) solver.printTimerReport() ############################################################################### ############################ Generating Plots ############################# ############################################################################### log.py_printf('NORMAL', 'Plotting data...') plotter.plot_fluxes(geometry, solver, energy_groups=[1,2,3,4,5,6,7,8], gridsize=1000)
try: import openmoc.cuda as cuda # Check whether or not the machine contains a GPU if cuda.machineContainsGPU(): solver = cuda.GPUSolver(geometry, track_generator) else: solver = ThreadPrivateSolver(geometry, track_generator) solver.setNumThreads(num_threads) # If the openmoc.cuda module has not been built, use the standard solver except: solver = ThreadPrivateSolver(geometry, track_generator, cmfd) solver.setNumThreads(num_threads) solver.setSourceConvergenceThreshold(tolerance) solver.convergeSource(max_iters) solver.printTimerReport() ############################################################################### ############################## Plotting Data ############################## ############################################################################### log.py_printf('NORMAL', 'Plotting data...') plotter.plot_fluxes(geometry, solver, energy_groups=[1,7], gridsize=400) log.py_printf('TITLE', 'Finished')
log.py_printf('NORMAL', 'Initializing the track generator...') track_generator = TrackGenerator(geometry, num_azim, track_spacing) track_generator.setNumThreads(num_threads) track_generator.generateTracks() ############################################################################### ########################### Running a Simulation ########################## ############################################################################### solver = CPUSolver(geometry, track_generator) solver.setSourceConvergenceThreshold(tolerance) solver.setNumThreads(num_threads) solver.convergeSource(max_iters) solver.printTimerReport() ############################################################################### ############################ Generating Plots ############################# ############################################################################### log.py_printf('NORMAL', 'Plotting data...') plotter.plot_materials(geometry, gridsize=500) plotter.plot_cells(geometry, gridsize=500) plotter.plot_flat_source_regions(geometry, gridsize=500) plotter.plot_fluxes(geometry, solver, energy_groups=[1,2,3,4,5,6,7]) log.py_printf('TITLE', 'Finished')