molstar
Version:
A comprehensive macromolecular library.
87 lines (86 loc) • 3.22 kB
JavaScript
/**
* Copyright (c) 2021 mol* contributors, licensed under MIT, See LICENSE file for more info.
*
* @author David Sehnal <david.sehnal@gmail.com>
*/
import { Column, Table } from '../../mol-data/db';
import { Task } from '../../mol-task';
import { createModels } from './basic/parser';
import { BasicSchema, createBasic } from './basic/schema';
import { ComponentBuilder } from './common/component';
import { EntityBuilder } from './common/entity';
function getModels(mol, ctx) {
const { molecules } = mol;
let count = 0;
for (const m of molecules)
count += m.count;
const type_symbols = new Array(count);
const id = new Int32Array(count);
const x = new Float32Array(count);
const y = new Float32Array(count);
const z = new Float32Array(count);
const model_num = new Int32Array(count);
let offset = 0;
for (let i = 0; i < molecules.length; i++) {
const m = molecules[i];
for (let j = 0; j < m.count; j++) {
type_symbols[offset] = m.type_symbol.value(j);
x[offset] = m.x.value(j);
y[offset] = m.y.value(j);
z[offset] = m.z.value(j);
id[offset] = j;
model_num[offset] = i;
offset++;
}
}
const MOL = Column.ofConst('MOL', count, Column.Schema.str);
const A = Column.ofConst('A', count, Column.Schema.str);
const seq_id = Column.ofConst(1, count, Column.Schema.int);
const type_symbol = Column.ofStringArray(type_symbols);
const atom_site = Table.ofPartialColumns(BasicSchema.atom_site, {
auth_asym_id: A,
auth_atom_id: type_symbol,
auth_comp_id: MOL,
auth_seq_id: seq_id,
Cartn_x: Column.ofFloatArray(x),
Cartn_y: Column.ofFloatArray(y),
Cartn_z: Column.ofFloatArray(z),
id: Column.ofIntArray(id),
label_asym_id: A,
label_atom_id: type_symbol,
label_comp_id: MOL,
label_seq_id: seq_id,
label_entity_id: Column.ofConst('1', count, Column.Schema.str),
occupancy: Column.ofConst(1, count, Column.Schema.float),
type_symbol,
pdbx_PDB_model_num: Column.ofIntArray(model_num),
}, count);
const entityBuilder = new EntityBuilder();
entityBuilder.setNames([['MOL', 'Unknown Entity']]);
entityBuilder.getEntityId('MOL', 0 /* MoleculeType.Unknown */, 'A');
const componentBuilder = new ComponentBuilder(seq_id, type_symbol);
componentBuilder.setNames([['MOL', 'Unknown Molecule']]);
componentBuilder.add('MOL', 0);
const basic = createBasic({
entity: entityBuilder.getEntityTable(),
chem_comp: componentBuilder.getChemCompTable(),
atom_site
});
return createModels(basic, XyzFormat.create(mol), ctx);
}
//
export { XyzFormat };
var XyzFormat;
(function (XyzFormat) {
function is(x) {
return (x === null || x === void 0 ? void 0 : x.kind) === 'xyz';
}
XyzFormat.is = is;
function create(mol) {
return { kind: 'xyz', name: 'xyz', data: mol };
}
XyzFormat.create = create;
})(XyzFormat || (XyzFormat = {}));
export function trajectoryFromXyz(mol) {
return Task.create('Parse XYZ', ctx => getModels(mol, ctx));
}