~/jendhub

def lnb_sums

bendtorch@1.3.0.0 · main.bend:1761

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

1761# sa + dxh[j] and sb + dxh[j]·xh[j], j = jj..jj + k, each in order
1762def lnb_sums(k: Nat, +j: U32, +off: U32, +mu: F32, +r: F32, sa: F32, sb: F32, rx: Array<F32> & F32, rd: Array<F32> & F32, rg: Array<F32> & F32) -> R3SS:
1763  match k rx rd rg:
1764    case 0n Tuple{x, +a} Tuple{dy, +b} Tuple{g, +c}:
1765      R3SS{x, dy, g, (sa + (b * c : F32) : F32), (sb + ((b * c : F32) * ((a - mu : F32) * r : F32) : F32) : F32)}
1766    case 1n+q Tuple{x, +a} Tuple{dy, +b} Tuple{g, +c}:
1767      lnb_sums(q, (j + 1 : U32), off, mu, r, (sa + (b * c : F32) : F32), (sb + ((b * c : F32) * ((a - mu : F32) * r : F32) : F32) : F32), Array.get(F32, x, ((off + j : U32) + 1 : U32)), Array.get(F32, dy, ((off + j : U32) + 1 : U32)), Array.get(F32, g, (j + 1 : U32)))

Uses (12)

Used by (4)