debug ld clump

src/LD.cpp

@@ -74,6 +74,7 @@ std::tuple<Int2D, Int2D> get_target_snp_idx(const SNPld& snp_t,
     i = 0;
     while (snp_t.chr[i] != snp.chr[c]) {
       i++;
+      if (i == (int)snp_t.chr.size()) break;
     }
     ord.push_back(i);
   }

src/Utils.hpp

@@ -51,9 +51,9 @@ struct SNPld {
   std::vector<int> pos;          // pos of each SNP
   std::vector<int> end_pos;      // 0-based index for last snp pos
   std::vector<std::string> chr;  // chr sequences
-  std::vector<int> ws, we;           //  the snp index, i.e the index for lead SNP
-  // std::vector<int> we;  // the number of SNPs (including lead SNP) in a window
-  Double1D af;          // allele frequency
+  Double1D af;                   // allele frequency
+  std::vector<int> ws;           //  the snp index, i.e the index for lead SNP
+  std::vector<int> we;  // the number of SNPs (including lead SNP) in a window
 };
 
 std::string get_machine();