Skip to content

Commit

Permalink
include mcmcio file
Browse files Browse the repository at this point in the history
  • Loading branch information
shajoezhu committed Apr 16, 2020
1 parent 113078a commit f63516c
Show file tree
Hide file tree
Showing 3 changed files with 35 additions and 8 deletions.
1 change: 1 addition & 0 deletions Makefile.am
Original file line number Diff line number Diff line change
Expand Up @@ -33,6 +33,7 @@ common_src = src/random/fastfunc.cpp \
src/vcfReader.cpp \
src/variantIndex.cpp \
src/gzstream/gzstream.cpp \
src/export/mcmcIO.cpp \
src/export/dEploidIOExport.cpp \
src/export/dEploidIOExportPosteriorProb.cpp \
src/export/writeMcmcRelated.cpp \
Expand Down
33 changes: 33 additions & 0 deletions src/export/mcmcIO.cpp
Original file line number Diff line number Diff line change
@@ -0,0 +1,33 @@
/*
* dEploid is used for deconvoluting Plasmodium falciparum genome from
* mix-infected patient sample.
*
* Copyright (C) 2016-2017 University of Oxford
*
* Author: Sha (Joe) Zhu
*
* This file is part of dEploid.
*
* dEploid is free software: you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation, either version 3 of the License, or
* (at your option) any later version.
*
* This program is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with this program. If not, see <http://www.gnu.org/licenses/>.
*
*/

#include "mcmc.hpp"

void McmcMachinery::printArray(vector <double> array) {
for (auto const& value : array) {
cout << value << " ";
}
cout << endl;
}
9 changes: 1 addition & 8 deletions src/mcmc.hpp
Original file line number Diff line number Diff line change
Expand Up @@ -173,13 +173,7 @@ class McmcMachinery {
double calcLogPriorTitre(const vector <double> &tmpTitre);
double rBernoulli(double p);

void printArray(vector <double> array) {
for (auto const& value : array) {
cout << value << " ";
}
cout << endl;
}

void printArray(vector <double> array);
void sampleMcmcEvent(bool useIBD = false);
void recordMcmcMachinery();
bool recordingMcmcBool_;
Expand Down Expand Up @@ -231,4 +225,3 @@ class McmcMachinery {
};

#endif

0 comments on commit f63516c

Please sign in to comment.