diff --git a/miscprob.cc b/miscprob.cc index 276de51be2ece567ae29d17b5f483b9f90befcea..37396b1ff92e85273b6529e1ef55c9b1320afedb 100644 --- a/miscprob.cc +++ b/miscprob.cc @@ -171,6 +171,7 @@ ReturnMatrix mvnrnd(const RowVector& mu, const SymmetricMatrix& covar, int nsamp return mvn.next(nsamp); } + /* // Saad: Wishart and inverseWishart Random Generator ReturnMatrix wishrnd(const SymmetricMatrix& sigma,const int dof){ // compute cholesky factor for sigma @@ -216,6 +217,8 @@ ReturnMatrix iwishrnd(const SymmetricMatrix& sigma,const int dof){ return res; } + */ + ReturnMatrix perms(const int n){ if(n<=1){ Matrix P(1,1); diff --git a/miscprob.h b/miscprob.h index 127de4f45bfdb01b1414e1be7f8e825b48406366..d661adf165811816371aa5dae7e468818e90e09a 100644 --- a/miscprob.h +++ b/miscprob.h @@ -17,7 +17,7 @@ #include "stdlib.h" using namespace NEWMAT; -using namespace NEWRAN; +// using namespace NEWRAN; namespace MISCMATHS { @@ -49,8 +49,8 @@ namespace MISCMATHS { // returns n! * n matrix of all possible permutations ReturnMatrix perms(const int n); - ReturnMatrix wishrnd(const SymmetricMatrix&,const int); - ReturnMatrix iwishrnd(const SymmetricMatrix&,const int); + // ReturnMatrix wishrnd(const SymmetricMatrix&,const int); + // ReturnMatrix iwishrnd(const SymmetricMatrix&,const int); class Mvnormrandm {