[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: Order (N) Sparse Banded Linear Solver
- To: petsc-users@xxxxxxxxxxx
- Subject: Re: Order (N) Sparse Banded Linear Solver
- From: "Matthew Knepley" <knepley@xxxxxxxxx>
- Date: Thu, 4 Oct 2007 09:01:54 -0500
- Cc: "Tim Stitt" <timothy.stitt@xxxxxxxx>
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=beta; h=domainkey-signature:received:received:message-id:date:from:to:subject:cc:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; bh=wLwB1cA4Jswgh2eY6zGKy6VGk/Xlik3U/qugQk2Y5PA=; b=i4da76UenHG+EK2wbPTp9zmNtorhjqhoWB8QbTQCz1Ef+oI+epLrqhAu5Zkyu0GXk5rKLWt8R2uTeU6OXqmYGvDd2xqSIlMHYHQdqH4wE+bmbYAiJdq4hGWYkeO/S3kdkDB2eT/rHDkO18FEsapeKEG3upf9VvBa0kdJKQGVeH0=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=beta; h=received:message-id:date:from:to:subject:cc:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; b=eM+1HPuo8luTOelMUTseHUbH3RllssUDDBQLPNi0olXVgI74bVyVQQnv5pIxXKo8VrLl6ZYqxR8qcgEDZYtd6ypJSsH0+qrpuD8XZbdMoE54E4P/ugxi1AqX0WL4ZSVpQqiRLaaa9+AUiWd2RqTmRrU6pfDNOpQjZoJns2UeXsI=
- In-reply-to: <4704F0CD.1000803@ichec.ie>
- References: <4704CE87.70202@ichec.ie> <a9f269830710040515w4ebe838dubac12ed7b45b6b65@mail.gmail.com> <4704F0CD.1000803@ichec.ie>
- Reply-to: petsc-users@xxxxxxxxxxx
- Sender: owner-petsc-users@xxxxxxxxxxx
On 10/4/07, Dr. Timothy Stitt <timothy.stitt@xxxxxxxx> wrote:
> Matthew,
>
> Out of curiosity what would you term as very very large? The matrices I
> am using will probably be of order 10^3 x 10^3 elements.
I believe this is fairly small. SuperLU and MUMPS routinely report results for
matrices of order 10^5x10^5 (on large parallel systems).
Matt
> Thanks,
>
> Tim.
>
> Matthew Knepley wrote:
> > On 10/4/07, Dr. Timothy Stitt <timothy.stitt@xxxxxxxx> wrote:
> >
> >> Dear PETSc Community,
> >>
> >> I was wondering if anyone knew of any sparse/dense banded linear solvers
> >> that run in order (N)? Does PETSc provide such a solver or link to an
> >> appropriate external solver? There seems to be a lot of links online to
> >> theoretical approaches but we are having difficulties in sourcing an
> >> implementation? I am hoping PETSc might be our saviour ???
> >>
> >
> > Have you tried MUMPS or SuperLU, both available through PETSc's configure?
> > They are not truly O(N), but its hard to tell the difference until you get very
> > very large, or have a nice elliptic PDE on a regular grid. In the
> > latter case, you
> > should use the PETSc DMMG.
> >
> > Matt
> >
> >
> >> Thanks in advance,
> >>
> >> Tim.
> >>
> >> --
> >> Dr. Timothy Stitt <timothy_dot_stitt_at_ichec.ie>
> >> HPC Application Consultant - ICHEC (www.ichec.ie)
> >>
> >> Dublin Institute for Advanced Studies
> >> 5 Merrion Square - Dublin 2 - Ireland
> >>
> >> +353-1-6621333 (tel) / +353-1-6621477 (fax) / +353-874195427 (mobile)
> >>
> >>
> >>
> >
> >
> >
>
>
> --
> Dr. Timothy Stitt <timothy_dot_stitt_at_ichec.ie>
> HPC Application Consultant - ICHEC (www.ichec.ie)
>
> Dublin Institute for Advanced Studies
> 5 Merrion Square - Dublin 2 - Ireland
>
> +353-1-6621333 (tel) / +353-1-6621477 (fax) / +353-874195427 (mobile)
>
>
--
What most experimenters take for granted before they begin their
experiments is infinitely more interesting than any results to which
their experiments lead.
-- Norbert Wiener