context("NMreadInits")

### I had to run these with as.fun=data.table and skip testing of
### lines$text.before to make them work on githubs different
### arcitectures. I am not sure if the data.table step is needed, but
### I think it was. Inserted this to drop testing text.before.
##   res1$lines[,text.before:=""]
### the reg expressions must evaluate differently on different architectures


test_that("basic",{

    fileRef <- "testReference/NMreadInits_01.rds"
    file.mod <- "testData/nonmem/xgxr032.mod"
    
    NMdataConf(reset=T)
    NMdataConf(as.fun="data.table")

    lines <- readLines(file.mod)
    res1  <- NMreadInits(lines=lines,return="all")

    ## res1$lines[grepl("^ +$",text.before),text.before:=""]
    res1$lines[,text.before:=""]
    
    ## cat(paste(res1$lines$text.before,collapse=":"),"\n")
    ## cat(paste(readRDS(fileRef)$lines$text.before,collapse=":"),"\n")

    ## expect_equal_to_reference(res1,fileRef)
    expect_equal_to_reference(res1$pars,fnAppend(fileRef,"pars"))
    if(F){
        ref.pars <- readRDS(fnAppend(fileRef,"pars"))
        ref.pars
        res1$pars
    }

    expect_equal_to_reference(res1$elements,fnAppend(fileRef,"elems"))
    expect_equal_to_reference(res1$lines,fnAppend(fileRef,"lines"))


    if(FALSE){

        ref.elems <- readRDS(fnAppend(fileRef,"elems"))
        compareCols(ref.elems,res1$elements)
        ref.elems
        res1$elements

    }
    
})

test_that("with OMEGA block",{
    NMdataConf(reset=T)
    NMdataConf(as.fun="data.table")

    fileRef <- "testReference/NMreadInits_02.rds"

    file.mod <- "testData/nonmem/xgxr133.mod"
    res1 <- NMreadInits(file=file.mod,return="all")

    ## res1$lines[grepl("^ +$",text.before),text.before:=""]
    res1$lines[,text.before:=""]

    ## expect_equal_to_reference(res1,fileRef)
    expect_equal_to_reference(res1$pars,fnAppend(fileRef,"pars"))
    if(F){
        ref.pars <- readRDS(fnAppend(fileRef,"pars"))
        ref.pars
        res1$pars
    }

    expect_equal_to_reference(res1$elements,fnAppend(fileRef,"elems"))
    expect_equal_to_reference(res1$lines,fnAppend(fileRef,"lines"))


    if(FALSE){
        ref.elems <- readRDS(fnAppend(fileRef,"elems"))
        compareCols(ref.elems,res1$elements)
        ref.elems
        res1$elements
        ##expect_equal_to_reference(res1$elements[,!("parameter")],fnAppend(fileRef,"elems"))    
    }


})

test_that("OMEGA SAME",{
    NMdataConf(reset=T)
    NMdataConf(as.fun="data.table")

    fileRef <- "testReference/NMreadInits_03.rds"
    
    text <- c("
$THETA
(0,0.1) ; THE1      - 30) 1st theta
 (0,4.2) ; THE2        - 31) 2nd theta
$OMEGA  0.08   ;    IIV.TH1  ; 1  ;IIV
$OMEGA  BLOCK(1)
 0.547465  ; IOV.TH1  ; 2 ;IOV
$OMEGA  BLOCK(1) SAME
$OMEGA  BLOCK(1) SAME")

    res1 <- NMreadInits(lines=text,return="all")
    res1$lines[,text.before:=""]

    ## expect_equal_to_reference(res1,fileRef)
    expect_equal_to_reference(res1$pars,fnAppend(fileRef,"pars"))
    if(F){
        ref.pars <- readRDS(fnAppend(fileRef,"pars"))
        ref.pars
        res1$pars
    }

    expect_equal_to_reference(res1$elements,fnAppend(fileRef,"elems"))
    expect_equal_to_reference(res1$lines,fnAppend(fileRef,"lines"))


    if(FALSE){
        ref.elems <- readRDS(fnAppend(fileRef,"elems"))
        ref.elems
        res1$elements
        compareCols(ref.elems,res1$elements)
        ##expect_equal_to_reference(res1$elements[,!("parameter")],fnAppend(fileRef,"elems"))

        ref.pars <- readRDS(fnAppend(fileRef,"pars"))
        compareCols(ref.pars,res1$pars)
        ref.pars
        res1$pars
    }


})


test_that("OMEGA BLOCK FIX 1",{
    NMdataConf(reset=T)
    NMdataConf(as.fun="data.table")

    fileRef <- "testReference/NMreadInits_04.rds"
    
    text1 <- c("
$THETA
(0,0.1) ; THE1      - 30) 1st theta
$OMEGA  BLOCK(2) 1 FIX  .3 ; first omega
 0.547465  ; IOV.TH1  ; 2 ;IOV
")

    text2 <- c("
$THETA
(0,0.1) ; THE1      - 30) 1st theta
$OMEGA  BLOCK(2) FIX
1 .3 ; first omega
 0.547465  ; IOV.TH1  ; 2 ;IOV
")

    ## options(warn=0)
    res1 <- NMreadInits(lines=text1,return="all",as.fun="data.table")
    res2 <- NMreadInits(lines=text2,return="all",as.fun="data.table")

    ## res1$lines
    ## res2$lines

    expect_equal(res1$pars,
                 res2$pars)
    
})


test_that("BLOCK(2) SAME",{

    fileRef <- "testReference/NMreadInits_05.rds"

#### 

    NMdataConf(reset=T)
    NMdataConf(as.fun="data.table")


    text <- c("
$THETA
(0,0.1) ;  1 ; 1st theta
 (0,4.2) ; 2 ; 2nd theta
$OMEGA  0.08   ;    IIV.TH1  ; 1  ;IIV
;; A block(2) + SAME
$OMEGA  BLOCK(2)
 0.547465  .01 .3     ; 2+3
$OMEGA  BLOCK(2) SAME ; 4+5
")
    as.NMctl(text,lines=T)


    res1  <- NMreadInits(lines=text,return="all")

###  elemnum should not inc because there is only one observed element
###  (SAME only written once in control stream)
    res1

    expect_equal_to_reference(res1$pars,
                              fileRef)


})