Skip to content

Commit

Permalink
Fixup
Browse files Browse the repository at this point in the history
  • Loading branch information
jeromekelleher committed Jan 10, 2025
1 parent 2f48e3f commit b407ed5
Showing 1 changed file with 2 additions and 9 deletions.
11 changes: 2 additions & 9 deletions bio2zarr/vcf2zarr/vcz.py
Original file line number Diff line number Diff line change
Expand Up @@ -843,17 +843,10 @@ def encode_local_alleles_partition(self, partition_index):
store=self.wip_partition_array_path(partition_index, "call_genotype"),
mode="r",
)
alleles_array = zarr.open_array(
store=self.wip_partition_array_path(partition_index, "variant_allele"),
mode="r",
)
for chunk_index in range(gt_array.cdata_shape[0]):
A = alleles_array.blocks[chunk_index]
G = gt_array.blocks[chunk_index]
for alleles, var in zip(A, G):
for genotypes in gt_array.blocks[chunk_index]:
j = call_LAA.next_buffer_row()
# TODO we should probably compute LAAs by chunk for efficiency
call_LAA.buff[j] = compute_laa_field(var, alleles)
call_LAA.buff[j] = compute_laa_field(genotypes)

call_LAA.flush()
self.finalise_partition_array(partition_index, "call_LAA")
Expand Down

0 comments on commit b407ed5

Please sign in to comment.