# All Questions

1,616 questions with no upvoted or accepted answers
Filter by
Sorted by
Tagged with
116 views

### The numerical solution of a (very ugly) set of integro-diferential equations

I've been having fun playing around with a simple model for reaction in a co-flow fluidized bed (I'm afraid I'm a chemical engineer, not a computer scientist). Unfortunately, simple physical ...
106 views

### ADMM for Linear Program over graph

I want to use ADMM to solve a LP defined over a graph. According to Distributed optimization and statistical learning via the alternating direction method of multipliers S. Boyd, N. Parikh, E. ...
163 views

### Poisson equation in frequency domain

I need some help in numerically solving the nonlinear Poisson's equation for electrons in frequency domain. The steady-state equation is: \nabla.(\epsilon\nabla\varphi) = q\left(n_i\...
400 views

### Numerically calculating the divergence of a set of oriented points

Say I have a set of oriented points at locations $\vec{v_i}$ with each some direction $\vec{n_i}$, in practice they represent the normals of some surface that has been non uniformly sampled. How would ...
86 views

### Can I make this numerical integration continuously differentiable?

Suppose I have the discrete values $f(x_i)$ for every discrete value $x_i$ greater than some $\varepsilon$, and I want to numerically calculate the following integral: n = \int_\...
149 views

### Eigenvalue-style optimization with quadratic constraints

Suppose $A\in\mathbb{R}^{n\times n}$ is symmetric and positive definite and that we have several symmetric matrices $B_i\in\mathbb{R}^{n\times n}$ that are low-rank and indefinite. I need an ...
238 views

### Boundary conditions in the Finite Element Method at only one side of computational domain

I want to solve a Sturm-Liouville problem in 1D, i.e., \begin{align} [p(x)\ u'(x)]'+q(x)u(x) = f(x) \end{align} with boundary conditions \begin{align} u(0)=a \hspace{1cm} u'(0)=b \end{align} How do I ...
251 views

### What are the career options for computational scientists other than being in software Industry?

I would like to know the career options for a computational scientist other than getting into software industry. I wish to use the my skills of Numerical Linear Algebra,Parallel Programming, ...
119 views

### What causes periodic humps in residual plots?

When using many iterative methods, whether for solving linear systems, looking for steady-state convergence in CFD, etc., the semilog plot of the residual often shows "humps" as the residual decays. ...
70 views

### Is there a term for Goodhart's Law in the context of optimization?

Let's say I'm optimizing something. To pick an arbitrary example, let's say I'm choosing the shape of some part to maximize strength-to-weight ratio. So I get some FEM software, parametrize the shape, ...
79 views

### Graph optimization for parallel processing

Consider the following example structure of overlapping images marked A,B,C,D. The possible overlaps are marked by gray color: The structure can be represented by a weighted undirected graph (images ...
110 views

### Computing entries of a matrix in parallel using GPU

I am trying to solve $Ax=b$ a system of linear equation. The matrix $A=(a_{ij})$ is constructed as $$a_{ij}=\int_0^1\int_0^1 K(x,y)h_i(x)h_j(y) dxdy.$$ The evaluation of integral takes a lot of time, ...
449 views

### Large meshing with tetgen

So I have a point cloud that I am creating a 3D flat rectangular surface from. I'm then turning it into a hollow box and connecting the corners by just dropping this surface mesh down. I need it to ...
645 views

### Neural network performs worse when using more input variables

This question is based more on the theory of neural networks than my particular implementation. Therefore I will leave out my code unless requested. I'm working on a project in C# which can create ...
842 views

198 views

Consider the integral $$2\pi T = -\frac{1}{2} \int_0^{2\pi} A\frac{\partial B}{\partial \xi} \ d\xi$$ where $$A= \sum_{-N}^N i \ sign(n) \ B_n e^{-in\xi}; \quad \quad B= \sum_{-N}^N \ B_n e^{-in\xi}... 0answers 75 views ### Large residual when integrating 2nd order ode close to singularity with SciPy ode / ODEPACK I am trying to integrate a 2nd order ODE with a singularity at close to the initial condition. Why do I get large residuals when I plug-in the result of my integration back into the ODE? The equation ... 0answers 238 views ### Is it possible to predict the null space of a structure from contributing elements null spaces? I am trying to solve an almost incompressible problem with heterogeneous properties by domain decomposition. Solution with CG converges slowly or divergerces completely. My problem becomes ill-... 0answers 142 views ### How big a matrix can we row reduce in reasonable time? I have very large matrices that I would like to row reduce (I need to keep track of the steps and find a basis of the kernel/image, not just find the rank). The good news is that I work mod 2 and the ... 0answers 301 views ### Frozen coefficient method (von Neumann stability analysis) Earlier it was considered that frozen coefficients method for Neumann stability analysis for finite difference scheme is more heuristic than rigorous. But I have read some information in a book by ... 0answers 991 views ### Poisson equation with pure Neumann boundary conditions (using FEM) I am trying to derive the correct variational form for the Poisson equation with pure Neumann boundary conditions, and an additional contraint \int_{\Omega} u \, {\rm d} x = 0, as described in this ... 0answers 1k views ### Python - calculation time derivative and laplacien by finite differences I would like to determine a temporal derivative and a Laplacian by the finite differences method to solve the heat equation in a 1-dimensional case. My aim is to get the sources term that is why I ... 0answers 98 views ### Tracking for two meshes I'm dealing with physical simulation (position based dynamic). Now I'm trying to realize tracking for two meshes. In order to explain what does it mean, let's assume that we have two similar meshes: ... 0answers 134 views ### Methods for integrating black box functions on a non-uniform grid If i have some function expressed as points on a non-uniform grid (I'm specifically interested in logarithmic grids, but general results are also interesting), and I want to integrate it, I believe ... 0answers 426 views ### Finite volume method I have question connected with finite volume method. Consider equation$$\frac{\partial u}{\partial t}=\operatorname{div}A\nabla u +f . \quad x\in \overline{B}_{1} (0)\subset \mathbb{R}^3 -\text{unit ...
247 views

Quadrature-free DG method using nodal Lagrangian basis are computationally very efficient. I have seen many papers using this method for linear PDE but almost no literature for non-linear PDE like ...
85 views

### How do I perform chebyshev interpolation from a to b with custom angle range?

Typically Chebyshev interpolation from $-1$ to $1$ with angle from $0$ to $\pi$: $\xi_j=\cos \left ({\pi j \over N}\right )$ $x_j=(1+\xi_j) * {L \over 2}$ $w$: $w_0=\pi/(2N)$ $w_{1,...,N-1}=\pi/(N)$...
82 views

### Computing linear combinations of sines and cosines (phasors)

I have a finite series that looks like this: $f(t) = \sum^n_{i=0} A_i cos(\Theta_i + \omega_i t) + B_i sin(\Theta_i + \omega_i t)$ That is, a finite series of pairs of phasors. What's the state of ...
229 views

### Sound Waves Simulation in 3D Environment

I want to do a simulation of sound waves including wave propagation, absorption, and reflection in 3D space. I did some research and I found this question in stackoverflow but it talks about ...
577 views

112 views

### What are the governing equations solved in coupled atmosphere-ocean models?

In (hydrostatic) atmospheric general circulation models, for example the so-calle Primitive Equations, consisting of the horizontal momentum equation, the hydrostatic balance, the continuity equation, ...
184 views

218 views

174 views

### Find maximum distance between elements given constraints on some

I have a list of numbered elements 1 to N that fit into positions on a number line starting with 1. I also have constraints for these elements: The element 1 is in position 1, and element N must be ...
104 views

### Existing software/scripts for spiral graphs?

I am looking for existing software or scripts to generate spiral graphs from cyclical (time) data, as presented by Webber and Muller. The graphs shown in the paper look like a great means of ...
607 views

### Perron-Frobenius theorem on general real symmetric matrices

From the Perron-Frobenius theorem, it might be concluded that the spectral radius is the largest eigenvalue for positive matrices, ie, matrices with strictly positive entries. In other words, the ...