Most Popular

1500 questions
7
votes
2 answers

How to compute the wavelet approximation of a function?

For the function $f(x)=x$, how to compute the wavelet approximation using Haar basis? I'm new to wavelet, I'm looking for a package which will do something like this from mpmath import * mp.dps = 15; mp.pretty = True nb_coeff = 3 interval = [0,…
user729
7
votes
4 answers

Introductory book on computational physics

I'm currently working on my MS in CS and have developed an interest in astrophysics. Luckily one of my professors is a astrophysicist and is currently doing research through computational physics and he agreed to help me get started. I'm currently…
inzombiak
  • 115
  • 4
7
votes
5 answers

Open source codes for 2D instationary Navier Stokes equations

What open source tool can be recommended for solving 2D instationary Navier Stokes equations (in simple geometries, but with high Reynolds numbers)? Most packages I found, I'm not very lucky with. OpenFOAM seems to be very powerful, but I think it…
Thomas W.
  • 810
  • 5
  • 16
7
votes
2 answers

Feasibility and software for acoustic simulation

I'm looking at doing a finite-element simulation of air flow essentially for the purposes of approximating the response to an external audio impulse of a smallish (~10-30 cm scale), stationary 3D-printed structure containing a microphone. I'm mostly…
Sami Liedes
  • 171
  • 3
7
votes
1 answer

Exploiting parallelism within a single heterogeneous computing node

I'm looking for a library to exploit parallelism within a single heterogeneous computing node (possibly using Accelerators like Xeon Phi or nVidia's GPGPU's) in a C++ FV/DG code using hierarchical octree-like grids. It should support multiple…
gnzlbg
  • 1,075
  • 10
  • 18
7
votes
2 answers

What heuristics can be used to minimize the asymptotic matrix bandwidth of a 5-point Laplacian discretization?

I can see that there are multiple heuristics to achieve a matrix with minimum bandwidth. As heuristics, they can't guarantee an optimal solution in polynomial time (after all, the problem is NP-complete). I wonder, then, for the finite difference…
Paul
  • 12,045
  • 7
  • 56
  • 129
7
votes
1 answer

Generation of variable with given auto-correlation function

How can I generate realizations of random complex variable $x(t)$ with a given autocorrelation function $C(s)$, defined by $$C(s) = \langle x(s) x(0) \rangle$$ and obeying the condition $C(-s) = C^*(s)$? Any link related to similar random variable…
7
votes
4 answers

Need Fortran 77 compiler

Does anyone know a compiler for Fortran 77 available as a free download? I have pre-written 77 code from a source published in the early 90's that I just need to compile, build, and run. But I don't know enough about Fortran to update it myself.
user5246
  • 71
  • 1
  • 2
7
votes
1 answer

Numerically solving systems of about 100 ODEs

I am looking to solve large systems of non-linear ODEs. There appears to be a very large list of methods available varying in complexity, and I have a hard time searching through them and picking one. Are any of these methods preferred for large…
RobVerheyen
  • 529
  • 1
  • 3
  • 11
7
votes
5 answers

Computational science contests. Why arent there any?

I was wondering why there are no online or offline computational science contests? At least I couldn't find much by googling. I mean, like a topcoder for computational sciences. I assume one reason is providing massive computational power for few…
7
votes
3 answers

A problem in 1D linear finite element method

When applying Galerkin method, we have two conventions, i.e. multiply the test function $v$ at left/right, $(v,u)/(u,v)$. Both ways won't matter for a simple problem like Poisson's equation, since the stiffness matrix $K$ and Mass matrix $M$ are…
lorniper
  • 593
  • 4
  • 13
7
votes
2 answers

Examples of numerical solution of stochastic differential equation(SDE)?

I want to simulate a nonlinear stochastic differential equation $$ {\rm d}X_t = f(X_t) {\rm d}t + g(X_t){\rm d}B_t $$ where $f,g \in C^{\infty}({\mathbb R}^n ,{\mathbb R})$ and $B_t$ is one-dimensional standard brownian motion. How can I do it by…
7
votes
3 answers

Conforming mesh refinement for quads/hex elements

The context - I'm working with a spectral FE (higher order interpolation at GLL nodes) code on conforming hexahedral meshes, and our PI is interested in improving mesh quality, possibly with adaptive refinement. However, the only local refinement…
Jesse Chan
  • 3,132
  • 1
  • 13
  • 17
7
votes
3 answers

Finite-size effects in molecular dynamics simulation

I will be blunt: this is an old exam question of a course that covers Molecular Dynamics (MD) and Monte Carlo (MC) simulations. Explain how to determine the severity of finite-size effects. My idea: Since we are talking about MD and MC,…
seb
  • 994
  • 1
  • 9
  • 22
7
votes
2 answers

Convergence/stagnation of BiCGStab(l)

I am solving 3D time-harmonic Maxwell FDFD problems (which result in huge sparse linear systems) using BiCGStab(l). I have tried out a bunch of different methods and for my specific use case, it seems like this algorithm outperforms the other common…
Costis
  • 1,330
  • 11
  • 18