def tdot

bendtorch-gpu@1.3.0.0 · main.bend:658

import Base
import bendtorch@1.3.0.0/main.bend as T

658# acc + a[ia]·b[ib] + .. down column j of a head's packed triangle (row i + 1's cell is i + 1 after row i's) and b's
659# rows (sb apart), from row i on: dot0's sum, but with the step in a growing by one a row
660def tdot(k: Nat, +ia: U32, +ib: U32, +i: U32, +sb: U32, acc: F32, ra: Array<F32> & F32, rb: Array<F32> & F32) -> T.D2:
661  match k ra rb:
662    case 0n Tuple{a, x} Tuple{b, y}: T.D2{a, b, acc}
663    case 1n+p Tuple{a, +x} Tuple{b, +y}: tdot(p, (ia + (i + 1 : U32) : U32), (ib + sb : U32), (i + 1 : U32), sb, (acc + (x * y : F32) : F32), Array.get(F32, a, (ia + (i + 1 : U32) : U32)), Array.get(F32, b, (ib + sb : U32)))