/** * Copyright (c) 2018 Mol* contributors, licensed under MIT, See LICENSE file for more info. * * @author David Sehnal <david.sehnal@gmail.com> */ import { MolScriptSymbolTable as MolScript } from '../../language/symbol-table'; import { DefaultQueryRuntimeTable, QuerySymbolRuntime, QueryRuntimeArguments } from './compiler'; import { Queries, StructureProperties, StructureElement, QueryContext } from 'mol-model/structure'; import { ElementSymbol } from 'mol-model/structure/model/types'; import { SetUtils } from 'mol-util/set'; import toUpperCase from 'mol-util/upper-case'; import { VdwRadius, AtomWeight, AtomNumber } from 'mol-model/structure/model/properties/atomic'; import { cantorPairing } from 'mol-data/util'; import C = QuerySymbolRuntime.Const import D = QuerySymbolRuntime.Dynamic const symbols = [ // ============= TYPES ============= C(MolScript.core.type.bool, (ctx, v) => !!v[0](ctx)), C(MolScript.core.type.num, (ctx, v) => +v[0](ctx)), C(MolScript.core.type.str, (ctx, v) => '' + v[0](ctx)), C(MolScript.core.type.list, (ctx, xs) => QueryRuntimeArguments.forEachEval(xs, ctx, (v, i, list) => list[i] = v, [] as any[])), C(MolScript.core.type.set, (ctx, xs) => QueryRuntimeArguments.forEachEval(xs, ctx, (v, i, set) => set.add(v), new Set<any>())), C(MolScript.core.type.regex, (ctx, v) => new RegExp(v[0](ctx), (v[1] && v[1](ctx)) || '')), C(MolScript.core.type.bitflags, (ctx, v) => +v[0](ctx)), C(MolScript.core.type.compositeKey, (ctx, xs) => QueryRuntimeArguments.forEachEval(xs, ctx, (v, i, list) => list[i] = '' + v, [] as string[]).join('-')), // ============= LOGIC ================ C(MolScript.core.logic.not, (ctx, v) => !v[0](ctx)), C(MolScript.core.logic.and, (ctx, xs) => { if (typeof xs.length === 'number') { for (let i = 0, _i = xs.length; i < _i; i++) if (!xs[i](ctx)) return false; } else { for (const k of Object.keys(xs)) if (!xs[k](ctx)) return false; } return true; }), C(MolScript.core.logic.or, (ctx, xs) => { if (typeof xs.length === 'number') { for (let i = 0, _i = xs.length; i < _i; i++) if (xs[i](ctx)) return true; } else { for (const k of Object.keys(xs)) if (xs[k](ctx)) return true; } return false; }), // ============= RELATIONAL ================ C(MolScript.core.rel.eq, (ctx, v) => v[0](ctx) === v[1](ctx)), C(MolScript.core.rel.neq, (ctx, v) => v[0](ctx) !== v[1](ctx)), C(MolScript.core.rel.lt, (ctx, v) => v[0](ctx) < v[1](ctx)), C(MolScript.core.rel.lte, (ctx, v) => v[0](ctx) <= v[1](ctx)), C(MolScript.core.rel.gr, (ctx, v) => v[0](ctx) > v[1](ctx)), C(MolScript.core.rel.gre, (ctx, v) => v[0](ctx) >= v[1](ctx)), C(MolScript.core.rel.inRange, (ctx, v) => { const x = v[0](ctx); return x >= v[1](ctx) && x <= v[2](ctx); }), // ============= ARITHMETIC ================ C(MolScript.core.math.add, (ctx, xs) => { let ret = 0; if (typeof xs.length === 'number') { for (let i = 0, _i = xs.length; i < _i; i++) ret += xs[i](ctx); } else { for (const k of Object.keys(xs)) ret += xs[k](ctx); } return ret; }), C(MolScript.core.math.sub, (ctx, xs) => { let ret = 0; if (typeof xs.length === 'number') { if (xs.length === 1) return -xs[0](ctx); ret = xs[0](ctx) || 0; for (let i = 1, _i = xs.length; i < _i; i++) ret -= xs[i](ctx); } else { const keys = Object.keys(xs); if (keys.length === 1) ret = xs[keys[0]](ctx) || 0; for (let i = 1, _i = keys.length; i < _i; i++) ret -= xs[keys[i]](ctx); } return ret; }), C(MolScript.core.math.mult, (ctx, xs) => { let ret = 1; if (typeof xs.length === 'number') { for (let i = 0, _i = xs.length; i < _i; i++) ret *= xs[i](ctx); } else { for (const k of Object.keys(xs)) ret *= xs[k](ctx); } return ret; }), C(MolScript.core.math.div, (ctx, v) => v[0](ctx) / v[1](ctx)), C(MolScript.core.math.pow, (ctx, v) => Math.pow(v[0](ctx), v[1](ctx))), C(MolScript.core.math.mod, (ctx, v) => v[0](ctx) % v[1](ctx)), C(MolScript.core.math.min, (ctx, xs) => { let ret = Number.POSITIVE_INFINITY; if (typeof xs.length === 'number') { for (let i = 0, _i = xs.length; i < _i; i++) ret = Math.min(xs[i](ctx), ret); } else { for (const k of Object.keys(xs)) ret = Math.min(xs[k](ctx), ret) } return ret; }), C(MolScript.core.math.max, (ctx, xs) => { let ret = Number.NEGATIVE_INFINITY; if (typeof xs.length === 'number') { for (let i = 0, _i = xs.length; i < _i; i++) ret = Math.max(xs[i](ctx), ret); } else { for (const k of Object.keys(xs)) ret = Math.max(xs[k](ctx), ret) } return ret; }), C(MolScript.core.math.floor, (ctx, v) => Math.floor(v[0](ctx))), C(MolScript.core.math.ceil, (ctx, v) => Math.ceil(v[0](ctx))), C(MolScript.core.math.roundInt, (ctx, v) => Math.round(v[0](ctx))), C(MolScript.core.math.abs, (ctx, v) => Math.abs(v[0](ctx))), C(MolScript.core.math.sqrt, (ctx, v) => Math.sqrt(v[0](ctx))), C(MolScript.core.math.sin, (ctx, v) => Math.sin(v[0](ctx))), C(MolScript.core.math.cos, (ctx, v) => Math.cos(v[0](ctx))), C(MolScript.core.math.tan, (ctx, v) => Math.tan(v[0](ctx))), C(MolScript.core.math.asin, (ctx, v) => Math.asin(v[0](ctx))), C(MolScript.core.math.acos, (ctx, v) => Math.acos(v[0](ctx))), C(MolScript.core.math.atan, (ctx, v) => Math.atan(v[0](ctx))), C(MolScript.core.math.sinh, (ctx, v) => Math.sinh(v[0](ctx))), C(MolScript.core.math.cosh, (ctx, v) => Math.cosh(v[0](ctx))), C(MolScript.core.math.tanh, (ctx, v) => Math.tanh(v[0](ctx))), C(MolScript.core.math.exp, (ctx, v) => Math.exp(v[0](ctx))), C(MolScript.core.math.log, (ctx, v) => Math.log(v[0](ctx))), C(MolScript.core.math.log10, (ctx, v) => Math.log10(v[0](ctx))), C(MolScript.core.math.atan2, (ctx, v) => Math.atan2(v[0](ctx), v[1](ctx))), // ============= STRING ================ C(MolScript.core.str.match, (ctx, v) => v[0](ctx).test(v[1](ctx))), C(MolScript.core.str.concat, (ctx, xs) => { let ret: string[] = []; if (typeof xs.length === 'number') { for (let i = 0, _i = xs.length; i < _i; i++) ret.push(xs[i](ctx).toString()); } else { for (const k of Object.keys(xs)) ret.push(xs[k](ctx).toString()); } return ret.join(''); }), // ============= LIST ================ C(MolScript.core.list.getAt, (ctx, v) => v[0](ctx)[v[1](ctx)]), // ============= SET ================ C(MolScript.core.set.has, (ctx, v) => v[0](ctx).has(v[1](ctx))), C(MolScript.core.set.isSubset, (ctx, v) => SetUtils.isSuperset(v[1](ctx) as Set<any>, v[0](ctx) as Set<any>)), // ============= FLAGS ================ C(MolScript.core.flags.hasAny, (ctx, v) => { const test = v[1](ctx); const tested = v[0](ctx); if (!test) return !!tested; return (tested & test) !== 0; }), C(MolScript.core.flags.hasAll, (ctx, v) => { const test = v[1](ctx); const tested = v[0](ctx); if (!test) return !tested; return (tested & test) === test; }), //////////////////////////////////// // Structure // ============= TYPES ================ C(MolScript.structureQuery.type.elementSymbol, (ctx, v) => ElementSymbol(v[0](ctx))), C(MolScript.structureQuery.type.atomName, (ctx, v) => toUpperCase(v[0](ctx))), // TODO: // C(MolScript.structureQuery.type.bondFlags, (ctx, v) => StructureRuntime.BondProperties.createFlags(env, v)), // C(MolScript.structureQuery.type.secondaryStructureFlags, (ctx, v) => StructureRuntime.AtomProperties.createSecondaryStructureFlags(env, v)), // C(MolScript.structureQuery.type.entityType, (ctx, v) => StructureRuntime.Common.entityType(v[0](ctx))), // C(MolScript.structureQuery.type.ringFingerprint, (ctx, v) => StructureRuntime.Common.ringFingerprint(env, v as any)), // C(MolScript.structureQuery.type.authResidueId, (ctx, v) => ResidueIdentifier.auth(v[0](ctx), v[1](ctx), v[2] && v[2](ctx))), // C(MolScript.structureQuery.type.labelResidueId, (ctx, v) => ResidueIdentifier.label(v[0](ctx), v[1](ctx), v[2](ctx), v[3] && v[3](ctx))), // ============= SLOTS ================ // TODO: slots might not be needed after all: reducer simply pushes/pops current element C(MolScript.structureQuery.slot.element, (ctx, _) => ctx.element), // C(MolScript.structureQuery.slot.elementSetReduce, (ctx, _) => ctx.element), // ============= GENERATORS ================ D(MolScript.structureQuery.generator.atomGroups, (ctx, xs) => Queries.generators.atoms({ entityTest: xs['entity-test'], chainTest: xs['chain-test'], residueTest: xs['residue-test'], atomTest: xs['atom-test'], groupBy: xs['group-by'] })(ctx)), D(MolScript.structureQuery.generator.all, (ctx) => Queries.generators.all(ctx)), D(MolScript.structureQuery.generator.empty, (ctx) => Queries.generators.none(ctx)), // ============= MODIFIERS ================ D(MolScript.structureQuery.modifier.includeSurroundings, (ctx, xs) => Queries.modifiers.includeSurroundings(xs[0] as any, { radius: xs['radius'](ctx), wholeResidues: !!(xs['as-whole-residues'] && xs['as-whole-residues'](ctx)), elementRadius: xs['atom-radius'] })(ctx)), D(MolScript.structureQuery.modifier.wholeResidues, (ctx, xs) => Queries.modifiers.wholeResidues(xs[0] as any)(ctx)), D(MolScript.structureQuery.modifier.union, (ctx, xs) => Queries.modifiers.union(xs[0] as any)(ctx)), D(MolScript.structureQuery.modifier.expandProperty, (ctx, xs) => Queries.modifiers.expandProperty(xs[0] as any, xs['property'])(ctx)), D(MolScript.structureQuery.modifier.exceptBy, (ctx, xs) => Queries.modifiers.exceptBy(xs[0] as any, xs['by'] as any)(ctx)), // ============= COMBINATORS ================ D(MolScript.structureQuery.combinator.merge, (ctx, xs) => Queries.combinators.merge(xs as any)(ctx)), // ============= ATOM PROPERTIES ================ // ~~~ CORE ~~~ D(MolScript.structureQuery.atomProperty.core.elementSymbol, atomProp(StructureProperties.atom.type_symbol)), D(MolScript.structureQuery.atomProperty.core.vdw, (ctx, _) => VdwRadius(StructureProperties.atom.type_symbol(ctx.element))), D(MolScript.structureQuery.atomProperty.core.mass, (ctx, _) => AtomWeight(StructureProperties.atom.type_symbol(ctx.element))), D(MolScript.structureQuery.atomProperty.core.atomicNumber, (ctx, _) => AtomNumber(StructureProperties.atom.type_symbol(ctx.element))), D(MolScript.structureQuery.atomProperty.core.x, atomProp(StructureProperties.atom.x)), D(MolScript.structureQuery.atomProperty.core.y, atomProp(StructureProperties.atom.y)), D(MolScript.structureQuery.atomProperty.core.z, atomProp(StructureProperties.atom.z)), D(MolScript.structureQuery.atomProperty.core.sourceIndex, atomProp(StructureProperties.atom.sourceIndex)), D(MolScript.structureQuery.atomProperty.core.operatorName, atomProp(StructureProperties.unit.operator_name)), D(MolScript.structureQuery.atomProperty.core.atomKey, (ctx, _) => cantorPairing(ctx.element.unit.id, ctx.element.element)), // TODO: // D(MolScript.structureQuery.atomProperty.core.bondCount, (ctx, _) => ), // ~~~ TOPOLOGY ~~~ // TODO // ~~~ MACROMOLECULAR ~~~ // TODO: // // identifiers // labelResidueId: prop((env, v) => ResidueIdentifier.labelOfResidueIndex(env.context.model, getAddress(env, v).residue)), // authResidueId: prop((env, v) => ResidueIdentifier.authOfResidueIndex(env.context.model, getAddress(env, v).residue)), // keys D(MolScript.structureQuery.atomProperty.macromolecular.residueKey, (ctx, _) => StructureElement.residueIndex(ctx.element)), D(MolScript.structureQuery.atomProperty.macromolecular.chainKey, (ctx, _) => StructureElement.chainIndex(ctx.element)), D(MolScript.structureQuery.atomProperty.macromolecular.entityKey, (ctx, _) => StructureElement.entityIndex(ctx.element)), // mmCIF D(MolScript.structureQuery.atomProperty.macromolecular.id, atomProp(StructureProperties.atom.id)), D(MolScript.structureQuery.atomProperty.macromolecular.isHet, (ctx, _) => StructureProperties.residue.group_PDB(ctx.element) !== 'ATOM'), D(MolScript.structureQuery.atomProperty.macromolecular.label_atom_id, atomProp(StructureProperties.atom.label_atom_id)), D(MolScript.structureQuery.atomProperty.macromolecular.label_alt_id, atomProp(StructureProperties.atom.label_alt_id)), D(MolScript.structureQuery.atomProperty.macromolecular.label_asym_id, atomProp(StructureProperties.chain.label_asym_id)), D(MolScript.structureQuery.atomProperty.macromolecular.label_comp_id, atomProp(StructureProperties.residue.label_comp_id)), D(MolScript.structureQuery.atomProperty.macromolecular.label_seq_id, atomProp(StructureProperties.residue.label_seq_id)), D(MolScript.structureQuery.atomProperty.macromolecular.label_entity_id, atomProp(StructureProperties.entity.id)), D(MolScript.structureQuery.atomProperty.macromolecular.auth_atom_id, atomProp(StructureProperties.atom.auth_atom_id)), D(MolScript.structureQuery.atomProperty.macromolecular.auth_asym_id, atomProp(StructureProperties.chain.auth_asym_id)), D(MolScript.structureQuery.atomProperty.macromolecular.auth_comp_id, atomProp(StructureProperties.residue.auth_comp_id)), D(MolScript.structureQuery.atomProperty.macromolecular.auth_seq_id, atomProp(StructureProperties.residue.auth_seq_id)), D(MolScript.structureQuery.atomProperty.macromolecular.pdbx_PDB_ins_code, atomProp(StructureProperties.residue.pdbx_PDB_ins_code)), D(MolScript.structureQuery.atomProperty.macromolecular.pdbx_formal_charge, atomProp(StructureProperties.atom.pdbx_formal_charge)), D(MolScript.structureQuery.atomProperty.macromolecular.occupancy, atomProp(StructureProperties.atom.occupancy)), D(MolScript.structureQuery.atomProperty.macromolecular.B_iso_or_equiv, atomProp(StructureProperties.atom.B_iso_or_equiv)), D(MolScript.structureQuery.atomProperty.macromolecular.entityType, atomProp(StructureProperties.entity.type)), D(MolScript.structureQuery.atomProperty.macromolecular.isModified, (ctx, _) => ctx.element.unit.model.properties.modifiedResidues.parentId.has(StructureProperties.residue.label_comp_id(ctx.element))), D(MolScript.structureQuery.atomProperty.macromolecular.modifiedParentName, (ctx, _) => { const id = StructureProperties.residue.label_comp_id(ctx.element); return ctx.element.unit.model.properties.modifiedResidues.parentId.get(id) || id }) // TODO // MolScript.structureQuery.atomProperty.macromolecular.secondaryStructureKey // MolScript.structureQuery.atomProperty.macromolecular.secondaryStructureFlags // ============= BOND PROPERTIES ================ ]; function atomProp(p: (e: StructureElement) => any): (ctx: QueryContext, _: any) => any { return (ctx, _) => p(ctx.element); } (function () { for (const s of symbols) { DefaultQueryRuntimeTable.addSymbol(s); } })();