[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
Re: Something weird with SNES convergence reason
- To: petsc-users@xxxxxxxxxxx
- Subject: Re: Something weird with SNES convergence reason
- From: "Matthew Knepley" <knepley@xxxxxxxxx>
- Date: Thu, 15 May 2008 07:28:40 -0500
- 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=on3Du1gihVkQepOsTjBzqwnLNwZ1AAqSOySr7fnThPI=; b=tEhwswuEX9p43X+Ggp6r3BhAIrS6ehDUtXXqEVKbz6w6n7EHCkqWEZo1Ir12tNThOQkhe2TESPNs2Lp7sk0LGFCqSZIh4AEV8Mt9e1aqe8q4U6BKUprXr36Qy21UhfoBXIZEf5XekJjRJ4opVAotdAGkrSVzwislrftDEDgPXXU=
- 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=JCcz/l5+jX7dn6Up41nRHbC5gug59rfpIrOyo5FM9pskIxYKE99dg4oSq8X5Pr0tVA3hRtDInpp9P1npGh9f6Bg+Nl/meoZtA/KVyltf36iCKDZSOaFExNq6lhqt8wos+as8CVqWt7TpLy7UtDmPg9b4XgvONflenvqHhm8ph8Q=
- In-reply-to: <3b6f83d40805142030j2c3d9f58i5dd85c41da52d0db@mail.gmail.com>
- References: <3b6f83d40805142030j2c3d9f58i5dd85c41da52d0db@mail.gmail.com>
- Reply-to: petsc-users@xxxxxxxxxxx
- Sender: owner-petsc-users@xxxxxxxxxxx
1) Are the linear systems really being solved in Newton?
2) What is the Bratu parameter? Turn it off and see that you
get convergence in 1 iteration.
Matt
On Wed, May 14, 2008 at 10:30 PM, Rafael Santos Coelho
<rafaelsantoscoelho@xxxxxxxxx> wrote:
> Hello everybody,
>
> I've coded a program which solves, in parallel, the three-dimensional Bratu
> problem. Afterwards, I've run tests in a cluster to see how it would go and,
> at first, it seemed ok to me, but then I've noticed that whenever I
> increased the number of processors (from 16 to 32, for example), the program
> started to diverge due to a failure in the Line Search Newton's Method. Here
> is what a monitoring function prints out:
>
> nonlinear iteration number = 1, norm(F(x)) = 1013.53, linear iterations =
> 16
> nonlinear iteration number = 2, norm(F(x)) = 1013.33, linear iterations =
> 32
> nonlinear iteration number = 3, norm(F(x)) = 1013.33, linear iterations =
> 48
> Nonlinear solve did not converge due to DIVERGED_LS_FAILURE
>
> Indeed, one can see that the method is really diverging (for smaller tests,
> though, say N = 8 * 8 * 8, it converges).
>
> What's wrong here? Is it something with my code? If yes, how can I fix it?
>
> Best regards,
>
> Rafael
>
>
--
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