forked from tensorflow/tfjs
-
Notifications
You must be signed in to change notification settings - Fork 0
/
Cumprod.ts
88 lines (74 loc) · 3.1 KB
/
Cumprod.ts
1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
/**
* @license
* Copyright 2022 Google LLC. All Rights Reserved.
* Licensed under the Apache License, Version 2.0 (the "License");
* you may not use this file except in compliance with the License.
* You may obtain a copy of the License at
*
* http://www.apache.org/licenses/LICENSE-2.0
*
* Unless required by applicable law or agreed to in writing, software
* distributed under the License is distributed on an "AS IS" BASIS,
* WITHOUT WARRANTIES OR CONDITIONS OF ANY KIND, either express or implied.
* See the License for the specific language governing permissions and
* limitations under the License.
* =============================================================================
*/
import {backend_util, Cumprod, CumprodAttrs, CumprodInputs, KernelConfig, KernelFunc, TensorInfo, TypedArray, upcastType, util} from '@tensorflow/tfjs-core';
import {MathBackendCPU} from '../backend_cpu';
import {assertNotComplex} from '../cpu_util';
import {transpose} from './Transpose';
export function cumprod(
args: {inputs: CumprodInputs, backend: MathBackendCPU,
attrs: CumprodAttrs}): TensorInfo {
const {inputs, backend, attrs} = args;
const {x} = inputs;
const {axis, exclusive, reverse} = attrs;
assertNotComplex(x, 'cumprod');
const permutation = backend_util.getAxesPermutation([axis], x.shape.length);
let $x = x;
if (permutation != null) {
$x = transpose({inputs: {x}, backend, attrs: {perm: permutation}});
}
const permutedAxis = backend_util.getInnerMostAxes(1, x.shape.length)[0];
if (permutedAxis !== $x.shape.length - 1) {
throw new Error(
`backend.cumprod in CPU expects an inner-most ` +
`axis=${$x.shape.length - 1} but got axis=${permutedAxis}`);
}
const resultDtype = upcastType($x.dtype, 'int32');
const vals = util.makeOnesTypedArray(
util.sizeFromShape($x.shape), resultDtype) as TypedArray;
const aVals = backend.data.get($x.dataId).values as TypedArray;
const finalDim = $x.shape[$x.shape.length - 1];
const indexAdjuster = reverse ?
(i: number, j: number) => i + finalDim - j - 1 :
(i: number, j: number) => i + j;
for (let i = 0; i < aVals.length; i += finalDim) {
for (let j = 0; j < finalDim; j++) {
const idx = indexAdjuster(i, j);
if (j === 0) {
vals[idx] = exclusive ? 1 : aVals[idx];
} else {
const prevIdx = indexAdjuster(i, j - 1);
vals[idx] = exclusive ? aVals[prevIdx] * vals[prevIdx] :
aVals[idx] * vals[prevIdx];
}
}
}
const result = backend.makeTensorInfo($x.shape, resultDtype, vals);
if (permutation != null) {
const reversePermutation = backend_util.getUndoAxesPermutation(permutation);
const reverseTransposedResult = transpose(
{inputs: {x: result}, backend, attrs: {perm: reversePermutation}});
backend.disposeIntermediateTensorInfo(result);
backend.disposeIntermediateTensorInfo($x);
return reverseTransposedResult;
}
return result;
}
export const cumprodConfig: KernelConfig = {
kernelName: Cumprod,
backendName: 'cpu',
kernelFunc: cumprod as unknown as KernelFunc
};