Questions tagged [numerical]

This tag is for questions concerning problems using numbers which either cannot be exactly solved, or where the exact solution may be much more difficult to acquire than by using numerical methods.

This tag is for questions concerning problems using numbers which either cannot be exactly solved, or where the exact solution may be much more difficult to acquire than by using numerical methods.

It includes questions about algorithms, precision and accuracy of underlying numbers, advantages and disadvantages of several methods, error propagation, problems with the underlying numbers and data structures.

508 questions
4
votes
1 answer

Numeric Keyboard?

I have several instances in my application that I would prefer to use a numerical keyboard. Does one exist in WP7 and if so, how would I go about popping that one up for a text box instead of the normal one?
Hosemeyer
  • 1,264
  • 4
  • 17
  • 29
4
votes
1 answer

How do the Sho dll's from Microsoft Research compare to the open-source Math.NET numerics project

I'm considering Sho from the perspective of developing general .NET applications and something that has been difficult for me in the past is that there is NO standard math library for the .NET platform. Aside from the commercial implementations, two…
4
votes
5 answers

LU Decomposition from Numerical Recipes not working; what am I doing wrong?

I've literally copied and pasted from the supplied source code for Numerical Recipes for C for in-place LU Matrix Decomposition, problem is its not working. I'm sure I'm doing something stupid but would appreciate anyone being able to point me in…
Bolster
  • 7,460
  • 13
  • 61
  • 96
4
votes
2 answers

how does numpy.linalg.eigh vs numpy.linalg.svd?

problem description For a square matrix, one can obtain the SVD X= USV' decomposition, by using simply numpy.linalg.svd u,s,vh = numpy.linalg.svd(X) routine or numpy.linalg.eigh, to compute the eig decomposition on Hermitian matrix X'X and…
4
votes
2 answers

Java / Scala math library with elliptic integrals and bessel functions?

I am looking for a math library for scientific computing to use in Java / Scala. Especially I need complete elliptic integrals und modified Bessel functions. I would be nice if it is open source, but I guess I will have to take whatever is out…
Radfahrer
  • 41
  • 1
  • 2
4
votes
1 answer

UMFPACK and BOOST's uBLAS Sparse Matrix

I am using Boost's uBLAS in a numerical code and have a 'heavy' solver in place: http://www.crystalclearsoftware.com/cgi-bin/boost_wiki/wiki.pl?LU_Matrix_Inversion The code works excellently, however, it is painfully slow. After some research, I…
ccook
  • 5,869
  • 6
  • 56
  • 81
4
votes
5 answers

Code to solve determinant using Python without using scipy.linalg.det

Description(this is a hwk question): I am not sure where to start here. I plan to use Laplace's Expansion but I am not sure how to implement it for nxn matrices. Any help would be appreciated. Note: I already have a function to generate random…
theAngryPhysicist
  • 129
  • 1
  • 5
  • 12
4
votes
3 answers

Numerical problems in Matlab: Same input, same code -> different output?

I am experiencing problems when I compare results from different runs of my Matlab software with the same input. To narrow the problem, I did the following: save all relevant variables using Matlab's save() method call a method which calculates…
Philipp
  • 11,549
  • 8
  • 66
  • 126
4
votes
2 answers

n-dimensional integration in R with limits as functions

I'm aware of R packages such as cubature which are able to perform numerical integration for many integrals. However, in my situation the limits of integration are not integers, but functions of the other variables. For instance, suppose I have…
nathanesau
  • 1,681
  • 16
  • 27
4
votes
1 answer

Bizarre numerical errors

I am implementing a simple two-state particle filter. If you don't know what a particle filter is, that's fine - the short version is that I need to compute weighted means with weights between 0 and 1, and values between 0 and 1. Each particle has…
user650261
  • 2,115
  • 5
  • 24
  • 47
4
votes
3 answers

Numerical bug in intersecting the equation of ray and torus when the camera is far from the torus

I am trying to ray trace a torus without triangulating the torus and just by intersecting the ray and torus analytic equation. I did that with the following code: void circularTorusIntersectFunc(const CircularTorus* circularToruses, RTCRay& ray,…
mmostajab
  • 1,937
  • 1
  • 16
  • 37
4
votes
5 answers

Numerical Pattern Matching

A project I'm researching requires some numerical pattern matching. My searches haven't turned up many relevant hits since most results tend to be around text pattern matching. The idea is we'll have certain wave patterns we'll need to be watching…
Timothy Strimple
  • 22,920
  • 6
  • 69
  • 76
4
votes
1 answer

Error at lapack cgesv when matrix is not singular

This is my first post. I usually ask classmates for help, but they have a lot of work now and I'm too desperate to figure this out on my own :). I am working on a project for school and I have come to a point where I need to solve a system of linear…
Jan Malec
  • 96
  • 4
4
votes
3 answers

Numerical derivative of a vector

I have a problem with numerical derivative of a vector that is x: Nx1 with respect to another vector t (time) that is the same size of x. I do the following (x is chosen to be sine function as an example): t=t0:ts:tf; x=sin(t); xd=diff(x)/ts; but…
milad
  • 169
  • 1
  • 2
  • 10
4
votes
1 answer

convolution of positively supported functions in R

I want the convolution of two functions defined on [0,Inf), say f=function(x) (1+0.5*cos(2*pi*x))*(x>=0) and g=function(x) exp(-2*x)*(x>0) Using the integrate function of R I can do this, cfg=function(x) integrate(function(y)…