nmr-processing 9.3.7 → 9.3.8-pre.1681135907

This diff represents the content of publicly available package versions that have been released to one of the supported registries. The information contained in this diff is provided for informational purposes only and reflects changes between package versions as they appear in their respective public registries.
Files changed (84) hide show
  1. package/lib/apodization/apodization.d.ts +9 -4
  2. package/lib/apodization/apodization.js +3 -1
  3. package/lib/apodization/apodization.js.map +1 -1
  4. package/lib/apodization/applyWindow.js.map +1 -1
  5. package/lib/apodization/compose.js.map +1 -1
  6. package/lib/apodization/getFunction.js +3 -0
  7. package/lib/apodization/getFunction.js.map +1 -1
  8. package/lib/apodization/shapes/WindowFunctions.d.ts +6 -1
  9. package/lib/apodization/shapes/sineBell.d.ts +27 -0
  10. package/lib/apodization/shapes/sineBell.js +17 -0
  11. package/lib/apodization/shapes/sineBell.js.map +1 -0
  12. package/lib/apodization/utils/getData.js.map +1 -1
  13. package/lib/signals/signalsToFID.d.ts +3 -0
  14. package/lib/utilities/checkData2DFid.d.ts +3 -0
  15. package/lib/utilities/checkData2DFid.js +10 -0
  16. package/lib/utilities/checkData2DFid.js.map +1 -0
  17. package/lib/xyz/Data2D.d.ts +21 -0
  18. package/lib/xyz/Data2D.js +3 -0
  19. package/lib/xyz/Data2D.js.map +1 -0
  20. package/lib/xyz/quadrature.d.ts +8 -0
  21. package/lib/xyz/quadrature.js +54 -0
  22. package/lib/xyz/quadrature.js.map +1 -0
  23. package/lib/xyz/util/fft2d/digitalFilter.d.ts +7 -0
  24. package/lib/xyz/util/fft2d/digitalFilter.js +20 -0
  25. package/lib/xyz/util/fft2d/digitalFilter.js.map +1 -0
  26. package/lib/xyz/util/fft2d/fftDirectDimension.d.ts +24 -0
  27. package/lib/xyz/util/fft2d/fftDirectDimension.js +57 -0
  28. package/lib/xyz/util/fft2d/fftDirectDimension.js.map +1 -0
  29. package/lib/xyz/util/fft2d/fftIndirectDimension.d.ts +23 -0
  30. package/lib/xyz/util/fft2d/fftIndirectDimension.js +72 -0
  31. package/lib/xyz/util/fft2d/fftIndirectDimension.js.map +1 -0
  32. package/lib/xyz/util/fft2d/removeDCOffset.d.ts +5 -0
  33. package/lib/xyz/util/fft2d/removeDCOffset.js +19 -0
  34. package/lib/xyz/util/fft2d/removeDCOffset.js.map +1 -0
  35. package/lib/xyz/util/fft2d/zeroFilling.d.ts +10 -0
  36. package/lib/xyz/util/fft2d/zeroFilling.js +29 -0
  37. package/lib/xyz/util/fft2d/zeroFilling.js.map +1 -0
  38. package/lib/xyz/xyzBidimensionalFFT.d.ts +41 -0
  39. package/lib/xyz/xyzBidimensionalFFT.js +79 -0
  40. package/lib/xyz/xyzBidimensionalFFT.js.map +1 -0
  41. package/lib-esm/apodization/apodization.js +3 -1
  42. package/lib-esm/apodization/apodization.js.map +1 -1
  43. package/lib-esm/apodization/applyWindow.js.map +1 -1
  44. package/lib-esm/apodization/compose.js.map +1 -1
  45. package/lib-esm/apodization/getFunction.js +3 -0
  46. package/lib-esm/apodization/getFunction.js.map +1 -1
  47. package/lib-esm/apodization/shapes/sineBell.js +13 -0
  48. package/lib-esm/apodization/shapes/sineBell.js.map +1 -0
  49. package/lib-esm/apodization/utils/getData.js.map +1 -1
  50. package/lib-esm/utilities/checkData2DFid.js +6 -0
  51. package/lib-esm/utilities/checkData2DFid.js.map +1 -0
  52. package/lib-esm/xyz/Data2D.js +2 -0
  53. package/lib-esm/xyz/Data2D.js.map +1 -0
  54. package/lib-esm/xyz/quadrature.js +50 -0
  55. package/lib-esm/xyz/quadrature.js.map +1 -0
  56. package/lib-esm/xyz/util/fft2d/digitalFilter.js +16 -0
  57. package/lib-esm/xyz/util/fft2d/digitalFilter.js.map +1 -0
  58. package/lib-esm/xyz/util/fft2d/fftDirectDimension.js +50 -0
  59. package/lib-esm/xyz/util/fft2d/fftDirectDimension.js.map +1 -0
  60. package/lib-esm/xyz/util/fft2d/fftIndirectDimension.js +65 -0
  61. package/lib-esm/xyz/util/fft2d/fftIndirectDimension.js.map +1 -0
  62. package/lib-esm/xyz/util/fft2d/removeDCOffset.js +15 -0
  63. package/lib-esm/xyz/util/fft2d/removeDCOffset.js.map +1 -0
  64. package/lib-esm/xyz/util/fft2d/zeroFilling.js +25 -0
  65. package/lib-esm/xyz/util/fft2d/zeroFilling.js.map +1 -0
  66. package/lib-esm/xyz/xyzBidimensionalFFT.js +75 -0
  67. package/lib-esm/xyz/xyzBidimensionalFFT.js.map +1 -0
  68. package/package.json +17 -17
  69. package/src/apodization/apodization.ts +6 -4
  70. package/src/apodization/applyWindow.ts +0 -1
  71. package/src/apodization/compose.ts +1 -0
  72. package/src/apodization/getFunction.ts +3 -0
  73. package/src/apodization/shapes/WindowFunctions.ts +7 -1
  74. package/src/apodization/shapes/sineBell.ts +41 -0
  75. package/src/apodization/utils/getData.ts +0 -1
  76. package/src/utilities/checkData2DFid.ts +10 -0
  77. package/src/xyz/Data2D.ts +25 -0
  78. package/src/xyz/quadrature.ts +69 -0
  79. package/src/xyz/util/fft2d/digitalFilter.ts +24 -0
  80. package/src/xyz/util/fft2d/fftDirectDimension.ts +111 -0
  81. package/src/xyz/util/fft2d/fftIndirectDimension.ts +123 -0
  82. package/src/xyz/util/fft2d/removeDCOffset.ts +19 -0
  83. package/src/xyz/util/fft2d/zeroFilling.ts +39 -0
  84. package/src/xyz/xyzBidimensionalFFT.ts +123 -0
@@ -6,17 +6,22 @@ export interface DataReIm {
6
6
  /** Array of im values */
7
7
  im: DoubleArray;
8
8
  }
9
- interface ApodizationOptions {
9
+ export interface ApodizationOptions {
10
10
  /**
11
11
  * number of points at the end of the shape to apply the window function inverted
12
12
  * @default 0
13
13
  */
14
14
  pointsToShift?: number;
15
- compose: ShapeOptions;
15
+ compose?: ShapeOptions;
16
16
  }
17
- export declare function apodization(data: DataReIm, options: ApodizationOptions): {
17
+ export declare function apodization(data: DataReIm, options?: ApodizationOptions): {
18
+ windowData: never[];
19
+ /** Array of re values */
20
+ re: DoubleArray;
21
+ /** Array of im values */
22
+ im: DoubleArray;
23
+ } | {
18
24
  re: DoubleArray;
19
25
  im: DoubleArray;
20
26
  windowData: Float64Array;
21
27
  };
22
- export {};
@@ -3,8 +3,10 @@ Object.defineProperty(exports, "__esModule", { value: true });
3
3
  exports.apodization = void 0;
4
4
  const applyWindow_1 = require("./applyWindow");
5
5
  const compose_1 = require("./compose");
6
- function apodization(data, options) {
6
+ function apodization(data, options = {}) {
7
7
  const { compose: composeOptions, pointsToShift } = options;
8
+ if (!composeOptions)
9
+ return { ...data, windowData: [] };
8
10
  const windowData = (0, compose_1.compose)(composeOptions);
9
11
  const applyWindowOptions = {
10
12
  windowData,
@@ -1 +1 @@
1
- {"version":3,"file":"apodization.js","sourceRoot":"","sources":["../../src/apodization/apodization.ts"],"names":[],"mappings":";;;AAEA,+CAA4C;AAC5C,uCAAkD;AAiBlD,SAAgB,WAAW,CAAC,IAAc,EAAE,OAA2B;IACrE,MAAM,EAAE,OAAO,EAAE,cAAc,EAAE,aAAa,EAAE,GAAG,OAAO,CAAC;IAC3D,MAAM,UAAU,GAAG,IAAA,iBAAO,EAAC,cAAc,CAAC,CAAC;IAE3C,MAAM,kBAAkB,GAAG;QACzB,UAAU;QACV,aAAa;KACd,CAAC;IAEF,MAAM,EAAE,GAAG,IAAA,yBAAW,EAAC,IAAI,CAAC,EAAE,EAAE,kBAAkB,CAAC,CAAC;IACpD,MAAM,EAAE,GAAG,IAAA,yBAAW,EAAC,IAAI,CAAC,EAAE,EAAE,kBAAkB,CAAC,CAAC;IAEpD,OAAO,EAAE,EAAE,EAAE,EAAE,EAAE,UAAU,EAAE,CAAC;AAChC,CAAC;AAbD,kCAaC"}
1
+ {"version":3,"file":"apodization.js","sourceRoot":"","sources":["../../src/apodization/apodization.ts"],"names":[],"mappings":";;;AAEA,+CAA4C;AAC5C,uCAAkD;AAiBlD,SAAgB,WAAW,CAAC,IAAc,EAAE,UAA8B,EAAE;IAC1E,MAAM,EAAE,OAAO,EAAE,cAAc,EAAE,aAAa,EAAE,GAAG,OAAO,CAAC;IAE3D,IAAI,CAAC,cAAc;QAAE,OAAO,EAAE,GAAG,IAAI,EAAE,UAAU,EAAE,EAAE,EAAE,CAAC;IAExD,MAAM,UAAU,GAAG,IAAA,iBAAO,EAAC,cAAc,CAAC,CAAC;IAE3C,MAAM,kBAAkB,GAAG;QACzB,UAAU;QACV,aAAa;KACd,CAAC;IACF,MAAM,EAAE,GAAG,IAAA,yBAAW,EAAC,IAAI,CAAC,EAAE,EAAE,kBAAkB,CAAC,CAAC;IACpD,MAAM,EAAE,GAAG,IAAA,yBAAW,EAAC,IAAI,CAAC,EAAE,EAAE,kBAAkB,CAAC,CAAC;IAEpD,OAAO,EAAE,EAAE,EAAE,EAAE,EAAE,UAAU,EAAE,CAAC;AAChC,CAAC;AAfD,kCAeC"}
@@ -1 +1 @@
1
- {"version":3,"file":"applyWindow.js","sourceRoot":"","sources":["../../src/apodization/applyWindow.ts"],"names":[],"mappings":";;;AA0BA;;GAEG;AACH,SAAgB,WAAW,CAAC,IAAiB,EAAE,OAA2B;IACxE,MAAM,UAAU,GAAG,IAAI,CAAC,MAAM,CAAC;IAE/B,MAAM,EACJ,UAAU,EACV,KAAK,GAAG,CAAC,EACT,MAAM,GAAG,UAAU,EACnB,aAAa,GAAG,CAAC,EACjB,MAAM,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,GAChC,GAAG,OAAO,CAAC;IAEZ,MAAM,aAAa,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,MAAM,EAAE,UAAU,GAAG,aAAa,CAAC,CAAC;IAC3E,KAAK,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,aAAa,EAAE,CAAC,EAAE,EAAE;QACjD,MAAM,CAAC,CAAC,CAAC,IAAI,UAAU,CAAC,CAAC,EAAE,CAAC,CAAC;KAC9B;IAED,KAAK,IAAI,CAAC,GAAG,UAAU,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,UAAU,GAAG,aAAa,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;QAC3E,MAAM,CAAC,CAAC,CAAC,IAAI,UAAU,CAAC,CAAC,EAAE,CAAC,CAAC;KAC9B;IAED,OAAO,MAAM,CAAC;AAChB,CAAC;AArBD,kCAqBC"}
1
+ {"version":3,"file":"applyWindow.js","sourceRoot":"","sources":["../../src/apodization/applyWindow.ts"],"names":[],"mappings":";;;AA0BA;;GAEG;AACH,SAAgB,WAAW,CAAC,IAAiB,EAAE,OAA2B;IACxE,MAAM,UAAU,GAAG,IAAI,CAAC,MAAM,CAAC;IAE/B,MAAM,EACJ,UAAU,EACV,KAAK,GAAG,CAAC,EACT,MAAM,GAAG,UAAU,EACnB,aAAa,GAAG,CAAC,EACjB,MAAM,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,GAChC,GAAG,OAAO,CAAC;IAEZ,MAAM,aAAa,GAAG,IAAI,CAAC,GAAG,CAAC,KAAK,GAAG,MAAM,EAAE,UAAU,GAAG,aAAa,CAAC,CAAC;IAC3E,KAAK,IAAI,CAAC,GAAG,KAAK,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,aAAa,EAAE,CAAC,EAAE,EAAE;QACjD,MAAM,CAAC,CAAC,CAAC,IAAI,UAAU,CAAC,CAAC,EAAE,CAAC,CAAC;KAC9B;IACD,KAAK,IAAI,CAAC,GAAG,UAAU,GAAG,CAAC,EAAE,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,UAAU,GAAG,aAAa,GAAG,CAAC,EAAE,CAAC,EAAE,EAAE;QAC3E,MAAM,CAAC,CAAC,CAAC,IAAI,UAAU,CAAC,CAAC,EAAE,CAAC,CAAC;KAC9B;IAED,OAAO,MAAM,CAAC;AAChB,CAAC;AApBD,kCAoBC"}
@@ -1 +1 @@
1
- {"version":3,"file":"compose.js","sourceRoot":"","sources":["../../src/apodization/compose.ts"],"names":[],"mappings":";;;AAAA,+CAA4C;AAC5C,+CAA4C;AAE5C,6CAA0C;AAyB1C,SAAgB,OAAO,CAAC,YAA0B;IAChD,MAAM,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,GAAG,EAAE,EAAE,GAAG,YAAY,CAAC;IACzD,IAAI,IAAI,GAAG,IAAI,YAAY,CAAC,UAAU,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAChD,MAAM,CAAC,OAAO,CAAC,CAAC,OAAO,EAAE,EAAE;QACzB,MAAM,EAAE,KAAK,EAAE,KAAK,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC;QACrC,MAAM,EAAE,OAAO,EAAE,YAAY,EAAE,GAAG,KAAK,CAAC;QACxC,MAAM,MAAM,GACV,QAAQ,IAAI,YAAY,CAAC,CAAC,CAAC,YAAY,CAAC,MAAM,CAAC,CAAC,CAAC,UAAU,GAAG,KAAK,CAAC;QAEtE,MAAM,UAAU,GAAG,IAAA,iBAAO,EAAC,IAAA,yBAAW,EAAC,KAAK,CAAC,EAAE,MAAM,CAAC,CAAC;QACvD,IAAA,yBAAW,EAAC,IAAI,EAAE;YAChB,UAAU;YACV,MAAM;YACN,KAAK;YACL,MAAM,EAAE,IAAI;SACb,CAAC,CAAC;IACL,CAAC,CAAC,CAAC;IACH,OAAO,IAAI,CAAC;AACd,CAAC;AAlBD,0BAkBC"}
1
+ {"version":3,"file":"compose.js","sourceRoot":"","sources":["../../src/apodization/compose.ts"],"names":[],"mappings":";;;AAAA,+CAA4C;AAC5C,+CAA4C;AAE5C,6CAA0C;AAyB1C,SAAgB,OAAO,CAAC,YAA0B;IAChD,MAAM,EAAE,MAAM,EAAE,UAAU,EAAE,MAAM,GAAG,EAAE,EAAE,GAAG,YAAY,CAAC;IACzD,IAAI,IAAI,GAAG,IAAI,YAAY,CAAC,UAAU,CAAC,CAAC,IAAI,CAAC,CAAC,CAAC,CAAC;IAChD,MAAM,CAAC,OAAO,CAAC,CAAC,OAAO,EAAE,EAAE;QACzB,MAAM,EAAE,KAAK,EAAE,KAAK,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC;QACrC,MAAM,EAAE,OAAO,EAAE,YAAY,EAAE,GAAG,KAAK,CAAC;QACxC,MAAM,MAAM,GACV,QAAQ,IAAI,YAAY,CAAC,CAAC,CAAC,YAAY,CAAC,MAAM,CAAC,CAAC,CAAC,UAAU,GAAG,KAAK,CAAC;QAEtE,MAAM,UAAU,GAAG,IAAA,iBAAO,EAAC,IAAA,yBAAW,EAAC,KAAK,CAAC,EAAE,MAAM,CAAC,CAAC;QAEvD,IAAA,yBAAW,EAAC,IAAI,EAAE;YAChB,UAAU;YACV,MAAM;YACN,KAAK;YACL,MAAM,EAAE,IAAI;SACb,CAAC,CAAC;IACL,CAAC,CAAC,CAAC;IACH,OAAO,IAAI,CAAC;AACd,CAAC;AAnBD,0BAmBC"}
@@ -3,6 +3,7 @@ Object.defineProperty(exports, "__esModule", { value: true });
3
3
  exports.getFunction = void 0;
4
4
  const exponential_1 = require("./shapes/exponential");
5
5
  const lorentzToGauss_1 = require("./shapes/lorentzToGauss");
6
+ const sineBell_1 = require("./shapes/sineBell");
6
7
  function getFunction(shape) {
7
8
  const { kind, options } = shape;
8
9
  switch (kind) {
@@ -10,6 +11,8 @@ function getFunction(shape) {
10
11
  return (0, exponential_1.exponential)(options);
11
12
  case 'lorentzToGauss':
12
13
  return (0, lorentzToGauss_1.lorentzToGauss)(options);
14
+ case 'sineBell':
15
+ return (0, sineBell_1.sineBell)(options);
13
16
  default:
14
17
  throw Error(`Unknown distribution ${kind}`);
15
18
  }
@@ -1 +1 @@
1
- {"version":3,"file":"getFunction.js","sourceRoot":"","sources":["../../src/apodization/getFunction.ts"],"names":[],"mappings":";;;AACA,sDAAmD;AACnD,4DAAyD;AAEzD,SAAgB,WAAW,CAAC,KAAsB;IAChD,MAAM,EAAE,IAAI,EAAE,OAAO,EAAE,GAAG,KAAK,CAAC;IAChC,QAAQ,IAAI,EAAE;QACZ,KAAK,aAAa;YAChB,OAAO,IAAA,yBAAW,EAAC,OAAO,CAAC,CAAC;QAC9B,KAAK,gBAAgB;YACnB,OAAO,IAAA,+BAAc,EAAC,OAAO,CAAC,CAAC;QACjC;YACE,MAAM,KAAK,CAAC,wBAAwB,IAAc,EAAE,CAAC,CAAC;KACzD;AACH,CAAC;AAVD,kCAUC"}
1
+ {"version":3,"file":"getFunction.js","sourceRoot":"","sources":["../../src/apodization/getFunction.ts"],"names":[],"mappings":";;;AACA,sDAAmD;AACnD,4DAAyD;AACzD,gDAA6C;AAE7C,SAAgB,WAAW,CAAC,KAAsB;IAChD,MAAM,EAAE,IAAI,EAAE,OAAO,EAAE,GAAG,KAAK,CAAC;IAChC,QAAQ,IAAI,EAAE;QACZ,KAAK,aAAa;YAChB,OAAO,IAAA,yBAAW,EAAC,OAAO,CAAC,CAAC;QAC9B,KAAK,gBAAgB;YACnB,OAAO,IAAA,+BAAc,EAAC,OAAO,CAAC,CAAC;QACjC,KAAK,UAAU;YACb,OAAO,IAAA,mBAAQ,EAAC,OAAO,CAAC,CAAC;QAC3B;YACE,MAAM,KAAK,CAAC,wBAAwB,IAAc,EAAE,CAAC,CAAC;KACzD;AACH,CAAC;AAZD,kCAYC"}
@@ -1,5 +1,6 @@
1
1
  import type { ExponentialOptions } from './exponential';
2
2
  import type { LorentzToGaussOptions } from './lorentzToGauss';
3
+ import { SineBellOptions } from './sineBell';
3
4
  interface Exponential {
4
5
  kind: 'exponential';
5
6
  options: ExponentialOptions;
@@ -8,5 +9,9 @@ interface LorentToGauss {
8
9
  kind: 'lorentzToGauss';
9
10
  options: LorentzToGaussOptions;
10
11
  }
11
- export type WindowFunctions = Exponential | LorentToGauss;
12
+ interface SineBell {
13
+ kind: 'sineBell';
14
+ options: SineBellOptions;
15
+ }
16
+ export type WindowFunctions = Exponential | LorentToGauss | SineBell;
12
17
  export {};
@@ -0,0 +1,27 @@
1
+ export interface SineBellOptions {
2
+ /**
3
+ * Specifies the exponent of the sine-bell; Non-integer values are allowed.
4
+ * Common values are 1.0 (for ordinary sine-bell) and 2.0 (for squared-bell functions).
5
+ * @default 1
6
+ */
7
+ exponent?: number;
8
+ /**
9
+ * Specifies the starting point of the sine-bell in units of pi radians.
10
+ * Common values are 0.0 (for a sine window which starts height at 0.0) and
11
+ * 0.5 (for a cosine window, which starts at height 1.0).
12
+ * @default 0
13
+ */
14
+ offset?: number;
15
+ /**
16
+ * Specifies the ending point of the sine-bell in units of pi radians.
17
+ * Common values are 1.0 (for a window which goes to 0.0 height at the last point) and
18
+ * 0.95 (for a window which doesn't go all the way to 0.0).
19
+ * @default 1
20
+ */
21
+ end?: number;
22
+ /**
23
+ * Specifies the number of points in the window function.
24
+ */
25
+ length: number;
26
+ }
27
+ export declare function sineBell(options: SineBellOptions): (i: number) => number;
@@ -0,0 +1,17 @@
1
+ "use strict";
2
+ Object.defineProperty(exports, "__esModule", { value: true });
3
+ exports.sineBell = void 0;
4
+ function sineBell(options) {
5
+ const { exponent = 1, offset = 0, end = 1, length } = options;
6
+ if (length === undefined) {
7
+ throw new Error('length options should exists for sineBell shape');
8
+ }
9
+ if (end < 0 || end > 1) {
10
+ throw new Error(`the end parameter should be [0-1]`);
11
+ }
12
+ const c1 = Math.PI * offset;
13
+ const c2 = (Math.PI * (end - offset)) / (length - 1);
14
+ return (i) => Math.pow(Math.sin(c1 + c2 * i), exponent);
15
+ }
16
+ exports.sineBell = sineBell;
17
+ //# sourceMappingURL=sineBell.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"sineBell.js","sourceRoot":"","sources":["../../../src/apodization/shapes/sineBell.ts"],"names":[],"mappings":";;;AA2BA,SAAgB,QAAQ,CAAC,OAAwB;IAC/C,MAAM,EAAE,QAAQ,GAAG,CAAC,EAAE,MAAM,GAAG,CAAC,EAAE,GAAG,GAAG,CAAC,EAAE,MAAM,EAAE,GAAG,OAAO,CAAC;IAE9D,IAAI,MAAM,KAAK,SAAS,EAAE;QACxB,MAAM,IAAI,KAAK,CAAC,iDAAiD,CAAC,CAAC;KACpE;IAED,IAAI,GAAG,GAAG,CAAC,IAAI,GAAG,GAAG,CAAC,EAAE;QACtB,MAAM,IAAI,KAAK,CAAC,mCAAmC,CAAC,CAAC;KACtD;IACD,MAAM,EAAE,GAAG,IAAI,CAAC,EAAE,GAAG,MAAM,CAAC;IAC5B,MAAM,EAAE,GAAG,CAAC,IAAI,CAAC,EAAE,GAAG,CAAC,GAAG,GAAG,MAAM,CAAC,CAAC,GAAG,CAAC,MAAM,GAAG,CAAC,CAAC,CAAC;IACrD,OAAO,CAAC,CAAS,EAAE,EAAE,CAAC,IAAI,CAAC,GAAG,CAAC,IAAI,CAAC,GAAG,CAAC,EAAE,GAAG,EAAE,GAAG,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC;AAClE,CAAC;AAbD,4BAaC"}
@@ -1 +1 @@
1
- {"version":3,"file":"getData.js","sourceRoot":"","sources":["../../../src/apodization/utils/getData.ts"],"names":[],"mappings":";;;AAAA;;GAEG;AACH,SAAgB,OAAO,CAAC,IAA2B,EAAE,MAAc;IACjE,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,CAAC;IACtC,IAAI,GAAG,GAAG,MAAM,CAAC,gBAAgB,CAAC;IAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE;QAC/B,MAAM,KAAK,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;QACtB,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;QAChB,IAAI,KAAK,GAAG,GAAG;YAAE,GAAG,GAAG,KAAK,CAAC;KAC9B;IAED,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE;QAAE,IAAI,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC;IAChD,OAAO,IAAI,CAAC;AACd,CAAC;AAXD,0BAWC"}
1
+ {"version":3,"file":"getData.js","sourceRoot":"","sources":["../../../src/apodization/utils/getData.ts"],"names":[],"mappings":";;;AAAA;;GAEG;AACH,SAAgB,OAAO,CAAC,IAA2B,EAAE,MAAc;IACjE,MAAM,IAAI,GAAG,IAAI,YAAY,CAAC,MAAM,CAAC,CAAC;IACtC,IAAI,GAAG,GAAG,MAAM,CAAC,gBAAgB,CAAC;IAClC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE;QAC/B,MAAM,KAAK,GAAG,IAAI,CAAC,CAAC,CAAC,CAAC;QACtB,IAAI,CAAC,CAAC,CAAC,GAAG,KAAK,CAAC;QAChB,IAAI,KAAK,GAAG,GAAG;YAAE,GAAG,GAAG,KAAK,CAAC;KAC9B;IACD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE;QAAE,IAAI,CAAC,CAAC,CAAC,IAAI,GAAG,CAAC;IAChD,OAAO,IAAI,CAAC;AACd,CAAC;AAVD,0BAUC"}
@@ -39,4 +39,7 @@ export interface SignalsToFIDOptions {
39
39
  export declare function signalsToFID(signals: NMRSignal1D[], options: SignalsToFIDOptions): {
40
40
  re: import("cheminfo-types").DoubleArray;
41
41
  im: import("cheminfo-types").DoubleArray;
42
+ } | {
43
+ re: import("cheminfo-types").DoubleArray;
44
+ im: import("cheminfo-types").DoubleArray;
42
45
  };
@@ -0,0 +1,3 @@
1
+ import { Data2D, Data2DFid } from '../xyz/Data2D';
2
+ export type Data2DFidComplete = Required<Data2DFid>;
3
+ export declare function checkData2DFid(data: Data2D): asserts data is Data2DFidComplete;
@@ -0,0 +1,10 @@
1
+ "use strict";
2
+ Object.defineProperty(exports, "__esModule", { value: true });
3
+ exports.checkData2DFid = void 0;
4
+ function checkData2DFid(data) {
5
+ if (!('re' in data && 'im' in data)) {
6
+ throw new Error('imaginary (im) data should exists');
7
+ }
8
+ }
9
+ exports.checkData2DFid = checkData2DFid;
10
+ //# sourceMappingURL=checkData2DFid.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"checkData2DFid.js","sourceRoot":"","sources":["../../src/utilities/checkData2DFid.ts"],"names":[],"mappings":";;;AAGA,SAAgB,cAAc,CAC5B,IAAY;IAEZ,IAAI,CAAC,CAAC,IAAI,IAAI,IAAI,IAAI,IAAI,IAAI,IAAI,CAAC,EAAE;QACnC,MAAM,IAAI,KAAK,CAAC,mCAAmC,CAAC,CAAC;KACtD;AACH,CAAC;AAND,wCAMC"}
@@ -0,0 +1,21 @@
1
+ import { DoubleArray } from 'cheminfo-types';
2
+ export interface MinMaxContent {
3
+ z: DoubleArray[];
4
+ minX: number;
5
+ minY: number;
6
+ minZ: number;
7
+ maxX: number;
8
+ maxY: number;
9
+ maxZ: number;
10
+ }
11
+ export type Data2D = Data2DFid | Data2DFt;
12
+ export interface Data2DFid {
13
+ re: MinMaxContent;
14
+ im?: MinMaxContent;
15
+ }
16
+ export interface Data2DFt {
17
+ rr: MinMaxContent;
18
+ ri?: MinMaxContent;
19
+ ir?: MinMaxContent;
20
+ ii?: MinMaxContent;
21
+ }
@@ -0,0 +1,3 @@
1
+ "use strict";
2
+ Object.defineProperty(exports, "__esModule", { value: true });
3
+ //# sourceMappingURL=Data2D.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"Data2D.js","sourceRoot":"","sources":["../../src/xyz/Data2D.ts"],"names":[],"mappings":""}
@@ -0,0 +1,8 @@
1
+ import { DoubleArray } from 'cheminfo-types';
2
+ import { Data2D } from './Data2D';
3
+ export declare function quadrature(data: Data2D, options: {
4
+ fnMode: string;
5
+ }): {
6
+ re: DoubleArray[];
7
+ im: DoubleArray[];
8
+ };
@@ -0,0 +1,54 @@
1
+ "use strict";
2
+ Object.defineProperty(exports, "__esModule", { value: true });
3
+ exports.quadrature = void 0;
4
+ const ml_spectra_processing_1 = require("ml-spectra-processing");
5
+ const checkData2DFid_1 = require("../utilities/checkData2DFid");
6
+ function quadrature(data, options) {
7
+ (0, checkData2DFid_1.checkData2DFid)(data);
8
+ const { fnMode } = options;
9
+ switch (fnMode) {
10
+ case 'ttpi':
11
+ return ttpiQuadrature(data);
12
+ case 'states':
13
+ return shrQuadrature(data);
14
+ case 'echo':
15
+ return echoQuadrature(data);
16
+ default:
17
+ return {
18
+ re: data.re.z,
19
+ im: data.im.z,
20
+ };
21
+ }
22
+ }
23
+ exports.quadrature = quadrature;
24
+ function echoQuadrature(data) {
25
+ const { re: { z: reMatrix }, im: { z: imMatrix }, } = data;
26
+ const newRe = [];
27
+ const newIm = [];
28
+ for (let i = 0; i < reMatrix.length; i += 2) {
29
+ let reCosine = reMatrix[i];
30
+ let imCosine = imMatrix[i];
31
+ let reSine = reMatrix[i + 1];
32
+ let imSine = imMatrix[i + 1];
33
+ newRe.push((0, ml_spectra_processing_1.xSubtract)((0, ml_spectra_processing_1.xMultiply)(imSine, -1), imCosine));
34
+ newRe.push((0, ml_spectra_processing_1.xSubtract)(reSine, reCosine));
35
+ newIm.push((0, ml_spectra_processing_1.xAdd)(reSine, reCosine));
36
+ newIm.push((0, ml_spectra_processing_1.xSubtract)(imSine, imCosine));
37
+ }
38
+ return { re: newRe, im: newIm };
39
+ }
40
+ function shrQuadrature(data) {
41
+ const { re, im } = data;
42
+ return { re: re.z, im: im.z };
43
+ }
44
+ function ttpiQuadrature(data) {
45
+ const newRe = [];
46
+ const newIm = [];
47
+ const { re: { z: reMatrix }, im: { z: imMatrix }, } = data;
48
+ for (let i = 0; i < reMatrix.length; i++) {
49
+ newRe.push(i % 2 ? reMatrix[i].slice() : (0, ml_spectra_processing_1.xMultiply)(reMatrix[i], -1));
50
+ newIm.push(i % 2 ? imMatrix[i].slice() : (0, ml_spectra_processing_1.xMultiply)(imMatrix[i], -1));
51
+ }
52
+ return { re: newRe, im: newIm };
53
+ }
54
+ //# sourceMappingURL=quadrature.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"quadrature.js","sourceRoot":"","sources":["../../src/xyz/quadrature.ts"],"names":[],"mappings":";;;AACA,iEAAmE;AAEnE,gEAAgF;AAGhF,SAAgB,UAAU,CAAC,IAAY,EAAE,OAA2B;IAClE,IAAA,+BAAc,EAAC,IAAI,CAAC,CAAC;IAErB,MAAM,EAAE,MAAM,EAAE,GAAG,OAAO,CAAC;IAE3B,QAAQ,MAAM,EAAE;QACd,KAAK,MAAM;YACT,OAAO,cAAc,CAAC,IAAI,CAAC,CAAC;QAC9B,KAAK,QAAQ;YACX,OAAO,aAAa,CAAC,IAAI,CAAC,CAAC;QAC7B,KAAK,MAAM;YACT,OAAO,cAAc,CAAC,IAAI,CAAC,CAAC;QAC9B;YACE,OAAO;gBACL,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;gBACb,EAAE,EAAE,IAAI,CAAC,EAAE,CAAC,CAAC;aACd,CAAC;KACL;AACH,CAAC;AAlBD,gCAkBC;AAED,SAAS,cAAc,CAAC,IAAuB;IAC7C,MAAM,EACJ,EAAE,EAAE,EAAE,CAAC,EAAE,QAAQ,EAAE,EACnB,EAAE,EAAE,EAAE,CAAC,EAAE,QAAQ,EAAE,GACpB,GAAG,IAAI,CAAC;IAET,MAAM,KAAK,GAAkB,EAAE,CAAC;IAChC,MAAM,KAAK,GAAkB,EAAE,CAAC;IAChC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC,IAAI,CAAC,EAAE;QAC3C,IAAI,QAAQ,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;QAC3B,IAAI,QAAQ,GAAG,QAAQ,CAAC,CAAC,CAAC,CAAC;QAC3B,IAAI,MAAM,GAAG,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAC7B,IAAI,MAAM,GAAG,QAAQ,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC;QAE7B,KAAK,CAAC,IAAI,CAAC,IAAA,iCAAS,EAAC,IAAA,iCAAS,EAAC,MAAM,EAAE,CAAC,CAAC,CAAC,EAAE,QAAQ,CAAC,CAAC,CAAC;QACvD,KAAK,CAAC,IAAI,CAAC,IAAA,iCAAS,EAAC,MAAM,EAAE,QAAQ,CAAC,CAAC,CAAC;QACxC,KAAK,CAAC,IAAI,CAAC,IAAA,4BAAI,EAAC,MAAM,EAAE,QAAQ,CAAC,CAAC,CAAC;QACnC,KAAK,CAAC,IAAI,CAAC,IAAA,iCAAS,EAAC,MAAM,EAAE,QAAQ,CAAC,CAAC,CAAC;KACzC;IAED,OAAO,EAAE,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,EAAE,CAAC;AAClC,CAAC;AAED,SAAS,aAAa,CAAC,IAAuB;IAC5C,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,IAAI,CAAC;IACxB,OAAO,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC,EAAE,CAAC;AAChC,CAAC;AAED,SAAS,cAAc,CAAC,IAAuB;IAC7C,MAAM,KAAK,GAAkB,EAAE,CAAC;IAChC,MAAM,KAAK,GAAkB,EAAE,CAAC;IAEhC,MAAM,EACJ,EAAE,EAAE,EAAE,CAAC,EAAE,QAAQ,EAAE,EACnB,EAAE,EAAE,EAAE,CAAC,EAAE,QAAQ,EAAE,GACpB,GAAG,IAAI,CAAC;IACT,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,QAAQ,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;QACxC,KAAK,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,IAAA,iCAAS,EAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;QACrE,KAAK,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,QAAQ,CAAC,CAAC,CAAC,CAAC,KAAK,EAAE,CAAC,CAAC,CAAC,IAAA,iCAAS,EAAC,QAAQ,CAAC,CAAC,CAAC,EAAE,CAAC,CAAC,CAAC,CAAC,CAAC;KACtE;IAED,OAAO,EAAE,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,EAAE,CAAC;AAClC,CAAC"}
@@ -0,0 +1,7 @@
1
+ import { DataReIm } from 'ml-spectra-processing';
2
+ export declare function digitalFilter(data: DataReIm, options: {
3
+ digitalFilterValue: number;
4
+ }): {
5
+ re: Float64Array;
6
+ im: Float64Array;
7
+ };
@@ -0,0 +1,20 @@
1
+ "use strict";
2
+ Object.defineProperty(exports, "__esModule", { value: true });
3
+ exports.digitalFilter = void 0;
4
+ function digitalFilter(data, options) {
5
+ let { digitalFilterValue = 0 } = options;
6
+ let re = new Float64Array(data.re);
7
+ let im = new Float64Array(data.im);
8
+ let pointsToShift = Math.floor(digitalFilterValue);
9
+ const skip = 0;
10
+ pointsToShift += 0;
11
+ const newRe = new Float64Array(re.length);
12
+ const newIm = new Float64Array(im.length);
13
+ newRe.set(re.slice(pointsToShift));
14
+ newRe.set(re.slice(skip, pointsToShift), re.length - pointsToShift);
15
+ newIm.set(im.slice(pointsToShift));
16
+ newIm.set(im.slice(skip, pointsToShift), im.length - pointsToShift);
17
+ return { re: newRe, im: newIm };
18
+ }
19
+ exports.digitalFilter = digitalFilter;
20
+ //# sourceMappingURL=digitalFilter.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"digitalFilter.js","sourceRoot":"","sources":["../../../../src/xyz/util/fft2d/digitalFilter.ts"],"names":[],"mappings":";;;AAEA,SAAgB,aAAa,CAC3B,IAAc,EACd,OAAuC;IAEvC,IAAI,EAAE,kBAAkB,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC;IACzC,IAAI,EAAE,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IACnC,IAAI,EAAE,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAEnC,IAAI,aAAa,GAAG,IAAI,CAAC,KAAK,CAAC,kBAAkB,CAAC,CAAC;IAEnD,MAAM,IAAI,GAAG,CAAC,CAAC;IACf,aAAa,IAAI,CAAC,CAAC;IAEnB,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC;IAC1C,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC;IAC1C,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,aAAa,CAAC,CAAC,CAAC;IACnC,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,IAAI,EAAE,aAAa,CAAC,EAAE,EAAE,CAAC,MAAM,GAAG,aAAa,CAAC,CAAC;IACpE,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,aAAa,CAAC,CAAC,CAAC;IACnC,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,IAAI,EAAE,aAAa,CAAC,EAAE,EAAE,CAAC,MAAM,GAAG,aAAa,CAAC,CAAC;IAEpE,OAAO,EAAE,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,EAAE,CAAC;AAClC,CAAC;AArBD,sCAqBC"}
@@ -0,0 +1,24 @@
1
+ import { DoubleArray } from 'cheminfo-types';
2
+ import Matrix from 'ml-matrix';
3
+ import { ApodizationOptions } from '../../../apodization/apodization';
4
+ import { ZeroFillingOptions } from './zeroFilling';
5
+ interface PhaseCorrection {
6
+ mode?: string;
7
+ ph0?: number;
8
+ ph1?: number;
9
+ }
10
+ export interface DirectDimensionOptions {
11
+ fnMode?: string;
12
+ digitalFilterValue?: number;
13
+ zeroFilling?: ZeroFillingOptions;
14
+ apodization?: ApodizationOptions;
15
+ phaseCorrection?: PhaseCorrection;
16
+ }
17
+ export declare function fftDirectDimension(data: {
18
+ re: DoubleArray[];
19
+ im: DoubleArray[];
20
+ }, options: DirectDimensionOptions): {
21
+ re: Matrix;
22
+ im: Matrix;
23
+ };
24
+ export {};
@@ -0,0 +1,57 @@
1
+ "use strict";
2
+ var __importDefault = (this && this.__importDefault) || function (mod) {
3
+ return (mod && mod.__esModule) ? mod : { "default": mod };
4
+ };
5
+ Object.defineProperty(exports, "__esModule", { value: true });
6
+ exports.fftDirectDimension = void 0;
7
+ const ml_matrix_1 = __importDefault(require("ml-matrix"));
8
+ const ml_spectra_processing_1 = require("ml-spectra-processing");
9
+ const apodization_1 = require("../../../apodization/apodization");
10
+ const digitalFilter_1 = require("./digitalFilter");
11
+ const removeDCOffset_1 = require("./removeDCOffset");
12
+ const zeroFilling_1 = require("./zeroFilling");
13
+ function fftDirectDimension(data, options) {
14
+ const { re: reData, im: imData } = data;
15
+ const nbRows = reData.length;
16
+ const nbColumns = reData[0].length;
17
+ const { zeroFilling: zeroFillingF2, digitalFilterValue = 0 } = options;
18
+ const { factor: f2Factor = 2 } = zeroFillingF2 || {};
19
+ const nbPointsF2 = (zeroFillingF2 === null || zeroFillingF2 === void 0 ? void 0 : zeroFillingF2.nbPoints) || 2 ** Math.round(Math.log2(nbColumns * f2Factor));
20
+ const reFFT = new ml_matrix_1.default(nbRows, nbPointsF2);
21
+ const imFFT = new ml_matrix_1.default(nbRows, nbPointsF2);
22
+ const { apodization: apodizationF2 } = options;
23
+ let { ph0: ph0F2, ph1: ph1F2, mode: pcModeF2 = 'no', } = (options === null || options === void 0 ? void 0 : options.phaseCorrection) || {};
24
+ for (let i = 0; i < nbRows; i++) {
25
+ const newData = (0, digitalFilter_1.digitalFilter)({ re: reData[i], im: imData[i] }, { digitalFilterValue });
26
+ Object.assign(newData, (0, removeDCOffset_1.removeDCOffset)(newData, { digitalFilterValue }));
27
+ if (apodizationF2) {
28
+ const { re: aRe, im: aIm } = (0, apodization_1.apodization)(newData, apodizationF2);
29
+ Object.assign(newData, { re: aRe, im: aIm });
30
+ }
31
+ Object.assign(newData, (0, zeroFilling_1.zeroFilling)(newData, { digitalFilterValue, nbPoints: nbPointsF2 }));
32
+ Object.assign(newData, (0, ml_spectra_processing_1.reimFFT)(newData, { applyZeroShift: true }));
33
+ const dfResidual = (digitalFilterValue - Math.round(digitalFilterValue)) * Math.PI * 2;
34
+ if (dfResidual > 0) {
35
+ Object.assign(newData, (0, ml_spectra_processing_1.reimPhaseCorrection)(newData, 0, -dfResidual));
36
+ }
37
+ if (i === 0 &&
38
+ pcModeF2 === 'pk' &&
39
+ ph0F2 === undefined &&
40
+ ph1F2 === undefined) {
41
+ const phased = (0, ml_spectra_processing_1.reimAutoPhaseCorrection)(newData, { magnitudeMode: false });
42
+ ph0F2 = phased.ph0;
43
+ ph1F2 = phased.ph1;
44
+ }
45
+ if (pcModeF2 === 'pk' && ph0F2 !== undefined && ph1F2 !== undefined) {
46
+ Object.assign(newData, (0, ml_spectra_processing_1.reimPhaseCorrection)(newData, (ph0F2 * Math.PI) / 180, (ph1F2 * Math.PI) / 180));
47
+ }
48
+ reFFT.setRow(i, newData.re);
49
+ imFFT.setRow(i, newData.im);
50
+ }
51
+ return {
52
+ re: reFFT,
53
+ im: imFFT,
54
+ };
55
+ }
56
+ exports.fftDirectDimension = fftDirectDimension;
57
+ //# sourceMappingURL=fftDirectDimension.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"fftDirectDimension.js","sourceRoot":"","sources":["../../../../src/xyz/util/fft2d/fftDirectDimension.ts"],"names":[],"mappings":";;;;;;AACA,0DAA+B;AAC/B,iEAI+B;AAE/B,kEAG0C;AAE1C,mDAAgD;AAChD,qDAAkD;AAClD,+CAAgE;AAgBhE,SAAgB,kBAAkB,CAChC,IAA8C,EAC9C,OAA+B;IAE/B,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,EAAE,EAAE,MAAM,EAAE,GAAG,IAAI,CAAC;IACxC,MAAM,MAAM,GAAG,MAAM,CAAC,MAAM,CAAC;IAC7B,MAAM,SAAS,GAAG,MAAM,CAAC,CAAC,CAAC,CAAC,MAAM,CAAC;IAEnC,MAAM,EAAE,WAAW,EAAE,aAAa,EAAE,kBAAkB,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC;IAEvE,MAAM,EAAE,MAAM,EAAE,QAAQ,GAAG,CAAC,EAAE,GAAG,aAAa,IAAI,EAAE,CAAC;IACrD,MAAM,UAAU,GACd,CAAA,aAAa,aAAb,aAAa,uBAAb,aAAa,CAAE,QAAQ,KAAI,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,SAAS,GAAG,QAAQ,CAAC,CAAC,CAAC;IAE9E,MAAM,KAAK,GAAG,IAAI,mBAAM,CAAC,MAAM,EAAE,UAAU,CAAC,CAAC;IAC7C,MAAM,KAAK,GAAG,IAAI,mBAAM,CAAC,MAAM,EAAE,UAAU,CAAC,CAAC;IAE7C,MAAM,EAAE,WAAW,EAAE,aAAa,EAAE,GAAG,OAAO,CAAC;IAC/C,IAAI,EACF,GAAG,EAAE,KAAK,EACV,GAAG,EAAE,KAAK,EACV,IAAI,EAAE,QAAQ,GAAG,IAAI,GACtB,GAAG,CAAA,OAAO,aAAP,OAAO,uBAAP,OAAO,CAAE,eAAe,KAAI,EAAE,CAAC;IAEnC,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,MAAM,EAAE,CAAC,EAAE,EAAE;QAC/B,MAAM,OAAO,GAAG,IAAA,6BAAa,EAC3B,EAAE,EAAE,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,EAAE,EAAE,MAAM,CAAC,CAAC,CAAC,EAAE,EAChC,EAAE,kBAAkB,EAAE,CACvB,CAAC;QAEF,MAAM,CAAC,MAAM,CAAC,OAAO,EAAE,IAAA,+BAAc,EAAC,OAAO,EAAE,EAAE,kBAAkB,EAAE,CAAC,CAAC,CAAC;QAExE,IAAI,aAAa,EAAE;YACjB,MAAM,EAAE,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,IAAA,yBAAW,EAAC,OAAO,EAAE,aAAa,CAAC,CAAC;YACjE,MAAM,CAAC,MAAM,CAAC,OAAO,EAAE,EAAE,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,CAAC,CAAC;SAC9C;QAED,MAAM,CAAC,MAAM,CACX,OAAO,EACP,IAAA,yBAAW,EAAC,OAAO,EAAE,EAAE,kBAAkB,EAAE,QAAQ,EAAE,UAAU,EAAE,CAAC,CACnE,CAAC;QAEF,MAAM,CAAC,MAAM,CAAC,OAAO,EAAE,IAAA,+BAAO,EAAC,OAAO,EAAE,EAAE,cAAc,EAAE,IAAI,EAAE,CAAC,CAAC,CAAC;QAEnE,MAAM,UAAU,GACd,CAAC,kBAAkB,GAAG,IAAI,CAAC,KAAK,CAAC,kBAAkB,CAAC,CAAC,GAAG,IAAI,CAAC,EAAE,GAAG,CAAC,CAAC;QACtE,IAAI,UAAU,GAAG,CAAC,EAAE;YAClB,MAAM,CAAC,MAAM,CAAC,OAAO,EAAE,IAAA,2CAAmB,EAAC,OAAO,EAAE,CAAC,EAAE,CAAC,UAAU,CAAC,CAAC,CAAC;SACtE;QAED,IACE,CAAC,KAAK,CAAC;YACP,QAAQ,KAAK,IAAI;YACjB,KAAK,KAAK,SAAS;YACnB,KAAK,KAAK,SAAS,EACnB;YACA,MAAM,MAAM,GAAG,IAAA,+CAAuB,EAAC,OAAO,EAAE,EAAE,aAAa,EAAE,KAAK,EAAE,CAAC,CAAC;YAC1E,KAAK,GAAG,MAAM,CAAC,GAAG,CAAC;YACnB,KAAK,GAAG,MAAM,CAAC,GAAG,CAAC;SACpB;QACD,IAAI,QAAQ,KAAK,IAAI,IAAI,KAAK,KAAK,SAAS,IAAI,KAAK,KAAK,SAAS,EAAE;YACnE,MAAM,CAAC,MAAM,CACX,OAAO,EACP,IAAA,2CAAmB,EACjB,OAAO,EACP,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,EACvB,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CACxB,CACF,CAAC;SACH;QAED,KAAK,CAAC,MAAM,CAAC,CAAC,EAAE,OAAO,CAAC,EAAE,CAAC,CAAC;QAC5B,KAAK,CAAC,MAAM,CAAC,CAAC,EAAE,OAAO,CAAC,EAAE,CAAC,CAAC;KAC7B;IAED,OAAO;QACL,EAAE,EAAE,KAAK;QACT,EAAE,EAAE,KAAK;KACV,CAAC;AACJ,CAAC;AA/ED,gDA+EC"}
@@ -0,0 +1,23 @@
1
+ import { DoubleArray } from 'cheminfo-types';
2
+ import Matrix from 'ml-matrix';
3
+ import { ApodizationOptions } from '../../../apodization/apodization';
4
+ import { ZeroFillingOptions } from './zeroFilling';
5
+ interface PhaseCorrection {
6
+ mode?: string;
7
+ ph0?: number;
8
+ ph1?: number;
9
+ }
10
+ export interface IndirectDimensionOptions {
11
+ fnMode?: string;
12
+ zeroFilling?: ZeroFillingOptions;
13
+ apodization?: ApodizationOptions;
14
+ phaseCorrection?: PhaseCorrection;
15
+ }
16
+ export declare function fftIndirectDimension(data: {
17
+ re: DoubleArray[] | Matrix;
18
+ im: DoubleArray[] | Matrix;
19
+ }, options: IndirectDimensionOptions): {
20
+ re: Matrix;
21
+ im: Matrix;
22
+ };
23
+ export {};
@@ -0,0 +1,72 @@
1
+ "use strict";
2
+ var __importDefault = (this && this.__importDefault) || function (mod) {
3
+ return (mod && mod.__esModule) ? mod : { "default": mod };
4
+ };
5
+ Object.defineProperty(exports, "__esModule", { value: true });
6
+ exports.fftIndirectDimension = void 0;
7
+ const ml_matrix_1 = __importDefault(require("ml-matrix"));
8
+ const ml_spectra_processing_1 = require("ml-spectra-processing");
9
+ const apodization_1 = require("../../../apodization/apodization");
10
+ const zeroFilling_1 = require("./zeroFilling");
11
+ function fftIndirectDimension(data, options) {
12
+ var _a;
13
+ const reData = ml_matrix_1.default.checkMatrix(data.re);
14
+ const imData = ml_matrix_1.default.checkMatrix(data.im);
15
+ const nbRows = reData.rows;
16
+ const nbColumns = reData.columns;
17
+ const { zeroFilling: zeroFillingF1 } = options;
18
+ const { factor: f1Factor = 2 } = zeroFillingF1 || {};
19
+ const nbPointsF1 = (zeroFillingF1 === null || zeroFillingF1 === void 0 ? void 0 : zeroFillingF1.nbPoints) || 2 ** Math.round(Math.log2(nbRows * f1Factor));
20
+ const { apodization: apodizationF1 } = options;
21
+ let { ph0: ph0F1, ph1: ph1F1, mode: pcModeF1, } = options.phaseCorrection || {};
22
+ const reResult = new ml_matrix_1.default(nbPointsF1, nbColumns);
23
+ const imResult = new ml_matrix_1.default(nbPointsF1, nbColumns);
24
+ for (let i = 0; i < nbColumns; i++) {
25
+ const { re: aRe, im: aIm } = (0, apodization_1.apodization)({
26
+ re: reData.getColumn(i),
27
+ im: imData.getColumn(i),
28
+ }, apodizationF1);
29
+ const newRe = [];
30
+ const newIm = [];
31
+ if (options.fnMode === 'echo') {
32
+ for (let j = 0; j < aRe.length; j += 2) {
33
+ newRe.push(aRe[j]);
34
+ newIm.push(aRe[j + 1]);
35
+ }
36
+ }
37
+ else {
38
+ for (let j = 0; j < aRe.length; j++) {
39
+ newRe.push(aRe[j]);
40
+ newIm.push(-aIm[j]);
41
+ }
42
+ }
43
+ let transformed = (0, ml_spectra_processing_1.reimFFT)((0, zeroFilling_1.zeroFilling)({
44
+ re: newRe,
45
+ im: newIm,
46
+ }, { digitalFilterValue: 0, nbPoints: nbPointsF1 }), { applyZeroShift: true });
47
+ let { re, im } = transformed;
48
+ if (((_a = options.phaseCorrection) === null || _a === void 0 ? void 0 : _a.mode) === 'mc') {
49
+ re = (0, ml_spectra_processing_1.reimAbsolute)({ re, im });
50
+ im = new Float64Array(im.length);
51
+ }
52
+ if (pcModeF1 === 'pk') {
53
+ if (i === 0 && ph0F1 === undefined && ph1F1 === undefined) {
54
+ //here we can improve it by checking if the experiment is phase sensitive.
55
+ const phased = (0, ml_spectra_processing_1.reimAutoPhaseCorrection)({ re, im }, { magnitudeMode: false });
56
+ ph0F1 = phased.ph0;
57
+ ph1F1 = phased.ph1;
58
+ }
59
+ if (ph0F1 !== undefined && ph1F1 !== undefined) {
60
+ Object.assign(transformed, (0, ml_spectra_processing_1.reimPhaseCorrection)(transformed, (ph0F1 * Math.PI) / 180, (ph1F1 * Math.PI) / 180));
61
+ }
62
+ }
63
+ reResult.setColumn(i, re);
64
+ imResult.setColumn(i, im);
65
+ }
66
+ return {
67
+ re: reResult,
68
+ im: imResult,
69
+ };
70
+ }
71
+ exports.fftIndirectDimension = fftIndirectDimension;
72
+ //# sourceMappingURL=fftIndirectDimension.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"fftIndirectDimension.js","sourceRoot":"","sources":["../../../../src/xyz/util/fft2d/fftIndirectDimension.ts"],"names":[],"mappings":";;;;;;AACA,0DAA+B;AAC/B,iEAK+B;AAE/B,kEAG0C;AAE1C,+CAAgE;AAehE,SAAgB,oBAAoB,CAClC,IAAgE,EAChE,OAAiC;;IAEjC,MAAM,MAAM,GAAG,mBAAM,CAAC,WAAW,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAC3C,MAAM,MAAM,GAAG,mBAAM,CAAC,WAAW,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IAE3C,MAAM,MAAM,GAAG,MAAM,CAAC,IAAI,CAAC;IAC3B,MAAM,SAAS,GAAG,MAAM,CAAC,OAAO,CAAC;IAEjC,MAAM,EAAE,WAAW,EAAE,aAAa,EAAE,GAAG,OAAO,CAAC;IAC/C,MAAM,EAAE,MAAM,EAAE,QAAQ,GAAG,CAAC,EAAE,GAAG,aAAa,IAAI,EAAE,CAAC;IACrD,MAAM,UAAU,GACd,CAAA,aAAa,aAAb,aAAa,uBAAb,aAAa,CAAE,QAAQ,KAAI,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,MAAM,GAAG,QAAQ,CAAC,CAAC,CAAC;IAC3E,MAAM,EAAE,WAAW,EAAE,aAAa,EAAE,GAAG,OAAO,CAAC;IAC/C,IAAI,EACF,GAAG,EAAE,KAAK,EACV,GAAG,EAAE,KAAK,EACV,IAAI,EAAE,QAAQ,GACf,GAAG,OAAO,CAAC,eAAe,IAAI,EAAE,CAAC;IAElC,MAAM,QAAQ,GAAG,IAAI,mBAAM,CAAC,UAAU,EAAE,SAAS,CAAC,CAAC;IACnD,MAAM,QAAQ,GAAG,IAAI,mBAAM,CAAC,UAAU,EAAE,SAAS,CAAC,CAAC;IACnD,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,SAAS,EAAE,CAAC,EAAE,EAAE;QAClC,MAAM,EAAE,EAAE,EAAE,GAAG,EAAE,EAAE,EAAE,GAAG,EAAE,GAAG,IAAA,yBAAW,EACtC;YACE,EAAE,EAAE,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC;YACvB,EAAE,EAAE,MAAM,CAAC,SAAS,CAAC,CAAC,CAAC;SACxB,EACD,aAAa,CACd,CAAC;QAEF,MAAM,KAAK,GAAG,EAAE,CAAC;QACjB,MAAM,KAAK,GAAG,EAAE,CAAC;QACjB,IAAI,OAAO,CAAC,MAAM,KAAK,MAAM,EAAE;YAC7B,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,IAAI,CAAC,EAAE;gBACtC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;gBACnB,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC;aACxB;SACF;aAAM;YACL,KAAK,IAAI,CAAC,GAAG,CAAC,EAAE,CAAC,GAAG,GAAG,CAAC,MAAM,EAAE,CAAC,EAAE,EAAE;gBACnC,KAAK,CAAC,IAAI,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;gBACnB,KAAK,CAAC,IAAI,CAAC,CAAC,GAAG,CAAC,CAAC,CAAC,CAAC,CAAC;aACrB;SACF;QAED,IAAI,WAAW,GAAG,IAAA,+BAAO,EACvB,IAAA,yBAAW,EACT;YACE,EAAE,EAAE,KAAK;YACT,EAAE,EAAE,KAAK;SACV,EACD,EAAE,kBAAkB,EAAE,CAAC,EAAE,QAAQ,EAAE,UAAU,EAAE,CAChD,EACD,EAAE,cAAc,EAAE,IAAI,EAAE,CACzB,CAAC;QACF,IAAI,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,WAAW,CAAC;QAE7B,IAAI,CAAA,MAAA,OAAO,CAAC,eAAe,0CAAE,IAAI,MAAK,IAAI,EAAE;YAC1C,EAAE,GAAG,IAAA,oCAAY,EAAC,EAAE,EAAE,EAAE,EAAE,EAAE,CAAC,CAAC;YAC9B,EAAE,GAAG,IAAI,YAAY,CAAC,EAAE,CAAC,MAAM,CAAC,CAAC;SAClC;QAED,IAAI,QAAQ,KAAK,IAAI,EAAE;YACrB,IAAI,CAAC,KAAK,CAAC,IAAI,KAAK,KAAK,SAAS,IAAI,KAAK,KAAK,SAAS,EAAE;gBACzD,0EAA0E;gBAC1E,MAAM,MAAM,GAAG,IAAA,+CAAuB,EACpC,EAAE,EAAE,EAAE,EAAE,EAAE,EACV,EAAE,aAAa,EAAE,KAAK,EAAE,CACzB,CAAC;gBACF,KAAK,GAAG,MAAM,CAAC,GAAG,CAAC;gBACnB,KAAK,GAAG,MAAM,CAAC,GAAG,CAAC;aACpB;YACD,IAAI,KAAK,KAAK,SAAS,IAAI,KAAK,KAAK,SAAS,EAAE;gBAC9C,MAAM,CAAC,MAAM,CACX,WAAW,EACX,IAAA,2CAAmB,EACjB,WAAW,EACX,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,EACvB,CAAC,KAAK,GAAG,IAAI,CAAC,EAAE,CAAC,GAAG,GAAG,CACxB,CACF,CAAC;aACH;SACF;QAED,QAAQ,CAAC,SAAS,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;QAC1B,QAAQ,CAAC,SAAS,CAAC,CAAC,EAAE,EAAE,CAAC,CAAC;KAC3B;IAED,OAAO;QACL,EAAE,EAAE,QAAQ;QACZ,EAAE,EAAE,QAAQ;KACb,CAAC;AACJ,CAAC;AA7FD,oDA6FC"}
@@ -0,0 +1,5 @@
1
+ import { DataReIm } from 'ml-spectra-processing';
2
+ export declare function removeDCOffset(data: DataReIm, options: any): {
3
+ re: Float64Array;
4
+ im: Float64Array;
5
+ };
@@ -0,0 +1,19 @@
1
+ "use strict";
2
+ Object.defineProperty(exports, "__esModule", { value: true });
3
+ exports.removeDCOffset = void 0;
4
+ const ml_spectra_processing_1 = require("ml-spectra-processing");
5
+ function removeDCOffset(data, options) {
6
+ let { digitalFilterValue = 0 } = options;
7
+ const nbPoints = data.re.length;
8
+ const newRe = new Float64Array(data.re);
9
+ const newIm = new Float64Array(data.im);
10
+ const averageRe = (0, ml_spectra_processing_1.xMean)(data.re.slice((nbPoints * 0.75) >> 0, nbPoints - digitalFilterValue));
11
+ const averageIm = (0, ml_spectra_processing_1.xMean)(data.im.slice((nbPoints * 0.75) >> 0, nbPoints - digitalFilterValue));
12
+ for (let i = digitalFilterValue; i < nbPoints - digitalFilterValue; i++) {
13
+ newRe[i] -= averageRe;
14
+ newIm[i] -= averageIm;
15
+ }
16
+ return { re: newRe, im: newIm };
17
+ }
18
+ exports.removeDCOffset = removeDCOffset;
19
+ //# sourceMappingURL=removeDCOffset.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"removeDCOffset.js","sourceRoot":"","sources":["../../../../src/xyz/util/fft2d/removeDCOffset.ts"],"names":[],"mappings":";;;AAAA,iEAAwD;AAExD,SAAgB,cAAc,CAAC,IAAc,EAAE,OAAY;IACzD,IAAI,EAAE,kBAAkB,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC;IACzC,MAAM,QAAQ,GAAG,IAAI,CAAC,EAAE,CAAC,MAAM,CAAC;IAChC,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IACxC,MAAM,KAAK,GAAG,IAAI,YAAY,CAAC,IAAI,CAAC,EAAE,CAAC,CAAC;IACxC,MAAM,SAAS,GAAG,IAAA,6BAAK,EACrB,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,QAAQ,GAAG,kBAAkB,CAAC,CACrE,CAAC;IACF,MAAM,SAAS,GAAG,IAAA,6BAAK,EACrB,IAAI,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,QAAQ,GAAG,IAAI,CAAC,IAAI,CAAC,EAAE,QAAQ,GAAG,kBAAkB,CAAC,CACrE,CAAC;IACF,KAAK,IAAI,CAAC,GAAG,kBAAkB,EAAE,CAAC,GAAG,QAAQ,GAAG,kBAAkB,EAAE,CAAC,EAAE,EAAE;QACvE,KAAK,CAAC,CAAC,CAAC,IAAI,SAAS,CAAC;QACtB,KAAK,CAAC,CAAC,CAAC,IAAI,SAAS,CAAC;KACvB;IACD,OAAO,EAAE,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,EAAE,CAAC;AAClC,CAAC;AAhBD,wCAgBC"}
@@ -0,0 +1,10 @@
1
+ import { DataReIm } from 'ml-spectra-processing';
2
+ export interface ZeroFillingOptions {
3
+ nbPoints?: number;
4
+ factor?: number;
5
+ digitalFilterValue?: number;
6
+ }
7
+ export declare function zeroFilling(data: DataReIm, options: ZeroFillingOptions): {
8
+ re: Float64Array;
9
+ im: Float64Array;
10
+ };
@@ -0,0 +1,29 @@
1
+ "use strict";
2
+ Object.defineProperty(exports, "__esModule", { value: true });
3
+ exports.zeroFilling = void 0;
4
+ function zeroFilling(data, options) {
5
+ let { nbPoints, factor = 2, digitalFilterValue: grpdly = 0 } = options;
6
+ if (!nbPoints) {
7
+ nbPoints = 2 ** Math.round(Math.log2(data.re.length * factor));
8
+ }
9
+ let pointsToShift;
10
+ if (grpdly > 0) {
11
+ pointsToShift = Math.floor(grpdly);
12
+ }
13
+ else {
14
+ pointsToShift = 0;
15
+ }
16
+ const { re, im } = data;
17
+ let newRE = new Float64Array(nbPoints);
18
+ let newIM = new Float64Array(nbPoints);
19
+ const length = Math.min(nbPoints, re.length);
20
+ newRE.set(re.slice(0, length - pointsToShift));
21
+ newIM.set(im.slice(0, length - pointsToShift));
22
+ if (pointsToShift > 0 && pointsToShift < nbPoints) {
23
+ newRE.set(re.slice(re.length - pointsToShift), nbPoints - pointsToShift);
24
+ newIM.set(im.slice(re.length - pointsToShift), nbPoints - pointsToShift);
25
+ }
26
+ return { re: newRE, im: newIM };
27
+ }
28
+ exports.zeroFilling = zeroFilling;
29
+ //# sourceMappingURL=zeroFilling.js.map
@@ -0,0 +1 @@
1
+ {"version":3,"file":"zeroFilling.js","sourceRoot":"","sources":["../../../../src/xyz/util/fft2d/zeroFilling.ts"],"names":[],"mappings":";;;AAQA,SAAgB,WAAW,CAAC,IAAc,EAAE,OAA2B;IACrE,IAAI,EAAE,QAAQ,EAAE,MAAM,GAAG,CAAC,EAAE,kBAAkB,EAAE,MAAM,GAAG,CAAC,EAAE,GAAG,OAAO,CAAC;IAEvE,IAAI,CAAC,QAAQ,EAAE;QACb,QAAQ,GAAG,CAAC,IAAI,IAAI,CAAC,KAAK,CAAC,IAAI,CAAC,IAAI,CAAC,IAAI,CAAC,EAAE,CAAC,MAAM,GAAG,MAAM,CAAC,CAAC,CAAC;KAChE;IAED,IAAI,aAAa,CAAC;IAClB,IAAI,MAAM,GAAG,CAAC,EAAE;QACd,aAAa,GAAG,IAAI,CAAC,KAAK,CAAC,MAAM,CAAC,CAAC;KACpC;SAAM;QACL,aAAa,GAAG,CAAC,CAAC;KACnB;IAED,MAAM,EAAE,EAAE,EAAE,EAAE,EAAE,GAAG,IAAI,CAAC;IAExB,IAAI,KAAK,GAAG,IAAI,YAAY,CAAC,QAAQ,CAAC,CAAC;IACvC,IAAI,KAAK,GAAG,IAAI,YAAY,CAAC,QAAQ,CAAC,CAAC;IAEvC,MAAM,MAAM,GAAG,IAAI,CAAC,GAAG,CAAC,QAAQ,EAAE,EAAE,CAAC,MAAM,CAAC,CAAC;IAE7C,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,MAAM,GAAG,aAAa,CAAC,CAAC,CAAC;IAC/C,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,CAAC,EAAE,MAAM,GAAG,aAAa,CAAC,CAAC,CAAC;IAE/C,IAAI,aAAa,GAAG,CAAC,IAAI,aAAa,GAAG,QAAQ,EAAE;QACjD,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,MAAM,GAAG,aAAa,CAAC,EAAE,QAAQ,GAAG,aAAa,CAAC,CAAC;QACzE,KAAK,CAAC,GAAG,CAAC,EAAE,CAAC,KAAK,CAAC,EAAE,CAAC,MAAM,GAAG,aAAa,CAAC,EAAE,QAAQ,GAAG,aAAa,CAAC,CAAC;KAC1E;IAED,OAAO,EAAE,EAAE,EAAE,KAAK,EAAE,EAAE,EAAE,KAAK,EAAE,CAAC;AAClC,CAAC;AA9BD,kCA8BC"}