R/gmedc.R

Defines functions gmedc

Documented in gmedc

gmedc <-
function(A,B=diag(nrow(A)))
  {
   ss<-medc(B,function(x) udtd(sqrt(x)))
   ev <- eigen(ss%*%A%*%ss)
   return(list(hfinv=ss,gvalues=ev$values,gvectors=ss%*%ev$vectors))
 }

Try the RKUM package in your browser

Any scripts or data that you put into this service are public.

RKUM documentation built on June 22, 2022, 9:06 a.m.