inversion {salad}

Determinant and matrix inversion for dual matrices

Description

Methods extending to dual matrices the corresponding methods for numeric matrices.

Usage

det.dual(x, ...)

## S4 method for signature 'dual'
det(x, ...)

## S3 method for class 'dual'
determinant(x, logarithm = TRUE, ...)

## S4 method for signature 'dual,dual'
solve(a, b, ...)

## S4 method for signature 'dual,missing'
solve(a, b, ...)

## S4 method for signature 'numericOrArray,dual'
solve(a, b, ...)

## S4 method for signature 'dual,numericOrArray'
solve(a, b, ...)

Arguments

x

a dual matrix

...

extra parameters (ignored)

logarithm

if 'TRUE', get logarithm of modulus of determinant

a, b

dual or numerical arguments for 'solve'

Value

'det' returns a dual scalar, 'determinant' a list with components 'modulus' (which is a dual object) and 'sign', and 'solve' returns a dual object (vector or matrix).

Examples

x <- dual( matrix(c(1,2,1,3), 2, 2) )
det(x)
d(det(x), "x1.1")
solve(x)
d(solve(x), "x1.1")