# All Questions

7,839 questions
Filter by
Sorted by
Tagged with
99 views

### Inverting a matrix from LU decomposition

The LAPACK routines xGETRI compute the inverse of a matrix $A = PLU$ in its LU decomposed form by first computing $U^{-1}$, and then solving the system: $$(A^{-1} P) L = U^{-1}$$ My question is: ...
125 views

### Fast and free server for computing

I have to calculate a huge differential equation. With my laptop, it's going to be computed for several days. Is there a free (I need just for 3 days) fast server for scientific calculations? My ...
80 views

35 views

### Computing excited states using itensor (with DMRG)

I am trying to compute first few excited states of some Hamiltonian (I am using itensor and its DMRG algorithm). To do so, I am ...
22 views

### Reduce projection error while retaining similar amount of elements in CG-FEM

Based on the answers I got to my questions (Interpolation of function onto mesh gives different results, depending on mesh density and Solving a non-linear heat equation with the galerkin method gives ...
26 views

### How can I implement the Invaded Cluster Algorithm for a network of Ising spins?

My primary concern is about finding the percolating cluster for any given network. For a lattice it is straight forward : When the size of a cluster reaches the length of the lattice, then it is said ...
27 views

### Numerically computing deflection due to thermal expansion

Using linear elasticity formulation, I am attempting to numerically compute the displacement due to thermal expansion. This is done for a 3-D isotropic material. The governing equations are simply: ...
2k views

### Mathematically, why does mass matrix / load vector lumping work?

I know that people often replace consistent mass matrices with lumped diagonal matrices. In the past, I've also implemented a code where the load vector is assembled in a lumped fashion rather than ...
65 views

### Re-using LU factorization within iterative (?) setup for a sum of two matrices

So, I would love to make at least some use of my preexisting data, no matter how small, and just out of ideas. Maybe I am just a prisoner of a Kahneman-like theatre-ticket paradox, and don't know ...
41 views

### Symplectic linear multistep method?

I'm doing a gravitational n-body simulator and I'm thinking of implementing linear multistep methods like Adam-Bashforth. But is there any symplectic multistep methods?
43 views

### Isotropic thermal expansion

I frequently see the equation $$\sigma_t = E\alpha \Delta T$$ as the equation for thermal stress. Where $E$ is Young's modulus, $\alpha$ is the CTE, and $\Delta T$ is the change in temperature. ...
83 views

### What method of Finite difference is this?

I am reviewing Numerical Recipes method on solving ODEs via relaxation (Chapter 18.3 in the 3rd edition) and they chose a finite difference method I am unfamiliar with (Equation 18.3.2): \begin{...
212 views

26 views

### Obtaining integer digits using the GNU Multiprecision Arithmetic Library (gmplib)

I'm using the GNU Multiprecision Arithmetic Library (gmplib) for some experiments in computational mathematics. I want to extract, and manipulate, the base-b digits (with 2 <= b <= 10) of ...
61 views

### How to check experimental data against a theoretical curve? (Python)

I am trying to check the agreement of a dataset against a theoretical curve, specifically a bandstop filter in an RLC circuit. I have generated a function which describes the curve we expect from the ...
138 views

### Automatic timestep adjustment in a CFD solver

I have developed my own 3D Finite Volume Navier-Stokes solver based on projection method for nonuniform grid. I am looking to incorporate automatic timestep adjustment at each time step based on ...
I write a test program to integrate foward on $[0,T_f]$ and then backward on $[T_f,0]$ from the endpoint of the forward integration an Hamiltonian system:  \dot q(t) = \frac{\partial H}{\partial p}(...
Consider the following space $A = \{(x_1,x_2,x_3)\in \mathbb{R}^3|x_1+x_2+x_3 = 1\}.$ Then say that we want to minimize a function $J(y):\mathbb{R}^{3}\to \mathbb{R}$ subjected to the constraint that \$...