1
$\begingroup$

Do you know any repository of huge matrices of a linear system?, Or to tell a problem in which a lot of linear equations are needed?

I want to solve huge linear equations systems but have not had any luck when trying to find huge matrices.

  • 0
    I'm not sure I understand your request. Could you maybe indicate in what context you want to do this? is the idea to test your own solver and benchmark it against others? if so why not generate a random matrix and independent vector and solve in the mean-square sense (say) Ax=b? Also, do you want to consider *any* kind of system or structured systems?2012-04-23
  • 1
    Why not generate a huge matrix, generate a solution, construct the rhs, and then solve?2012-04-23
  • 0
    Ok, I understand, could you just tell a problem in which huge linear systems are used, like in economy or something...2012-04-23
  • 0
    Here is a nice one. [Interpolation using Vandermone matrix](http://en.wikipedia.org/wiki/Vandermonde_matrix#Applications). Generate a random polynomial of degree 10000000 or any bgillion degree you like. Warning though: numerically ill for very large degrees..2012-04-24
  • 1
    @J.D. "numerically ill" - indeed, but the growth rate very much depends on your choice of interpolation points. If the interpolation points are the roots of unity, the Vandermonde matrix is, up to scaling, unitary (and thus, well-conditioned)...2012-04-24
  • 0
    In numerical methods for solving partial differential equations - eg finite elements, finite differences, etc - you end out with arbitrarily large matrices as you make the mesh finer and finer.2012-04-24

3 Answers 3

4

You might try http://math.nist.gov/MatrixMarket/

  • 0
    ...and the [links to other sites](http://math.nist.gov/MatrixMarket/links.html) show other places to get nice examples!2012-04-23
2

Let $n\geq 1$ be an integer, and fix $n+1$ points in $\mathbb{R}^2$, say $(x_1,y_1)$, ... , $(x_n,y_n)$, $(x_{n+1},y_{n+1})$, such that $x_1<x_2<\ldots<x_n$.

Problem: find the coefficients of the unique polynomial $p(x)\in \mathbb{R}[x]$ of degree $n$ that interpolates all $n+1$ points, i.e., find coefficients $a_0,\ldots,a_n\in\mathbb{R}$ such that the polynomial $$p(x)=a_0+a_1x+\cdots+a_nx^n$$ satisfies $p(x_k)=y_k$ for all $k=1,\ldots,n+1$.

  • 0
    Oops, I just noticed that J.D. pointed this one out in a comment.2012-04-24
1

Vast systems of linear equations are solved as part of some integer factorization algorithms, such as the Quadratic Sieve, and the Number Field Sieves (both special and general). Much literature can be found by searching on those terms.