43
votes
Good examples of "two is easy, three is hard" in computational sciences
One example that appears in many areas of physics, and in particular classical mechanics and quantum physics, is the two-body problem. The two-body problem here means the task of calculating the ...
- 2,557
36
votes
Good examples of "two is easy, three is hard" in computational sciences
In one and two dimensions, all roads lead to Rome, but not in three dimensions.
Specifically, given a random walk (equally likely to move in any direction) on the integers in one or two dimensions, ...
- 2,012
30
votes
Good examples of "two is easy, three is hard" in computational sciences
A famous example is the boolean satisfiability problem (SAT). 2-SAT is not complicated to solve in polynomial time, but 3-SAT is NP-complete.
- 8,610
27
votes
Good examples of "two is easy, three is hard" in computational sciences
In social choice theory, designing an election scheme with two candidates is easy (majority rules), but designing an election scheme with three or more candidates necessarily involves making trade-...
- 371
27
votes
Accepted
Conserving Energy in Physics Simulation with imperfect Numerical Solver
There are a few ways to conserve energy during ODE integration.
Method 1: Symplectic Integration
The cheapest way that is to use a symplectic integrator. A symplectic integrator solves the ODE on a ...
- 11.6k
26
votes
Good examples of "two is easy, three is hard" in computational sciences
Here's one close to the hearts of the contributors at SciComp.SE:
The Navier–Stokes existence and smoothness problem
The three-dimensional version is of course a famous open problem and the subject ...
- 2,405
24
votes
What language should I use when teaching an undergraduate course in computer programming?
In 2014, I would've said Python. In 2017, I wholeheartedly believe that the language to teach undergraduates is Julia.
Teaching is always about a tradeoff. On one hand, you want to choose something ...
- 11.6k
15
votes
Accepted
Why should I renormalize physical variables?
This question is very closely related to (and possibly a duplicate of) Is variable scaling essential when solving some PDE problems numerically?.
There are still good practical reasons to ...
- 29.8k
15
votes
Accepted
Why is leapfrog integration symplectic and RK4 not, if the latter is more accurate?
TL;DR: It depends on what kind of accuracy you need.
Energy conservation does not automatically equal accuracy. Suppose, you want to simulate the solar system, and you are using a solver that – to ...
- 1,839
15
votes
Good examples of "two is easy, three is hard" in computational sciences
Simultaneous diagonalization of two matrices $A_1$ and $A_2$:
$$
U_1^T A_1 V = \Sigma_1,\quad U_2^TA_2V=\Sigma_2
$$
is covered by existing generalized singular value decomposition.
However, when the ...
- 8,382
13
votes
Accepted
What are the differences between CFD simulations and realistic ocean/atmosphere model simulations?
Atmosphere and ocean have highly-stratified flows in which the Coriolis force is a major source of dynamics. Maintaining geostrophic balance is extremely important and many numerical schemes are ...
- 25.3k
13
votes
Accepted
Galerkin method: Test functions vs. Basis functions
Suppose that the solution $u$ of the PDE lives in some function space $X$. We'll write the PDE as a bilinear form
$A(u, v) = f(v)$
for all $v$ in $X$, where $f$ is some element of the dual space $X^*...
- 8,142
12
votes
Why is it not computationally possible to accurately predict the weather that would occur after 14 days?
Further to Chris' answer:
Yes, weather (or the equations describing it) is extremely sensitive to the initial conditions. The fact that the weather system contains phenomena at pretty much all time ...
- 2,106
12
votes
Finite-difference software for solving custom equations
I'm going to assume since you mention electrodynamics that you're interested in PDEs.
You've already mentioned FEniCS in your comment. FEniCS offers a domain-specific language (DSL) called UFL.
This ...
- 8,142
10
votes
Good examples of "two is easy, three is hard" in computational sciences
There are plenty of examples in quantum computing, although I've been out of this for a while and so don't remember many. One major one is that bipartite entanglement (entanglement between two ...
- 209
9
votes
Accepted
Computer Build for Scientific Computing
Many of us in scientific computing simply have well-equipped laptops for regular software development tasks, some multicore workstations for smaller-scale testing, and access to clusters for larger ...
- 50.7k
9
votes
Why is it not computationally possible to accurately predict the weather that would occur after 14 days?
Errors grow exponentially in a chaotic system, and most people believe weather is chaotic. So even if you get a fairly exact numerical approximation, the fact that your input data (temperatures, etc.) ...
- 11.6k
9
votes
Accepted
Consumer hardware for scientific computing?
One issue that you should be aware of is that NVIDIA has a market segmentation strategy in which it sells relatively inexpensive GPU's to the gaming and graphics workstation markets (GeForce and ...
- 17.6k
8
votes
In Matlab, how can I be consistent with units?
Just simply by being consistent in all of my code?
Yes this is the only way. Matlab or any other programming language does not know about units. They only know about numbers.
As an example consider ...
- 2,983
8
votes
Accepted
In Matlab, how can I be consistent with units?
I would say that you have, mainly, two methods:
Being consistent in all your code, as already suggested in another answer. For that purpose, I always keep a table like this one with me, since it ...
- 8,006
8
votes
Why there are people that still prefer fortran 77 over new versions?
I think it's generally true that there are no advantages of Fortran 77 over either newer versions of Fortran or in fact any number of other programming languages that are widely used in scientific ...
- 50.7k
8
votes
Efficiently finding all (x,y,z) points within certain distance of point P
You can use Morton keying to sort the coordinate locations by binning them into cubes of some specified size $d$. This is an $\mathcal{O}(N\log N)$ operation. Then, given any point P, you can use its ...
- 974
8
votes
Finite-difference software for solving custom equations
You might want to check out DifferentialEquations.jl. It supports ODEs, PDEs, stochastic equations, delay equations, and basically everything else. It also has really good automatic sparsity detection ...
- 211
8
votes
Accepted
Why are fluid simulations so hard?
At their core, simulation of realistic fluid behavior is difficult because of the length scales involved. If you have turbulent flow, as is common in hydrodynamic or magneto-hydrodynamic situations, ...
- 50.7k
7
votes
How can an engineering student become a computational scinece expert in a short time
There is no shortcut. Just like there is no shortcut to becoming an "engineering expert in a short time".
The thing is that to be an expert in civil engineering, you need to understand load analysis, ...
- 50.7k
7
votes
Calculating partial trace of array in NumPy
I didn't follow exactly your notation so I can't say for sure what this would look like for your example, but two suggestions:
if you really want your life made simple, check out ...
- 191
7
votes
Good examples of "two is easy, three is hard" in computational sciences
Angle bisection with straightedge and compass is simple, angle trisection is in general impossible.
- 2,557
7
votes
Accepted
Why the magnetisation shows abrupt behaviour for this 3D ising spin system
Your lattice consists of 5 x 5 x 5 = 125 spins, so your number of Montecarlo steps to reach equilibrium should be >> 125, because you randomly picking a site and flipping it, so random numbers should ...
- 276
6
votes
Accepted
Solve implicit ODE numerically in orbit simulation
(1) Using the previous value of $\ddot{r}_j$ is like adding an error term to the r.h.s. of your equation of magnitude $\mathit{const}\times(\ddot{r}_j(t+\delta t) - \ddot{r}_j(t))$, meaning your ...
- 11.4k
Only top scored, non community-wiki answers of a minimum length are eligible