Newest Questions
                11,828 questions
            
            
            
                3
            
            votes
        
        
            
                0
            
            answers
        
        
            
                36
            
            views
        
        
            
            
        AMD Epyc F or X variant for MODFLOW?
                    I'm trying to choose a processor type for an HPC cluster for calculations with MODFLOW 6 (extended version). I'm looking at the F and X variants of AMD Epyc processors. The F variant has a higher ...
                
            
       
        
            
                0
            
            votes
        
        
            
                1
            
            answer
        
        
            
                52
            
            views
        
        
            
            
        Help implementing the Hessenberg QR algorithm using the Keras Ops API
                    I am trying to implement the Complex and Real Schur Decompositions using only the Keras Ops API, so I can utilize it to impose some eigenvalue constraints on a weight matrix. I am also using the JAX ...
                
            
       
        
            
                1
            
            vote
        
        
            
                3
            
            answers
        
        
            
                83
            
            views
        
        
            
            
            
        Solve Projection Problem with Linear Equality and Box Constraints
                    Solve the following problem:
$$
\begin{alignat*}{3}
\arg \min_{ \boldsymbol{x} } & \quad & \frac{1}{2} \left\| \boldsymbol{x} - \boldsymbol{y} \right\|_{2}^{2} \\
\text{subject to} & \quad ...
                
            
       
        
            
                3
            
            votes
        
        
            
                1
            
            answer
        
        
            
                59
            
            views
        
        
            
        M2M in python using SciPy.sph_harm_y
                    I’m implementing a 3D Laplace Fast Multipole Method in Python and I’m trying to validate my P2M, potential evaluation, and M2M steps against direct summation. P2M + evaluation agree with direct ...
                
            
       
        
            
                1
            
            vote
        
        
            
                0
            
            answers
        
        
            
                196
            
            views
        
        
            
            
        Numerical Simulation of Path Integral
                    I'm trying to code the path integral simulation described in this article. I followed all the instructions, but no matter how hard I try, my program won't work. Could someone please explain what I am ...
                
            
       
        
            
                0
            
            votes
        
        
            
                0
            
            answers
        
        
            
                46
            
            views
        
        
            
        What does the k and l in top-k and bottom-l usually stand for?
                    I've seen a lot of papers with the terms top-k, bottom-l, avg-k, etc e.g. in this paper. Is there a usual meaning for the k and l that I'm unaware of? Maybe I didn't pay due attention during ...
                
            
       
        
            
                4
            
            votes
        
        
            
                0
            
            answers
        
        
            
                54
            
            views
        
        
            
        What are the current best practices for reproducibility when simulations depend on chaotic systems?
                    I’m running computational fluid dynamics simulations where small perturbations lead to divergent trajectories (chaotic systems). Even with fixed seeds, compiler differences and rounding lead to ...
                
            
       
        
            
                0
            
            votes
        
        
            
                0
            
            answers
        
        
            
                58
            
            views
        
        
            
            
        More accurate ways to approximate a Finite Fourier Sine Transform?
                    I'm reading through this research paper, and from section 2, I have deduced that the Finite Fourier Sine Transform along $0\le x\le L_x$ is formulated as
$$F(m) = \frac{2}{L_x}\int_0^{L_x} f(x) \sin\...
                
            
       
        
            
                1
            
            vote
        
        
            
                4
            
            answers
        
        
            
                450
            
            views
        
        
            
            
        Relativistic correction when integrating equations of motion for charged particles in static electromagnetic fields?
                    I'm numerically integrating electron trajectories through a combination of static electric and magnetic fields.
The classical equations of motion are trivial, but since kinetic energy may reach 20 keV ...
                
            
       
        
            
                0
            
            votes
        
        
            
                0
            
            answers
        
        
            
                45
            
            views
        
        
            
        Is SUPG/VMS very easy to diverge at flow region that transitions from hypersonic flow to subsonic flow happens?
                    I am coding with FEnicsx by using variational multiscale method and add an extra shock capturing method, about a high pressure gas ball(1cm radius with pressure 65 atm and 1400K), in the initial case, ...
                
            
       
        
            
                0
            
            votes
        
        
            
                0
            
            answers
        
        
            
                28
            
            views
        
        
            
        Force oscillations during vc-relax of ZnCu₂SnS₄ with quantum espresso
                    I am running a vc-relax calculation in Quantum ESPRESSO for ZnCu₂SnS₄.
Cutoff energy: 75 Ry
k-point mesh: 6×6×6
SCF converges well (total SCF corrections ~1e-5 Ry)
Issue:
Total forces oscillate ...
                
            
       
        
            
                4
            
            votes
        
        
            
                1
            
            answer
        
        
            
                128
            
            views
        
        
            
            
        How do large-scale simulations maintain numerical stability when chaotic systems amplify rounding errors?
                    Chaotic systems like weather models are extremely sensitive to initial conditions, and computers can only store finite precision. What algorithms or statistical techniques are used to keep results ...
                
            
       
        
            
                1
            
            vote
        
        
            
                1
            
            answer
        
        
            
                80
            
            views
        
        
            
            
        Oscillations in 1D compressible Navier–Stokes simulation
                    I am currently struggling with the implementation of the compressible Navier–Stokes equations in my research work. To simplify the problem as much as possible, I reformulated it as a textbook-like ...
                
            
       
        
            
                3
            
            votes
        
        
            
                0
            
            answers
        
        
            
                81
            
            views
        
        
            
        Does the Aubin-Nitsche trick apply for time-dependent PDEs?
                    I have the occasion to be reviewing some of the "basics" of the finite element method. In particular, I am interested in several technical details related to the combination of the FEM with ...
                
            
       
        
            
                1
            
            vote
        
        
            
                1
            
            answer
        
        
            
                116
            
            views
        
        
            
            
            
        Relationship between eigenvalues of the inertia tensor and the principal moments
                    Given an inertia tensor, $I$:
$$
\mathbf{I} = \begin{bmatrix}
I_{xx} & I_{xy} & I_{xz} \\
I_{yx} & I_{yy} & I_{yz} \\
I_{zx} & I_{zy} & I_{zz}
\end{bmatrix}
$$
The principal ...
                
            
       
         
         
         
         
        