#####################################
blas_set_num_threads(1)
include("Boot_code.jl")
n1=50
n2=50
p=200 
m0 = ifloor((n1+n2-2)*.5)
m = ifloor([.5*m0, m0, 1.5*m0, 2*m0-5])
k = [0,1]

kap0 = [.1,1,10]
mu=1

pct = [1.0,.5,.95]
parm = hcat(map(x->vcat(x...),product(n1,n2,p,k,m,kap0,pct,mu))...)'  ### all parameters combinasion
#parm = hcat(map(x->vcat(x...),product(n1,n2,p,k,m,kap0,.01,0.0))...)'  ### all parameters combinasion
length(m)
size(parm)

## This part works pretty well no problem

out = BayesFac_Block_noloc(parm[1,:])
out

## But the part involving pmap does not work work at all

println("Start here ... ") 
  addprocs(2)
  M = { parm[i,:] for i in 1:2} # size(parm,1)}
  require("Boot_code.jl")
  ot = pmap(BayesFac_Block_noloc,M)
  println("Done!!") 

 ot 
 outf=ot[1]
 for i in 2:(size(ot,1))
 outf = vcat(outf,ot[i])
 end 
 data = convert(DataFrame,outf)
 names!(data,[:n1,:n2,:p,:k,:m,:kap0,:pct,:mu,:Log_bay,:pos,:stron,:ver_strg])