haplotypes_phase needs to be set to NULL if not initialized

bmvdgeijn · Sep 13, 2018 · de82c31 · de82c31
1 parent 3f6a327
commit de82c31
Showing 1 changed file with 7 additions and 5 deletions.
diff --git a/snp2h5/snp2h5.c b/snp2h5/snp2h5.c
@@ -1169,14 +1169,16 @@ void parse_vcf(Arguments *args, Chromosome *all_chroms, int n_chrom,
 					 vcf->sample_names, vcf->n_sample);
 	sample_tab_free(samp_tab);
 
-  /* initialize the phase matrix */
-  haplotypes_phase = my_malloc(vcf->n_sample * sizeof(char));
-  init_h5matrix(&haplotype_phase_info, haplotype_info->n_row,
-          vcf->n_sample,
-          HAPLOTYPE_DATATYPE, util_str_concat("phase_", chrom->name, NULL));
+	/* initialize the phase matrix */
+	haplotypes_phase = my_malloc(vcf->n_sample * sizeof(char));
+	init_h5matrix(&haplotype_phase_info, haplotype_info->n_row,
+		      vcf->n_sample,
+		      HAPLOTYPE_DATATYPE,
+		      util_str_concat("phase_", chrom->name, NULL));
       }
     } else {
       haplotypes = NULL;
+      haplotypes_phase = NULL;
     }
     if(args->snp_index_file) {
       /* SNP index vector is same length as chromosome,