 |
Peano
|
|
Peano
|
Variables | |
| int | numPart = 65 |
| float | gamma = 5.0 / 3.0 |
| float | rho0 = 1.0 |
| float | P0 = 1.0e-6 |
| float | E0 = 1.0 |
| int | N_inject = 3 |
| str | fileName = "sedov.hdf5" |
| coords = zeros((numPart, 3)) | |
| h = zeros(numPart) | |
| float | vol = 1.0 |
| v = zeros((numPart, 3)) | |
| ids = linspace(1, numPart, numPart) | |
| m = zeros(numPart) | |
| u = zeros(numPart) | |
| r = zeros(numPart) | |
| index = argsort(r) | |
| file = h5py.File(fileName, "w") | |
| grp = file.create_group("/Header") | |
| data | |
| dtype | |
| int | numPart_L = 800 |
| float | x_min = -1.0 |
| float | x_max = 1.0 |
| float | rho_L = 1.0 |
| float | rho_R = 0.125 |
| float | v_L = 0.0 |
| float | v_R = 0.0 |
| float | P_L = 1.0 |
| float | P_R = 0.1 |
| float | boxSize = x_max - x_min |
| numPart_R = int(numPart_L * (rho_R / rho_L)) | |
| tuple | delta_L = (boxSize / 2) / numPart_L |
| tuple | delta_R = (boxSize / 2) / numPart_R |
| tuple | offset_L = delta_L / 2 |
| tuple | offset_R = delta_R / 2 |
| int | i = numPart_L + j |
| parser = argparse.ArgumentParser(description="SPH benchmarking script") | |
| dest | |
| required | |
| type | |
| help | |
| default | |
| args = parser.parse_args() | |
| npart = args.particle_number | |
| seed = args.random_seed | |
| max_displacement = args.max_displacement | |
| float | dx = boxSize / npart |
| rng = np.random.RandomState(seed=seed) | |
| tuple | displacement = (2.0 * rng.random(coords.shape) - 1.0) * max_displacement * dx |
| action | |
| int | ind = 0 |
| tuple | x = (i + 0.5) * dx |
| tuple | y = (j + 0.5) * dx |
| float | scale = 0.8 |
| s | |
| marker | |
| tuple | z = (k + 0.5) * dx |
| tuple makeIC.displacement = (2.0 * rng.random(coords.shape) - 1.0) * max_displacement * dx |
| makeIC.parser = argparse.ArgumentParser(description="SPH benchmarking script") |
1.10.0