/** * Copyright (c) 2018 mol* contributors, licensed under MIT, See LICENSE file for more info. * * @author David Sehnal */ import { createRangeArray, makeBuckets } from '../../../mol-data/util'; import { Column, Table } from '../../../mol-data/db'; import { RuntimeContext } from '../../../mol-task'; import { AtomSite } from './schema'; import { arrayIsIdentity } from '../../../mol-util/array'; export type SortedAtomSite = { atom_site: AtomSite sourceIndex: Column } export async function sortAtomSite(ctx: RuntimeContext, atom_site: AtomSite, start: number, end: number): Promise { const indices = createRangeArray(start, end - 1); const { label_entity_id, label_asym_id, label_seq_id } = atom_site; const entityBuckets = makeBuckets(indices, label_entity_id.value); if (ctx.shouldUpdate) await ctx.update(); for (let ei = 0, _eI = entityBuckets.length - 1; ei < _eI; ei++) { const chainBuckets = makeBuckets(indices, label_asym_id.value, { start: entityBuckets[ei], end: entityBuckets[ei + 1] }); for (let cI = 0, _cI = chainBuckets.length - 1; cI < _cI; cI++) { const aI = chainBuckets[cI]; // are we in HETATM territory? if (label_seq_id.valueKind(aI) !== Column.ValueKinds.Present) continue; makeBuckets(indices, label_seq_id.value, { sort: true, start: aI, end: chainBuckets[cI + 1] }); if (ctx.shouldUpdate) await ctx.update(); } if (ctx.shouldUpdate) await ctx.update(); } if (arrayIsIdentity(indices) && indices.length === atom_site._rowCount) { return { atom_site, sourceIndex: Column.ofIntArray(indices) }; } return { atom_site: Table.view(atom_site, atom_site._schema, indices) as AtomSite, sourceIndex: Column.ofIntArray(indices) }; }