sweeper: Beaton's sweep function

In mnbem: Matrix Normal Block EM

Usage

sweeper(a, indi, rev = FALSE)

Arguments

a
indi
rev

Examples

##---- Should be DIRECTLY executable !! ----
##-- ==>  Define data, use random,
##--	or do  help(data=index)  for the standard data sets.

## The function is currently defined as
function( a, indi, rev=FALSE ){
    d = dim( a )
    if( d[1]!=d[2] ) stop( "Only works on square matrix" )
    if( rev ){ rv = -1 } else { rv = 1 }
    z <- .Call( "SWEEP",
                as.integer( indi-1 ),
                as.double( a ),
                as.integer( d[1] ),
                as.double( rv )
              );
    return( array( z, d ) );
  }

mnbem documentation built on May 2, 2019, 6:16 p.m.