Yun Wang

Does `dDynamicArray` allow us to remove multiple elements at the same time?

2

Hi, @leonfoks I'm testing the functionalities of ```idDynamicArray_Class```. And I know that I can remove one element of an array by using ```idDynamicArray%remove(index)```, where ```index``` is an integer. What I...

Failed to compile the library in Linux with gfortran version using gfortran version 6.3.1

4

Hello, @leonfoks Following your instructions, I typed ```cmake ../src``` in the command line, and everything is OK. Then I typed ```make``` to compile the library. However, it failed and returned...

Issues about the module "m_random"

4

Hi, I found three possible issues when I tested "coretran" with some simple programs. 1. Outputs of ```rngInteger``` are constants For testing the functionalities of ```rngInteger```, I write the code...

"Overflow Error: value too large to convert to int" When I used "ScatterArray()" to scatter large numpy arrays

1

Hi, @rainwoodman For performing 'FOF' parallel, I want to load all coordinates of 1820**3 particles if `comm.rank == 0`. And then scatter these coordinates to other ranks. My code snippet...

Hi, The snapshot data of IllustrisTNG at a given redshift is split into many hdf5 chunks. For example, the snapshot at redshift 1.0 was split into 448 chunks, from 'snapshot-99.0.hdf5'...

How to convert particles' velocities to the gridded velocity field?

4

Hello, I have 1820^3 particles and their spatial velocities ( unit: km/s, shape: (1820^3, 3) ). Then I want to convert these 'scattered' velocity data to the 1024^3 mesh data,...

Is there any way to parallelize the computation of multiple convolutions?

8

Hi, First, I have a 3D density mesh (Nmesh=1536^3) obtained by using nbodykit, stored as a bigfile on my disk. Then I want to convolve the mesh with a band-pass...

```HDFCatalog.to_mesh``` generate negative values for the density field 'delta+1'.

6

By using ```HDFCatalog.to_mesh```, I want to generate a mesh of 1024^3 grids from the discrete particles of TNG100 simulation. ``` in_mesh = HDFCatalog.to_mesh(Nmesh=1024,BoxSize=75,dtype='f8', interlaced=True, compensated=True, resampler='cic', position='Coordinates',weight='Masses') ``` However, I...

How to load array from the disk, such as in the format of '.npy' ?

1

Hi, I find that mpi4py-fft is an excellent tool for performing the FFT in parallel, which greatly improves computational efficiency. However, I don't see anything in the documentation about how...