{- |
Module      :  ELynx.Import.MarkovProcess.EDMModelPhylobayesSpec
Copyright   :  (c) Dominik Schrempf 2019
License     :  GPL-3

Maintainer  :  dominik.schrempf@gmail.com
Stability   :  unstable
Portability :  portable

Creation date: Tue Jan 29 12:26:49 2019.

-}

module ELynx.Import.MarkovProcess.EDMModelPhylobayesSpec
  (spec) where

import           Numeric.LinearAlgebra                         (fromList)
import           Test.Hspec

import           ELynx.Import.MarkovProcess.EDMModelPhylobayes
import           ELynx.Tools.InputOutput                       (parseFileWith)
import           Files                                         (componentsFilePhylobayes)

testComponents :: [EDMComponent]
testComponents =
  [ (0.3690726059430788 , fromList
    [ 0.0746116859730418, 0.0069967050701822, 0.0792378063875672, 0.1103113280153337
    , 0.0084753541579630, 0.0420976864270153, 0.0359995156499732, 0.0176053658394712
    , 0.1140305648643406, 0.0306656386818212, 0.0143923751236750, 0.0805907035564168
    , 0.0136556223619922, 0.0854829959418090, 0.0907856579607629, 0.0955855590305961
    , 0.0585231905811694, 0.0276404860443685, 0.0024024774767569, 0.0109092808557433 ])
  , (0.28019163862430846, fromList
    [ 0.0379912838025543, 0.0085917039746034, 0.0016977963056655, 0.0029780993213938
    , 0.0369086689176780, 0.0051768704164023, 0.0027496797757250, 0.2650202435999584
    , 0.0035114356263963, 0.2705514660371282, 0.0634852940497682, 0.0036512748214048
    , 0.0030666766077233, 0.0046946399741140, 0.0039700770877914, 0.0137110414055372
    , 0.0286159990065166, 0.2351259204392804, 0.0018121721806795, 0.0066896566496796 ])
  , (0.23699880225859807, fromList
    [ 0.2699288544670601, 0.0262695884363299, 0.0119030286238640, 0.0143694121662917
    , 0.0113341089668980, 0.1028437949687162, 0.0097376980401260, 0.0250662077996463
    , 0.0137176088824149, 0.0302792626622677, 0.0151697996735645, 0.0255041240872322
    , 0.0194279535358810, 0.0155053062932145, 0.0146597721555602, 0.2261838872893110
    , 0.1040919090386891, 0.0534187313840882, 0.0024899475832214, 0.0080990039456234 ])
  , (0.11373695317401472, fromList
    [ 0.0325183213724537, 0.0153142424565790, 0.0072297704620506, 0.0101526783552517
    , 0.2682049848838039, 0.0107298905683658, 0.0420960604309369, 0.0414751553493385
    , 0.0131985876449809, 0.1152120909184748, 0.0300072241461393, 0.0173787316286536
    , 0.0062421326435569, 0.0173124967325860, 0.0181776303531701, 0.0328778350287987
    , 0.0250382111285986, 0.0437257936268565, 0.0258103804738129, 0.2272977817955918 ]) ]

spec :: Spec
spec = describe "phylobayes" $
  it "parses a text file with stationary distributions in phylobayes format" $ do
  cs <- parseFileWith phylobayes componentsFilePhylobayes
  cs `shouldBe` testComponents