#    source("RFtest.model.R")

if (EXTENDED.TESTING <- file.exists("source.R")) { source("source.R")
} else if (file.exists(f <- "~/R/RF/RandomFields/tests/source.R")) source(f) 

x <- function(...) {
  print(unlist(list(...)))
  str(PrepareModel(...))
  cat("--------------------------------\n")
  str(m <- PrepareModel(...))
  invisible(m)
}
#x <- function(...) str(PrepareModel(...))

# stop("")
x(model="nugg",param=c(0,2,3,4))
x(model="nugg",param=c(0,0,3,4))
x(model="nugg",param=c(0,3,0,4))
x(model="nugg",param=c(0,3,0,4),me="ci")

x(model="sphe",param=c(1,2,0,4))
x(model="sphe",param=c(1,2,3,4))
x(model="sphe",param=c(1,2,3,4), me="ci")

x(model="sphe",param=cbind(c(2,4),c(8,1)), me="ci")
x(model="sphe",param=cbind(c(2,4),c(3,0)), me="ci")
x(model="sphe",param=cbind(c(2,4),c(3,0)), trend=1, me="ci")
x(model="sphe",param=cbind(c(2,4),c(3,0)), trend=c(1,2),me="ci")#not simple
x(model="sphe",param=cbind(c(2,4),c(3,0),c(7,0)), me="ci")
x(model="sphe",param=cbind(c(2,4),c(3,0),c(8,1)), me="ci")
x(model="sphe",param=cbind(c(2,4),c(3,0),c(8,1)), trend=1)

x(model="whi",param=c(1,2,3,4,5))
x(model="whi",param=c(1,2,0,4,5))
x(model="whi",param=c(1,2,3,4,5), me="ci")
x(model="whi",param=c(1,2,0,4,5), me="ci")

x(model="whi",param=cbind(c(2,4,5),c(8,1,1)), me="ci")
x(model="whi",param=cbind(c(2,4,5),c(3,0,0)), me="ci", trend=9)
x(model="whi",param=cbind(c(2,4,5),c(3,0,0),c(7,0,0)), me="ci")
x(model="whi",param=cbind(c(2,4,5),c(3,0,0),c(8,1,1)), me="ci")

# list to old, vector
x(model=list(list(model="whi",k=5,var=2,s=4)), me="ci")
x(model=list(list(model="whi",k=5,var=2,s=4),"+",
    list(model="whi",k=1,var=3,s=0)
    ),
   me="ci")
## note difference to the next result !!!!!!!!!
xx <- PrepareModel(model= list(list(model="whi",k=5,var=2,s=4),
                     "+",
                     list(model="whi",k=1,var=3,s=0)),
                   me="ci")
x(model=xx$mo, pa=xx$pa,  me=xx$me)


# list to nested
x(model=list(list(model="whi",k=5,var=2,s=4),"+",
    list(model="nugg",var=3,s=0),"+",
    list(model="nugg",var=7,s=0)
    ),
  me="ci")

x(model=list(list(model="whi",k=5,var=2,s=4),"+",
    list(model="whi",k=1,var=6,s=1),"+",
    list(model="nugg",var=3,s=0)
   ),
  me="ci")

x(model=list(list(model="whi",k=5,var=2,s=4),"+",
    list(model="whi",k=1,var=6,s=1),"+",
    list(model="nugg",var=3,s=0)
   ),
  me="ci")


# multiplicative model
x(model=list(list(model="whi",k=5,var=2,s=4),"+",
    list(model="whi",k=9,var=6,s=1),"*",
    list(model="nugg",var=3,s=0)
   ),
  me="ci")

## anisotropy
x(model=list(list(model="whi",k=5,var=2,a=1:4),"+",
    list(model="whi",k=1,var=6,a=2:5),"+",
    list(model="nugg",var=3,a=3:6)
   ),
  me="ci")

# geht auch
x(model=list(list(model="whi",k=5,var=2,s=4),"+",
    list(model="whi",k=1,var=6,s=1),"*",
    list(model="nugg",var=3,s=0)
   ),
   me="ci")