mol-script.ts 3.8 KB

1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980818283848586878889909192939495
  1. import { MolScriptBuilder } from '../mol-script/language/builder';
  2. import { compile, QuerySymbolRuntime, DefaultQueryRuntimeTable } from '../mol-script/runtime/query/compiler';
  3. import { QueryContext, Structure, StructureQuery, CustomPropertyDescriptor } from '../mol-model/structure';
  4. import { readCifFile, getModelsAndStructure } from '../apps/structure-info/model';
  5. import { CustomPropSymbol } from '../mol-script/language/symbol';
  6. import Type from '../mol-script/language/type';
  7. import { parseMolScript } from '../mol-script/language/parser';
  8. import * as util from 'util'
  9. import { transpileMolScript } from '../mol-script/script/mol-script/symbols';
  10. import { formatMolScript } from '../mol-script/language/expression-formatter';
  11. import { StructureQualityReport, StructureQualityReportProvider } from '../mol-model-props/pdbe/structure-quality-report';
  12. import fetch from 'node-fetch';
  13. // import Examples from 'mol-script/script/mol-script/examples'
  14. // import { parseMolScript } from 'mol-script/script/mol-script/parser'
  15. // //import { compileAST } from 'mol-script/script/mol-script/compile';
  16. // for (const e of Examples) {
  17. // const expr = parseMolScript(e.value)[0];
  18. // console.log(e.name, util.inspect(expr, true, 10, true));
  19. // }
  20. // const exprs = parseMolScript(`(sel.atom.atom-groups
  21. // :residue-test (= atom.auth_comp_id ALA)
  22. // ;; ho ho ho
  23. // :atom-test (set.has { _C _N } atom.el)) ; comm
  24. // ;; this is a comment
  25. // ((hi) (ho))`);
  26. // ;; :residue-test (= atom.label_comp_id REA)
  27. const exprs = parseMolScript(`(sel.atom.atom-groups
  28. :residue-test (> pdbe.structure-quality.issue-count 0)
  29. :atom-test (= atom.el _C))`);
  30. const tsp = transpileMolScript(exprs[0]);
  31. // console.log(util.inspect(exprs, true, 10, true));
  32. console.log(util.inspect(tsp, true, 10, true));
  33. console.log(formatMolScript);
  34. console.log(formatMolScript(tsp));
  35. // //console.log(expr);
  36. const expr = MolScriptBuilder.core.math.add([1, 2, 3]);
  37. const compiled = compile<number>(expr);
  38. const result = compiled(new QueryContext(Structure.Empty));
  39. console.log(result);
  40. const CustomProp = CustomPropertyDescriptor({
  41. name: 'test_prop',
  42. cifExport: { prefix: '', categories: [ ]},
  43. symbols: {
  44. residueIndex: QuerySymbolRuntime.Dynamic(CustomPropSymbol('custom.test-prop', 'residue-index', Type.Num), ctx => {
  45. const e = ctx.element;
  46. // console.log(e.element, e.unit.model.atomicHierarchy.residueAtomSegments.index[e.element])
  47. return e.unit.model.atomicHierarchy.residueAtomSegments.index[e.element];
  48. })
  49. }
  50. });
  51. DefaultQueryRuntimeTable.addCustomProp(CustomProp);
  52. DefaultQueryRuntimeTable.addCustomProp(StructureQualityReportProvider.descriptor);
  53. export async function testQ() {
  54. const frame = await readCifFile('e:/test/quick/1cbs_updated.cif');
  55. const { structure } = await getModelsAndStructure(frame);
  56. const model = structure.models[0]
  57. const rawData = await fetch(`https://www.ebi.ac.uk/pdbe/api/validation/residuewise_outlier_summary/entry/${model.entryId.toLowerCase()}`, { timeout: 1500 });
  58. const data = StructureQualityReport.fromJson(model, await rawData.json());
  59. StructureQualityReportProvider.set(model, { serverUrl: '' }, data)
  60. let expr = MolScriptBuilder.struct.generator.atomGroups({
  61. 'atom-test': MolScriptBuilder.core.rel.eq([
  62. MolScriptBuilder.struct.atomProperty.core.elementSymbol(),
  63. MolScriptBuilder.es('C')
  64. ]),
  65. // 'residue-test': MolScriptBuilder.core.rel.eq([
  66. // MolScriptBuilder.struct.atomProperty.macromolecular.label_comp_id(),
  67. // 'REA'
  68. // ])
  69. 'residue-test': MolScriptBuilder.core.rel.inRange([CustomProp.symbols.residueIndex.symbol(), 1, 5])
  70. });
  71. expr = tsp;
  72. const compiled = compile<StructureQuery>(expr);
  73. const result = compiled(new QueryContext(structure));
  74. console.log(result);
  75. }
  76. testQ();