[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: sub_pc lu zero pivot
- To: petsc-users@xxxxxxxxxxx
- Subject: Re: sub_pc lu zero pivot
- From: "Matthew Knepley" <knepley@xxxxxxxxx>
- Date: Fri, 8 Feb 2008 12:02:45 -0600
- Dkim-signature: v=1; a=rsa-sha256; c=relaxed/relaxed; d=gmail.com; s=gamma; h=domainkey-signature:received:received:message-id:date:from:to:subject:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; bh=6radnUs6293L7ZFPwFas3qhZizZwezJ/P7uEROXtjWs=; b=aE+ll/c8f9y4rUqMwzKQMBuhzV6lAy1ecjEiCrgdjx32gKq84R71UzFGUJTJsKByrG/ZdRti58mwqv+HDI5XZyR2W6rKS8VQH+zDwKz1yCIHQuMOuvKFDbbcfgJy9NAFj74AZOULqRrN5leptM+HNXdP3hRGBypF1VeyWr1lWWQ=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=gamma; h=message-id:date:from:to:subject:in-reply-to:mime-version:content-type:content-transfer-encoding:content-disposition:references; b=u7fqETHWuw5C0YeO49s3x1SquZkMhDgjtRkeW5AJg4mXYK+0QMNgAsDeMw8YWuWVBr9wZlA0b4tI4SKcgC5nXy3He11y8FaCyBbbhLnXpB/969IjE2SMelnEluOoUOzEdDXQQCqQZpzLaTvUXDROKf7yv9wsP1aSWTChTt0xhsY=
- In-reply-to: <20080208104609.AAW99320@batman.int.colorado.edu>
- References: <20080208104609.AAW99320@batman.int.colorado.edu>
- Reply-to: petsc-users@xxxxxxxxxxx
- Sender: owner-petsc-users@xxxxxxxxxxx
We would like to see:
1) The complete error message
2) The output of -snes_view so we can see exactly what solver setup you have
Thanks,
Matt
On Feb 8, 2008 11:46 AM, Andrew T Barker <Andrew.Barker@xxxxxxxxxxxx> wrote:
>
> For my code on a single processor, using the options
>
> -ksp_type gmres -pc_type lu
>
> works fine, but
>
> -ksp_type gmres -pc_type asm -sub_pc_type lu
>
> produces a zero pivot. On one processor I would expect these two to be identical, and in fact when I print out the asm submatrices (there is only one) in the second case and the original matrix in the first case, they are identical. Even stranger is that
>
> -ksp_type gmres -pc_type lu -sub_pc_type lu
>
> also produces a zero pivot. I would expect it to ignore the sub_pc in this case.
>
> Thanks for any help,
>
> Andrew
>
> ---
> Andrew T. Barker
> andrew.barker@xxxxxxxxxxxx
> Applied Math Department
> University of Colorado, Boulder
> 526 UCB, Boulder, CO 80309-0526
>
>
--
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