Skip to content

Commit

Permalink
tiny update
Browse files Browse the repository at this point in the history
  • Loading branch information
Zilong-Li committed Feb 11, 2025
1 parent 6bba622 commit 40c6771
Show file tree
Hide file tree
Showing 3 changed files with 4 additions and 6 deletions.
4 changes: 2 additions & 2 deletions Makefile
Original file line number Diff line number Diff line change
Expand Up @@ -24,7 +24,7 @@ IOMP5 = 0

########################### end ###########################

VERSION=0.4.9
VERSION=0.5.0
# detect OS architecture and add flags
Platform := $(shell uname -s)

Expand Down Expand Up @@ -135,7 +135,7 @@ else ifeq ($(Platform),Darwin)
MYFLAGS += -DEIGEN_USE_BLAS
INC += -I${ACCELERATE_ROOT}/include
LPATHS += -L${ACCELERATE_ROOT}/lib
DLIBS += -Wl,-rpath,${ACCELERATE_ROOT}/lib -llapack -lblas -lgomp -lpthread
DLIBS += -Wl,-rpath,${ACCELERATE_ROOT}/lib -llapack -lblas -lpthread

else
DLIBS += -lomp -lpthread -L$(LIB_PATH)
Expand Down
4 changes: 1 addition & 3 deletions src/Data.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -171,9 +171,7 @@ void Data::write_eigs_files(const Mat1D &S, const Mat2D &U, const Mat2D &V) {
std::ofstream outs(params.fileout + ".eigvals");
std::ofstream outu(params.fileout + ".eigvecs");
Eigen::IOFormat fmt(6, Eigen::DontAlignCols, "\t", "\n");
if (outs.is_open()) {
outs << (S * params.ploidy).format(fmt) << '\n';
}
if (outs.is_open()) outs << (S * params.ploidy).format(fmt) << '\n';
if (outu.is_open()) outu << U.format(fmt) << '\n';
if (params.project == 0 && params.printv) {
save_snps_in_bim();
Expand Down
2 changes: 1 addition & 1 deletion src/Halko.cpp
Original file line number Diff line number Diff line change
Expand Up @@ -173,7 +173,7 @@ void FancyRsvdOpData::computeGandH(Mat2D& G, Mat2D& H, int pi) {
// blocksize: how many snps in each block
blocksize = (unsigned int)ceil((double)data->nsnps / data->params.bands);
if (blocksize < data->params.bands)
cao.warn("block size < window size. please consider the IRAM method with --svd 1");
cao.warn("block size < window size. please consider the IRAM method with --svd 0");
if (data->params.perm) {
cao.print(tick.date(), "permuting data matrix by columns in place");
PCAone::permute_matrix(data->G, data->perm);
Expand Down

0 comments on commit 40c6771

Please sign in to comment.