zachetienne / nrpytutorial Goto Github PK

Code:

import cmdline_helper as cmd
cmd.new_C_compile(Ccodesrootdir, os.path.join("output", "BrillLindquist_Playground"),
                  uses_free_parameters_h=True, compiler_opt_option="fast") # fastdebug or debug also supported

# Change to output directory
os.chdir(os.path.join(Ccodesrootdir, "output"))
# Clean up existing output files
cmd.delete_existing_files("out*.txt")
cmd.delete_existing_files("out*.png")
# Run executable with CFL_FACTOR = 1.0, which is allowed since
#   simulation is axisymmetric and all phi derivs are set to zero.
# cmd.Execute("BrillLindquist_Playground", "72 12 2 "+str(CFL_FACTOR))
# cmd.Execute("BrillLindquist_Playground", "800 16 2 "+str(CFL_FACTOR))
cmd.Execute("BrillLindquist_Playground", str(Nr) + " " + str(Nth) + " " + str(Nph) + " " + str(CFL_FACTOR))
os.chdir(os.path.join("..", ".."))

Then jupyter shows an error:

File D:\Anaconda3_2309\Lib\subprocess.py:1538, in Popen._execute_child(self, args, executable, preexec_fn, close_fds, pass_fds, cwd, env, startupinfo, creationflags, shell, p2cread, p2cwrite, c2pread, c2pwrite, errread, errwrite, unused_restore_signals, unused_gid, unused_gids, unused_uid, unused_umask, unused_start_new_session, unused_process_group)
   1536 # Start the process
   1537 try:
-> 1538     hp, ht, pid, tid = _winapi.CreateProcess(executable, args,
   1539                              # no special security
   1540                              None, None,
   1541                              int(not close_fds),
   1542                              creationflags,
   1543                              env,
   1544                              cwd,
   1545                              startupinfo)
   1546 finally:
   1547     # Child is launched. Close the parent's copy of those pipe
   1548     # handles that only the child should have open.  You need
   (...)
   1551     # pipe will not close when the child process exits and the
   1552     # ReadFile will hang.
   1553     self._close_pipe_fds(p2cread, p2cwrite,
   1554                          c2pread, c2pwrite,
   1555                          errread, errwrite)

OSError: [WinError 193] %1 不是有效的 Win32 应用程序。

https://notebooks.gesis.org/binder/jupyter/user/zachetienne-nrpytutorial-tqha3xtc/notebooks/Tutorial-Start_to_Finish-BSSNCurvilinear-ScalarField_Collapse.ipynb
It tries to download ffmpeg from http://astro.phys.wvu.edu/zetienne/ffmpeg-static-amd64-johnvansickle.tar.xz

Good evening all,

I encountered a possible issue with the OpenMP parallelization of NRPy+/SENR.
I am outputting C code for the BSSN RHSs in the way that it's done in the "Colliding black holes" tutorial. Currently, the OMP parallelization only happens in the exterior loop, over i2. The problem with this is that, in spherical- or cylindrical-like coordinates, this will be the azimuthal direction which, from my experience, does not require that many points (I typically use 16 for what I'm doing) to accurately sample the physical space. Therefore, if one tries to use more threads than the number of points in direction 2, there is no performance gain.

I tried circumventing this issue by using a custom OMP pragma, in the form "#pragma omp parallel for collapse(3)", but if rfm_precompute is enabled, there are lines of code inserted between the nested loops which result in errors, because with the collapse clause only works if the loops are perfectly nested.

Is there a good way of doing this in NRPy+ implemented right now?

Thank you :)

I use NRPyPN to set up my binary black hole simulations. It would be very useful for me to have an estimate of the time to merger given some initial data. This is so that I can have a sense of how long the simulations is going to be, and to set up horizon finders more precisely. Additionally, sometimes I want to run simulations for N orbit and I don't care about the initial separation. It would be lovely if NRPyPN could tell me the initial separation given the number of orbits that I want to simulate.

It runs "!cp ffmpeg-static-amd64-johnvansickle/ffmpeg ~/.local/bin/" , but when i manually extract ffmpeg-static-amd64-johnvansickle.tar.xz ffmpeg-5.0.1-amd64-static directory gets extracted.

In https://github.com/zachetienne/nrpytutorial/blob/master/Tutorial-Solving_the_Scalar_Wave_Equation_with_NumPy.ipynb the range over which to loop is given as for i in range(NGHOSTS, Nx+NGHOSTS) etc. While correct (Nx is the number of interior points) it may be more intuitive to write this as for i in range(NGHOSTS, NGHOSTS+Nx) to provide more of an indication that the range is "from NGHOSTS to Nx points past NGHOSTS" rather than "to NGHOSTS points more than Nx").

I think that there's a minor copy&paste error at line 853 of the file IllinoisGRMHD/doc/Tutorial-IllinoisGRMHD__the_conservative_to_primitive_algorithm.ipynb.

Extract of the current file (lines 852-853):

"          CCTK_REAL K_ppoly_tab     = eos.K_ppoly_tab[polytropic_index];\n",
"          CCTK_REAL Gamma_ppoly_tab = eos.K_ppoly_tab[polytropic_index];\n",

Shouldn't it be

"          CCTK_REAL K_ppoly_tab     = eos.K_ppoly_tab[polytropic_index];\n",
"          CCTK_REAL Gamma_ppoly_tab = eos.Gamma_ppoly_tab[polytropic_index];\n",

?