함수형 갱신: 직접 바꾸지 말고 새 값을 돌려주기

NumPy의 제자리 변경을 JAX의 함수형 갱신으로 옮기는 패턴은 매일 쓰게 될 거야.

핵심 도구가 .at[] 인덱서:

import jax.numpy as jnp

a = jnp.zeros((4, 4))

# .at[index].set(value) — 새 array 반환, 원본 안 변함
b = a.at[0, 0].set(1.0)
print(a)  # 여전히 0
print(b)  # (0,0)에 1

# 다른 update operation 들
a.at[0].add(5.0)        # +=
a.at[0].multiply(2.0)   # *=
a.at[0].divide(2.0)     # /=
a.at[0].power(2.0)      # **=
a.at[0].min(0.0)        # element-wise min
a.at[0].max(1.0)        # element-wise max
a.at[0].apply(jnp.sqrt) # 임의 함수 적용

모든 인덱싱 형태 지원:

a = jnp.zeros((5, 5))

# slice
b = a.at[1:3, :].set(jnp.ones((2, 5)))

# fancy indexing
b = a.at[[0, 2, 4]].set(1.0)

# boolean mask
mask = jnp.eye(5).astype(bool)
b = a.at[mask].set(99.0)

# 체이닝
b = a.at[0, 0].set(1).at[1, 1].set(2).at[2, 2].set(3)

scatter / gather 같은 좀 더 복잡한 갱신:

# 같은 index 에 여러 번 더하기 — 누적
indices = jnp.array([0, 0, 1, 2, 0])
values = jnp.array([1.0, 2.0, 3.0, 4.0, 5.0])
result = jnp.zeros(5).at[indices].add(values)
print(result)  # [8.0, 3.0, 4.0, 0.0, 0.0] — index 0 에 1+2+5=8

# scatter (segment_sum)
result2 = jax.ops.segment_sum(values, indices, num_segments=5)

NumPy 패턴 → JAX 패턴 매핑:

NumPy                       | JAX
----------------------------|---------------------------
a[i] = v                    | a = a.at[i].set(v)
a[i] += v                   | a = a.at[i].add(v)
a[i:j] = b                  | a = a.at[i:j].set(b)
a[mask] = 0                 | a = a.at[mask].set(0)
                            | (또는 jnp.where(mask, 0, a))
np.add.at(a, idx, vals)     | a.at[idx].add(vals)

💡 .at vs jnp.where

불리언 마스크로 갱신할 때, jnp.where(mask, new, old)가 보통 더 깔끔하고 빨라. .at은 정확한 인덱스를 알 때 (특히 스캐터) 강력해. 둘 다 jit + grad + vmap 친화적이라, 둘 중 어느 걸 써도 큰 문제 없어. 가독성을 우선하면 돼.

처음 .at 패턴을 보면 "왜 이렇게 길게 써야 하지?" 싶은데, 익으면 NumPy의 a[i] = v가 실제로 얼마나 위험한 추상화였는지 보이게 돼. 부수 효과를 한 줄짜리 문법으로 숨겨 놓은 거였거든. JAX는 그 비용을 눈에 보이게 드러낸 거야.

Code

import jax.numpy as jnp

# IMPERATIVE (NumPy) style — won't work in JAX
# array[i] = value
# dict[key] = value
# object.attribute = value

# FUNCTIONAL (JAX) style — creates new values
# new_array = array.at[i].set(value)
# new_dict = {**dict, key: value}
# new_params = params._replace(attribute=value)  # NamedTuple

import jax.numpy as jnp

x = jnp.zeros((5, 5))

# Set values at specific positions
x = x.at[0, 0].set(1.0)
x = x.at[2, 3].set(5.0)
x = x.at[4, 4].set(1.0)

# Increment values
counts = jnp.zeros(10, dtype=jnp.int32)
data = jnp.array([3, 1, 4, 1, 5, 9, 2, 6, 5, 3])
# Build a histogram — add 1 at each index
histogram = counts.at[data].add(1)
print(histogram)  # [0 2 1 2 1 2 0 0 0 1]
# Note: duplicate indices accumulate (both 1s add, both 3s add, etc.)

# Scatter operations
indices = jnp.array([0, 2, 4])
values = jnp.array([10.0, 20.0, 30.0])
x = jnp.zeros(6)
x = x.at[indices].set(values)
print(x)  # [10.  0. 20.  0. 30.  0.]

import jax.numpy as jnp

x = jnp.array([-3.0, -1.0, 0.0, 2.0, 5.0])

# ReLU: replace negatives with 0
relu = jnp.where(x > 0, x, 0.0)
print(relu)  # [0. 0. 0. 2. 5.]

# Clamp between -1 and 1
clamped = jnp.where(x > 1, 1.0, jnp.where(x < -1, -1.0, x))
print(clamped)  # [-1. -1. 0. 1. 1.]

# Conditional based on another array
mask = jnp.array([True, False, True, False, True])
result = jnp.where(mask, x, 0.0)
print(result)  # [-3. 0. 0. 0. 5.]

import jax
import jax.numpy as jnp

# IMPERATIVE: modify accumulator in a loop
def imperative_cumsum(arr):
    result = []
    total = 0
    for x in arr:
        total += x
        result.append(total)
    return result

# FUNCTIONAL: use jax.lax.scan (functional loop)
def functional_cumsum(arr):
    def step(carry, x):
        new_carry = carry + x
        return new_carry, new_carry  # (next_carry, output)

    _, cumulative = jax.lax.scan(step, 0.0, arr)
    return cumulative

x = jnp.array([1.0, 2.0, 3.0, 4.0, 5.0])
print(functional_cumsum(x))  # [1. 3. 6. 10. 15.]

💡 .at vs jnp.where

Code

External links

Exercise

Progress

댓글 0