[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
adding SNESSetLinearSolve()
- To: petsc-dev@xxxxxxxxxxx
- Subject: adding SNESSetLinearSolve()
- From: "Lisandro Dalcin" <dalcinl@xxxxxxxxx>
- Date: Mon, 29 Oct 2007 14:04:03 -0300
- 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:mime-version:content-type:content-transfer-encoding:content-disposition; bh=AtF3+0zkFHfDgjbkn4fBsFgQA2fBOi6Cz/secYuE6pQ=; b=ZpEIBrsSQU0lAcqPic1xmVWb2aZPus01mwhkFSrHBttPAKdpDoF+iRQXDgm7YwjJpYdcj6HyhSQbowyZJoVvYCrOIFVyH9pvGNWVLWWYSYGhTo6n3qgjFCaup+eycvmHIJw0tVR/FMAbjw9n2sdDXsjYNYDwi8L++oUq01qd5DE=
- Domainkey-signature: a=rsa-sha1; c=nofws; d=gmail.com; s=beta; h=received:message-id:date:from:to:subject:mime-version:content-type:content-transfer-encoding:content-disposition; b=ifJKIHUXPlJsIv/d8d7Y0aLpCJH+PAjJQ+DyF6fNRRnN0QOlfPpjkn+/qgHYU/XnKXlf8TsJKY4Bprqum5LNGnS1eVL/bHsXboz1H9viIDOSzy8mn3UxIc5z589LqGjk/7l/QNH08XzlE1++s0ZkBuU0jtIClxxihjx0TsmFVBM=
- Reply-to: petsc-dev@xxxxxxxxxxx
- Sender: owner-petsc-dev@xxxxxxxxxxx
Some time ago I made a request for adding SNESSetLinearSolve(), in
order to let users set a custom linear solve routine to SNES. I was
not clear to me if this idea was finally accepted, so I want to ask
again.
I really needs this feature in the near future, I I would like this to
go in the next petsc release. My use case is rather simple: I'm
helping a coworker to use SNES for solving a nonlinear optimization
problem related to mesh movement (that is, move nodes retaining
element connectivities). This problem uses a regularized functional,
wich introduces a global degree of freedom coupled with all others
dofs. To avoid assembling a sparse parallel AIJ matrix with a dense
column and row, we endup needing to solve two linear systems with the
same matrix for actually solving the Hessian of the regularized
functional.
--
Lisandro Dalcín
---------------
Centro Internacional de Métodos Computacionales en Ingeniería (CIMEC)
Instituto de Desarrollo Tecnológico para la Industria Química (INTEC)
Consejo Nacional de Investigaciones Científicas y Técnicas (CONICET)
PTLC - Güemes 3450, (3000) Santa Fe, Argentina
Tel/Fax: +54-(0)342-451.1594