[Date Prev][Date Next][Thread Prev][Thread Next][Date Index][Thread Index]
MATSHELL a bit broken...
- To: petsc-dev@xxxxxxxxxxx
- Subject: MATSHELL a bit broken...
- From: "Lisandro Dalcin" <dalcinl@xxxxxxxxx>
- Date: Thu, 25 Oct 2007 18:20:48 -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=BiFzcbrjw8+XDMhiXdaWj+kWLsYRgKaY1ydUU+sf24Q=; b=bqWaehKQrC8HlTImZiRwiH9WVvi3zkuMUzoWgtv+2aCfYAvuuoFkWZi90f7YEuSMJ8ZPmdea9M8VfnwvK8QBaqx694gnfIWdwACCRmwUg1tsfAmXDledCl3E1+uJCrkrfBeUD+ZAl0dPat4AiVfcXlHRr0qRhFIJeFmmQWPESHo=
- 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=EBUmbwVLn+X2v+Ma/94V55hQi0lMKOgMsq3cnluIO80HVTBlAY9y4k8vT40rA4U3B4RNZJlScm0LNRjbXHNOyte0dU8gzQIbo1ct2iy3EUwZl9rZ2PC1GKwYoAcu8LYUOgdOyo3s4MuSJKobQo1mJVl30Bz2hZEsoxFPMW0wn2s=
- Reply-to: petsc-dev@xxxxxxxxxxx
- Sender: owner-petsc-dev@xxxxxxxxxxx
It seems MATSHELL implementation is a bit broken. If I'm not wrong, if
a user set the operations MOP_ASSEMBLY_END, then the use of the
automatic MatScale/MatShift would lead to unexpected result...
Should this operation be managed the same way as
mult/multtranspose/get diagonal currently is, that is, by backupping
the user provided function pointer?
BTW, sorry if I disturb you asking for making changes that seems
trivial, but I really do not like to change any line of code after
knowing your opinion.
--
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