close Warning: Can't synchronize with repository "(default)" (/var/svn/tolp does not appear to be a Subversion repository.). Look in the Trac log for more information.

Opened 18 years ago

Last modified 14 years ago

#513 assigned enhancement

New function about LDLt transformation

Reported by: Víctor de Buen Remiro Owned by: Jorge
Priority: normal Milestone:
Component: Math Version: head
Severity: minor Keywords:
Cc:

Description

Propongo la creación de una nueva función TOL built-in LAPACK.DSYTRF que sirva de interfaz a la función DSYTRF de LAPACK cuya descripción se puede ver en

​http://javadocs.approximatrix.com/jlapack/org.netlib.lapack.DSYTRF.html

Esta función podría usarse cuando Choleski presenta problemas numéricos. La API evidentemente sería distinta pues algunos argumentos son internos y no le interesan al usuario

Set LAPACK.DSYTRF(Matrix A, Text UPLO)

Si UPLO=='L' devuelve la descomposición triangular inferior L,D,INFO?
Si UPLO=='L' devuelve la descomposición triangular superior U,D,INFO?

Si INFO==0 entonces A es definida y D es diagonal.

a.) Si D[1,1] < 0 entonces A es definida negativa
b.) Si D[1,1] > 0 entonces A es definida positiva

b.1.) Si UPLO=='L' la descomposición triangular inferior de Cholesky es L*D(1/2)
b.2.) Si UPLO=='U' la descomposición triangular superior de Cholesky es U*D(1/2)

Change History (3)

comment:1 Changed 18 years ago by Jorge

Status: new → assigned

will work on it soon

comment:2 Changed 16 years ago by Víctor de Buen Remiro

Component: Various → Math

comment:3 Changed 14 years ago by Pedro Gea

Priority: highest → normal
Note: See TracTickets for help on using tickets.