inst/test_objs/SFTEVAL_test.R

In Auburngrads/SMRD: Statistical Methods For Reliability Data

library(smrdfortran)
kdmod = 1
mu1 = c(0.5, 0.25, 0.015, 5)
sig1 = 0.25
mu2 = 1.1
sig2  = 0.55
rho = .5
df = -0.5
d0 = 0
sfact = 0
tf = 1
kprint = 0
debug1 = F

  max.length <- max(length(mu1) , length(sig1) , length(mu2),
                    length(sig2), length(rho)  , length(df),
                    length(d0)  , length(sfact), length(tf))
  mu1    <- smrdfortran:::expand.vec(mu1, max.length)
  sig1   <- smrdfortran:::expand.vec(sig1, max.length)
  mu2    <- smrdfortran:::expand.vec(mu2, max.length)
  sig2   <- smrdfortran:::expand.vec(sig2, max.length)
  rho    <- smrdfortran:::expand.vec(rho, max.length)
  df     <- smrdfortran:::expand.vec(df, max.length)
  d0     <- smrdfortran:::expand.vec(d0, max.length)
  sfact  <- smrdfortran:::expand.vec(sfact, max.length)
  tf     <- smrdfortran:::expand.vec(tf, max.length)
  zout <- .Fortran("sfteval", as.integer(kdmod), as.double(mu1),
                   as.double(sig1), as.double(mu2), as.double(sig2), as.double(rho),
                   as.double(df), as.double(d0), as.double(sfact), as.double(tf),
                   as.integer(max.length), answer = double(max.length),
                   ier = integer(max.length), as.integer(kprint))
  new <- SFTEVAL(as.integer(kdmod), 
                 as.double(mu1),
                 as.double(sig1), 
                 as.double(mu2), 
                 as.double(sig2), 
                 as.double(rho),
                 as.double(df), 
                 as.double(d0), 
                 as.double(sfact), 
                 as.double(tf),
                 as.integer(max.length), 
                 answer = double(max.length),
                 ier = integer(max.length), 
                 as.integer(kprint))

  sfteval_test <- function(kdmod = 1,
                           mu1 = 0.5,
                           sig1 = 0.25,
                           mu2 = 1,
                           sig2 = 0.5,
                           rho = 0.5,
                           df = -0.5,
                           d0 = 0,
                           sfact = 0,
                           tf = 1,
                           kprint = 0){

    max.length <- max(length(mu1) , length(sig1) , length(mu2),
                    length(sig2), length(rho)  , length(df),
                    length(d0)  , length(sfact), length(tf))
  mu1    <- smrdfortran:::expand.vec(mu1, max.length)
  sig1   <- smrdfortran:::expand.vec(sig1, max.length)
  mu2    <- smrdfortran:::expand.vec(mu2, max.length)
  sig2   <- smrdfortran:::expand.vec(sig2, max.length)
  rho    <- smrdfortran:::expand.vec(rho, max.length)
  df     <- smrdfortran:::expand.vec(df, max.length)
  d0     <- smrdfortran:::expand.vec(d0, max.length)
  sfact  <- smrdfortran:::expand.vec(sfact, max.length)
  tf     <- smrdfortran:::expand.vec(tf, max.length)

  new <- SMRD::sfteval(as.integer(kdmod),
                         as.double(mu1),
                         as.double(sig1),
                         as.double(mu2),
                         as.double(sig2),
                         as.double(rho),
                         as.double(df),
                         as.double(d0),
                         as.double(sfact),
                         as.double(tf),
                         as.integer(max.length),
                         answer = double(max.length),
                         ier = integer(max.length),
                         as.integer(kprint))
  old <- smrdfortran:::fteval(kdmod,
                       mu1,
                       sig1,
                       mu2,
                       sig2,
                       rho,
                       df,
                       d0,
                       sfact,
                       tf,
                       kprint,
                       debug1= F)

    return(abs(old - new[,1]))

                           }