| 1 | /**
|
| 2 | * @license
|
| 3 | * Copyright 2021 Google LLC. All Rights Reserved.
|
| 4 | * Licensed under the Apache License, Version 2.0 (the "License");
|
| 5 | * you may not use this file except in compliance with the License.
|
| 6 | * You may obtain a copy of the License at
|
| 7 | *
|
| 8 | * http://www.apache.org/licenses/LICENSE-2.0
|
| 9 | *
|
| 10 | * Unless required by applicable law or agreed to in writing, software
|
| 11 | * distributed under the License is distributed on an "AS IS" BASIS,
|
| 12 | * WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
|
| 13 | * See the License for the specific language governing permissions and
|
| 14 | * limitations under the License.
|
| 15 | * =============================================================================
|
| 16 | */
|
| 17 | import { ENGINE } from '../engine';
|
| 18 | import { Einsum } from '../kernel_names';
|
| 19 | import { convertToTensor } from '../tensor_util_env';
|
| 20 | import { op } from './operation';
|
| 21 | /**
|
| 22 | * Tensor contraction over specified indices and outer product.
|
| 23 | *
|
| 24 | * `einsum` allows defining Tensors by defining their element-wise computation.
|
| 25 | * This computation is based on
|
| 26 | * [Einstein summation](https://en.wikipedia.org/wiki/Einstein_notation).
|
| 27 | *
|
| 28 | * Some special cases include:
|
| 29 | *
|
| 30 | * Matrix multiplication:
|
| 31 | * ```js
|
| 32 | * const x = tf.tensor2d([[1, 2, 3], [4, 5, 6]]);
|
| 33 | * const y = tf.tensor2d([[0, 1], [2, 3], [4, 5]]);
|
| 34 | * x.print();
|
| 35 | * y.print();
|
| 36 | * tf.einsum('ij,jk->ik', x, y).print();
|
| 37 | * ```
|
| 38 | *
|
| 39 | * Dot product:
|
| 40 | * ```js
|
| 41 | * const x = tf.tensor1d([1, 2, 3]);
|
| 42 | * const y = tf.tensor1d([0, 1, 2]);
|
| 43 | * x.print();
|
| 44 | * y.print();
|
| 45 | * tf.einsum('i,i->', x, y).print();
|
| 46 | * ```
|
| 47 | *
|
| 48 | * Batch dot product:
|
| 49 | * ```js
|
| 50 | * const x = tf.tensor2d([[1, 2, 3], [4, 5, 6]]);
|
| 51 | * const y = tf.tensor2d([[0, 1, 2], [3, 4, 5]]);
|
| 52 | * x.print();
|
| 53 | * y.print();
|
| 54 | * tf.einsum('bi,bi->b', x, y).print();
|
| 55 | * ```
|
| 56 | *
|
| 57 | * Outer prouduct:
|
| 58 | * ```js
|
| 59 | * const x = tf.tensor1d([1, 3, 5]);
|
| 60 | * const y = tf.tensor1d([2, 4, 6]);
|
| 61 | * x.print();
|
| 62 | * y.print();
|
| 63 | * tf.einsum('i,j->ij', x, y).print();
|
| 64 | * ```
|
| 65 | *
|
| 66 | * Matrix transpose:
|
| 67 | * ```js
|
| 68 | * const x = tf.tensor2d([[1, 2], [3, 4]]);
|
| 69 | * x.print();
|
| 70 | * tf.einsum('ij->ji', x).print();
|
| 71 | * ```
|
| 72 | *
|
| 73 | * Batch matrix transpose:
|
| 74 | * ```js
|
| 75 | * const x = tf.tensor3d([[[1, 2], [3, 4]], [[-1, -2], [-3, -4]]]);
|
| 76 | * x.print();
|
| 77 | * tf.einsum('bij->bji', x).print();
|
| 78 | * ```
|
| 79 | *
|
| 80 | * Limitations:
|
| 81 | *
|
| 82 | * This implementation of einsum has the following limitations:
|
| 83 | *
|
| 84 | * - Does not support >2 input tensors.
|
| 85 | * - Does not support duplicate axes for any given input tensor. E.g., equation
|
| 86 | * 'ii->' is not suppoted.
|
| 87 | * - The `...` notation is not supported.
|
| 88 | *
|
| 89 | * @param equation a string describing the contraction, in the same format as
|
| 90 | * [numpy.einsum](https://numpy.org/doc/stable/reference/generated/numpy.einsum.html).
|
| 91 | * @param tensors the input(s) to contract (each one a Tensor), whose shapes
|
| 92 | * should be consistent with equation.
|
| 93 | * @returns The output tensor.
|
| 94 | *
|
| 95 | * @doc {heading: 'Tensors', subheading: 'Matrices'}
|
| 96 | */
|
| 97 | export function einsum_(equation, ...tensors) {
|
| 98 | const $tensors = tensors.map((t, i) => convertToTensor(t, `tensors${i}`, 'einsum'));
|
| 99 | const attrs = { equation };
|
| 100 | return ENGINE.runKernel(Einsum, $tensors, attrs);
|
| 101 | }
|
| 102 | export const einsum = op({ einsum_ });
|
| 103 | //# sourceMappingURL=data:application/json;base64,eyJ2ZXJzaW9uIjozLCJmaWxlIjoiZWluc3VtLmpzIiwic291cmNlUm9vdCI6IiIsInNvdXJjZXMiOlsiLi4vLi4vLi4vLi4vLi4vLi4vdGZqcy1jb3JlL3NyYy9vcHMvZWluc3VtLnRzIl0sIm5hbWVzIjpbXSwibWFwcGluZ3MiOiJBQUFBOzs7Ozs7Ozs7Ozs7Ozs7R0FlRztBQUVILE9BQU8sRUFBQyxNQUFNLEVBQUMsTUFBTSxXQUFXLENBQUM7QUFDakMsT0FBTyxFQUFDLE1BQU0sRUFBYyxNQUFNLGlCQUFpQixDQUFDO0FBSXBELE9BQU8sRUFBQyxlQUFlLEVBQUMsTUFBTSxvQkFBb0IsQ0FBQztBQUVuRCxPQUFPLEVBQUMsRUFBRSxFQUFDLE1BQU0sYUFBYSxDQUFDO0FBRS9COzs7Ozs7Ozs7Ozs7Ozs7Ozs7Ozs7Ozs7Ozs7Ozs7Ozs7Ozs7Ozs7Ozs7Ozs7Ozs7Ozs7Ozs7Ozs7Ozs7Ozs7Ozs7Ozs7Ozs7Ozs7R0EyRUc7QUFDSCxNQUFNLFVBQVUsT0FBTyxDQUFDLFFBQWdCLEVBQUUsR0FBRyxPQUFpQjtJQUM1RCxNQUFNLFFBQVEsR0FDVixPQUFPLENBQUMsR0FBRyxDQUFDLENBQUMsQ0FBQyxFQUFFLENBQUMsRUFBRSxFQUFFLENBQUMsZUFBZSxDQUFDLENBQUMsRUFBRSxVQUFVLENBQUMsRUFBRSxFQUFFLFFBQVEsQ0FBQyxDQUFDLENBQUM7SUFDdkUsTUFBTSxLQUFLLEdBQWdCLEVBQUMsUUFBUSxFQUFDLENBQUM7SUFDdEMsT0FBTyxNQUFNLENBQUMsU0FBUyxDQUNuQixNQUFNLEVBQUUsUUFBZ0MsRUFBRSxLQUEyQixDQUFDLENBQUM7QUFDN0UsQ0FBQztBQUVELE1BQU0sQ0FBQyxNQUFNLE1BQU0sR0FBRyxFQUFFLENBQUMsRUFBQyxPQUFPLEVBQUMsQ0FBQyxDQUFDIiwic291cmNlc0NvbnRlbnQiOlsiLyoqXG4gKiBAbGljZW5zZVxuICogQ29weXJpZ2h0IDIwMjEgR29vZ2xlIExMQy4gQWxsIFJpZ2h0cyBSZXNlcnZlZC5cbiAqIExpY2Vuc2VkIHVuZGVyIHRoZSBBcGFjaGUgTGljZW5zZSwgVmVyc2lvbiAyLjAgKHRoZSBcIkxpY2Vuc2VcIik7XG4gKiB5b3UgbWF5IG5vdCB1c2UgdGhpcyBmaWxlIGV4Y2VwdCBpbiBjb21wbGlhbmNlIHdpdGggdGhlIExpY2Vuc2UuXG4gKiBZb3UgbWF5IG9idGFpbiBhIGNvcHkgb2YgdGhlIExpY2Vuc2UgYXRcbiAqXG4gKiBodHRwOi8vd3d3LmFwYWNoZS5vcmcvbGljZW5zZXMvTElDRU5TRS0yLjBcbiAqXG4gKiBVbmxlc3MgcmVxdWlyZWQgYnkgYXBwbGljYWJsZSBsYXcgb3IgYWdyZWVkIHRvIGluIHdyaXRpbmcsIHNvZnR3YXJlXG4gKiBkaXN0cmlidXRlZCB1bmRlciB0aGUgTGljZW5zZSBpcyBkaXN0cmlidXRlZCBvbiBhbiBcIkFTIElTXCIgQkFTSVMsXG4gKiBXSVRIT1VUIFdBUlJBTlRJRVMgT1IgQ09ORElUSU9OUyBPRiBBTlkgS0lORCwgZWl0aGVyIGV4cHJlc3Mgb3IgaW1wbGllZC5cbiAqIFNlZSB0aGUgTGljZW5zZSBmb3IgdGhlIHNwZWNpZmljIGxhbmd1YWdlIGdvdmVybmluZyBwZXJtaXNzaW9ucyBhbmRcbiAqIGxpbWl0YXRpb25zIHVuZGVyIHRoZSBMaWNlbnNlLlxuICogPT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT09PT1cbiAqL1xuXG5pbXBvcnQge0VOR0lORX0gZnJvbSAnLi4vZW5naW5lJztcbmltcG9ydCB7RWluc3VtLCBFaW5zdW1BdHRyc30gZnJvbSAnLi4va2VybmVsX25hbWVzJztcbmltcG9ydCB7TmFtZWRBdHRyTWFwfSBmcm9tICcuLi9rZXJuZWxfcmVnaXN0cnknO1xuaW1wb3J0IHtUZW5zb3J9IGZyb20gJy4uL3RlbnNvcic7XG5pbXBvcnQge05hbWVkVGVuc29yTWFwfSBmcm9tICcuLi90ZW5zb3JfdHlwZXMnO1xuaW1wb3J0IHtjb252ZXJ0VG9UZW5zb3J9IGZyb20gJy4uL3RlbnNvcl91dGlsX2Vudic7XG5cbmltcG9ydCB7b3B9IGZyb20gJy4vb3BlcmF0aW9uJztcblxuLyoqXG4gKiBUZW5zb3IgY29udHJhY3Rpb24gb3ZlciBzcGVjaWZpZWQgaW5kaWNlcyBhbmQgb3V0ZXIgcHJvZHVjdC5cbiAqXG4gKiBgZWluc3VtYCBhbGxvd3MgZGVmaW5pbmcgVGVuc29ycyBieSBkZWZpbmluZyB0aGVpciBlbGVtZW50LXdpc2UgY29tcHV0YXRpb24uXG4gKiBUaGlzIGNvbXB1dGF0aW9uIGlzIGJhc2VkIG9uXG4gKiBbRWluc3RlaW4gc3VtbWF0aW9uXShodHRwczovL2VuLndpa2lwZWRpYS5vcmcvd2lraS9FaW5zdGVpbl9ub3RhdGlvbikuXG4gKlxuICogU29tZSBzcGVjaWFsIGNhc2VzIGluY2x1ZGU6XG4gKlxuICogTWF0cml4IG11bHRpcGxpY2F0aW9uOlxuICogYGBganNcbiAqIGNvbnN0IHggPSB0Zi50ZW5zb3IyZChbWzEsIDIsIDNdLCBbNCwgNSwgNl1dKTtcbiAqIGNvbnN0IHkgPSB0Zi50ZW5zb3IyZChbWzAsIDFdLCBbMiwgM10sIFs0LCA1XV0pO1xuICogeC5wcmludCgpO1xuICogeS5wcmludCgpO1xuICogdGYuZWluc3VtKCdpaixqay0+aWsnLCB4LCB5KS5wcmludCgpO1xuICogYGBgXG4gKlxuICogRG90IHByb2R1Y3Q6XG4gKiBgYGBqc1xuICogY29uc3QgeCA9IHRmLnRlbnNvcjFkKFsxLCAyLCAzXSk7XG4gKiBjb25zdCB5ID0gdGYudGVuc29yMWQoWzAsIDEsIDJdKTtcbiAqIHgucHJpbnQoKTtcbiAqIHkucHJpbnQoKTtcbiAqIHRmLmVpbnN1bSgnaSxpLT4nLCB4LCB5KS5wcmludCgpO1xuICogYGBgXG4gKlxuICogQmF0Y2ggZG90IHByb2R1Y3Q6XG4gKiBgYGBqc1xuICogY29uc3QgeCA9IHRmLnRlbnNvcjJkKFtbMSwgMiwgM10sIFs0LCA1LCA2XV0pO1xuICogY29uc3QgeSA9IHRmLnRlbnNvcjJkKFtbMCwgMSwgMl0sIFszLCA0LCA1XV0pO1xuICogeC5wcmludCgpO1xuICogeS5wcmludCgpO1xuICogdGYuZWluc3VtKCdiaSxiaS0+YicsIHgsIHkpLnByaW50KCk7XG4gKiBgYGBcbiAqXG4gKiBPdXRlciBwcm91ZHVjdDpcbiAqIGBgYGpzXG4gKiBjb25zdCB4ID0gdGYudGVuc29yMWQoWzEsIDMsIDVdKTtcbiAqIGNvbnN0IHkgPSB0Zi50ZW5zb3IxZChbMiwgNCwgNl0pO1xuICogeC5wcmludCgpO1xuICogeS5wcmludCgpO1xuICogdGYuZWluc3VtKCdpLGotPmlqJywgeCwgeSkucHJpbnQoKTtcbiAqIGBgYFxuICpcbiAqIE1hdHJpeCB0cmFuc3Bvc2U6XG4gKiBgYGBqc1xuICogY29uc3QgeCA9IHRmLnRlbnNvcjJkKFtbMSwgMl0sIFszLCA0XV0pO1xuICogeC5wcmludCgpO1xuICogdGYuZWluc3VtKCdpai0+amknLCB4KS5wcmludCgpO1xuICogYGBgXG4gKlxuICogQmF0Y2ggbWF0cml4IHRyYW5zcG9zZTpcbiAqIGBgYGpzXG4gKiBjb25zdCB4ID0gdGYudGVuc29yM2QoW1tbMSwgMl0sIFszLCA0XV0sIFtbLTEsIC0yXSwgWy0zLCAtNF1dXSk7XG4gKiB4LnByaW50KCk7XG4gKiB0Zi5laW5zdW0oJ2Jpai0+YmppJywgeCkucHJpbnQoKTtcbiAqIGBgYFxuICpcbiAqIExpbWl0YXRpb25zOlxuICpcbiAqIFRoaXMgaW1wbGVtZW50YXRpb24gb2YgZWluc3VtIGhhcyB0aGUgZm9sbG93aW5nIGxpbWl0YXRpb25zOlxuICpcbiAqIC0gRG9lcyBub3Qgc3VwcG9ydCA+MiBpbnB1dCB0ZW5zb3JzLlxuICogLSBEb2VzIG5vdCBzdXBwb3J0IGR1cGxpY2F0ZSBheGVzIGZvciBhbnkgZ2l2ZW4gaW5wdXQgdGVuc29yLiBFLmcuLCBlcXVhdGlvblxuICogICAnaWktPicgaXMgbm90IHN1cHBvdGVkLlxuICogLSBUaGUgYC4uLmAgbm90YXRpb24gaXMgbm90IHN1cHBvcnRlZC5cbiAqXG4gKiBAcGFyYW0gZXF1YXRpb24gYSBzdHJpbmcgZGVzY3JpYmluZyB0aGUgY29udHJhY3Rpb24sIGluIHRoZSBzYW1lIGZvcm1hdCBhc1xuICogW251bXB5LmVpbnN1bV0oaHR0cHM6Ly9udW1weS5vcmcvZG9jL3N0YWJsZS9yZWZlcmVuY2UvZ2VuZXJhdGVkL251bXB5LmVpbnN1bS5odG1sKS5cbiAqIEBwYXJhbSB0ZW5zb3JzIHRoZSBpbnB1dChzKSB0byBjb250cmFjdCAoZWFjaCBvbmUgYSBUZW5zb3IpLCB3aG9zZSBzaGFwZXNcbiAqICAgICBzaG91bGQgYmUgY29uc2lzdGVudCB3aXRoIGVxdWF0aW9uLlxuICogQHJldHVybnMgVGhlIG91dHB1dCB0ZW5zb3IuXG4gKlxuICogQGRvYyB7aGVhZGluZzogJ1RlbnNvcnMnLCBzdWJoZWFkaW5nOiAnTWF0cmljZXMnfVxuICovXG5leHBvcnQgZnVuY3Rpb24gZWluc3VtXyhlcXVhdGlvbjogc3RyaW5nLCAuLi50ZW5zb3JzOiBUZW5zb3JbXSk6IFRlbnNvciB7XG4gIGNvbnN0ICR0ZW5zb3JzID1cbiAgICAgIHRlbnNvcnMubWFwKCh0LCBpKSA9PiBjb252ZXJ0VG9UZW5zb3IodCwgYHRlbnNvcnMke2l9YCwgJ2VpbnN1bScpKTtcbiAgY29uc3QgYXR0cnM6IEVpbnN1bUF0dHJzID0ge2VxdWF0aW9ufTtcbiAgcmV0dXJuIEVOR0lORS5ydW5LZXJuZWwoXG4gICAgICBFaW5zdW0sICR0ZW5zb3JzIGFzIHt9IGFzIE5hbWVkVGVuc29yTWFwLCBhdHRycyBhcyB7fSBhcyBOYW1lZEF0dHJNYXApO1xufVxuXG5leHBvcnQgY29uc3QgZWluc3VtID0gb3Aoe2VpbnN1bV99KTtcbiJdfQ== |
| \ | No newline at end of file |