Most Popular

1500 questions
9
votes
3 answers

GPU libraries to use for Linear Algebra operations

I am looking for GPU libraries to accelerate a code I have, whose primary "high performance" regions have singular value decomposition, QR factorization and Eigen values, Eigenvectors computation. I have looked through Nvidia's website on its tools…
atmaere
  • 213
  • 1
  • 6
9
votes
5 answers

Is Discrete Exterior Calculus currently a focusing point in numerial computing world or simulation in industry,

I am just wondering if Discrete Exterior Calculus, as a new numerical method , is good at numericall solving problems in elasticity, fluids or other physical/real area.
user5226
  • 91
  • 2
9
votes
1 answer

How to efficiently determine the intersection of a vertical cutting plane with a mesh

I have a list of vertical cut planes, and I have a polygonal mesh ( it's a 2D+0.5D mesh, something like a 2D mesh with an extra dimension, $z$ attached to each point). One can assume that the mesh contains vertices $V$ , edges $E$ and face $F$. The…
Graviton
  • 488
  • 5
  • 16
9
votes
2 answers

Peculiar error when solving the Poisson equation on a non-uniform mesh (1D only) finite volume method

I have been trying to debug this error the last few days I wondered if anybody has advice on how to proceed. I am solving the Poisson equation for a step charge distribution (a common problem in electrostatics/semiconductor physics) on a non-uniform…
boyfarrell
  • 5,409
  • 3
  • 35
  • 67
9
votes
1 answer

Matlab solution for implicit finite difference heat equation with kinetic reactions

I am trying to model heat conduction within a wood cylinder using implicit finite difference methods. The general heat equation that I'm using for cylindrical and spherical shapes is: Where p is the shape factor, p = 1 for cylinder and p = 2 for…
wigging
  • 325
  • 1
  • 4
  • 11
9
votes
2 answers

Continuity of eigenvectors of parametric matrix

I have $n$-dimensional matrices $\mathrm{\hat{H}}(\vec{k})$ depending on vector parameter $\vec{k}$. Now, eigenvalue routines return eigenvalues in no particular order (they are usually sorted), but I want to trace eigenvalues $E_i$ as smooth…
tomic
  • 91
  • 2
9
votes
3 answers

Can my project survive without Object Orientation?

I am writing a small MATLAB package which will solve a certain class numerical problems. There are 3 stages of the algorithm and the user has 5 choices for each stage. I have implemented the entire problem using $\approx 20$ functions and 3 switch…
MATLABOOP
  • 93
  • 2
9
votes
1 answer

Is there some good mailing list for `Computational Science`?

I am wondering whether there is some very good mailing list or google groups for Computational Science, where we can discuss questions instead of only asking and replying questions. In fact, I am more interested parallel computing and numerical…
eccstartup
  • 367
  • 1
  • 10
9
votes
7 answers

Is there a standard rating system for scientific journal publications?

I have heard that some journals are rated more highly than others. Is this true? And if so, what are the criteria for judging the value of one peer reviewed journal over another? How do I find out its rating? Will my publication be of less…
Paul
  • 12,045
  • 7
  • 56
  • 129
9
votes
4 answers

How to start the Simplex method from a feasible internal point?

I have one algorithm that generates a feasible solution to a linear programming problem. However, it is very likely that this is not a corner point. This makes it not suitable for direct use as an initial feasible solution for a bounded Simplex…
Dylan
  • 192
  • 1
  • 5
9
votes
1 answer

Shortley-Weller finite difference method

can you give me a link for a good and simple explanation of the Shortley-Weller finite-difference scheme? I tried to google it but all I get is (inaccessible) academic publications. I also tried reading the dedicated chapter (4.8) in the Wolfgang…
Michael
  • 91
  • 2
9
votes
1 answer

Solving linear systems by fft

I read in a paper and also at wiki that we can solve the system $$Ax=B$$ by Fast Fourier Transform, where $A$ is a circulant matrix. The solution is $$x=\mathtt{ifft}(\mathtt{fft}(B)/\mathtt{fft}(a))$$ where $a$ is first column of $A$, ifft is the…
Ömer
  • 643
  • 6
  • 10
9
votes
0 answers

Simple turbulence model appropriate for buoyancy-driven cavity like problem

Which turbulence model is suitable for resolving incompressible buoyancy-driven flow of a fluid within an cylindrical ampoule? I prefer turbulence model which is sufficiently simple so that fully coupled (UFL) variational form of…
Jan Blechta
  • 927
  • 5
  • 13
9
votes
1 answer

Example of a PDE model with nonlinear Dirichlet boundary conditions

Is there any application for PDEs with nonlinear Dirichlet boundary conditions? That is, I am looking for an example of a partial differential equation for a state $u$ posed on a domain $\Omega$ with $g(u|_{\Gamma})=0$, where $\Gamma$ is the…
Jan
  • 3,418
  • 22
  • 37
9
votes
3 answers

Solving a sparse and highly ill-conditioned system

I intend to solve Ax = b where A is complex, sparse, unsymmetric and highly ill-conditioned (condition number ~ 1E+20) square or rectangular matrix. I have been able to solve the system with ZGELSS in LAPACK accurately. But as the degrees of freedom…
user1234
  • 167
  • 2
  • 7