---
title: 'Zero Coordinate Shift (ZCS): Operator Learning & Codes'
url: https://www.emergentmind.com/topics/zero-coordinate-shift-zcs
type: topic
---

# Zero Coordinate Shift (ZCS): Operator Learning & Codes

Zero Coordinate Shift (ZCS) refers to a class of algorithmic and mathematical constructions that, depending on research context, address either (1) optimally structured code sequences for communication with guaranteed zero correlation within specific shift zones, or (2) efficient automatic differentiation for operator learning in physics-informed deep learning. Despite sharing a name, these usages are distinct in scope and technical realization. The following entry details both interpretations as independently developed in the contemporary literature.

## 1. ZCS in Physics-Informed Operator Learning

### Motivation and Background

Zero Coordinate Shift (ZCS) in the context of physics-informed operator learning addresses the scalability bottleneck in automatic differentiation (AD) when computing high-order spatial or temporal derivatives with respect to collocation coordinates in Physics-Informed Neural Operators (PINOs) such as DeepONets. Standard AD frameworks require making all sampled coordinates leaf nodes, resulting in a many-roots-many-leaves (_∞∞) problem; this approach causes computational and memory requirements to scale poorly with the number of input functions $M$ and collocation points $N$ [2311.00860].

Existing attempts to circumvent this include explicit function loops (computing separate graphs for each input function) and data vectorization (treating the entire $(M \times N)$ field as a single vector), each incurring $O(M)$ memory and compute blowup. ZCS introduces a shared infinitesimal shift variable $z$ (per spatial dimension), collapsing these differentiation complexities into substantially smaller computational graphs.

### Mathematical Formulation

Let $u_{ij} = f_\phi(p_i, x_j)$ denote the network output for function $p_i$ at collocation point $x_j$. Standard AD computes derivatives $\partial u_{ij}/\partial x_j$ by making all coordinates $\{x_j\}$ leaves in the graph. ZCS introduces a scalar $z$ so that
$$
v_{ij}(z) = f_\phi(p_i, x_j + z),
$$
which enables the desired derivative to be obtained as
$$
\left. \frac{\partial u_{ij}}{\partial x_j} \right|_{x_j} = \left. \frac{\partial v_{ij}}{\partial z} \right|_{z=0}.
$$
Here, only $z$ is a leaf, shifting all reliance on the coordinate dimension to a single computational node.

To recover the full $(M \times N)$ tensor of derivatives efficiently, one introduces a root-collapsing weight tensor $a_{ij}$ and considers
$$
\omega = \sum_{i=1}^M \sum_{j=1}^N a_{ij} v_{ij}.
$$
The entries $v_{ij}$ are then available via $v_{ij} = \partial \omega / \partial a_{ij}$, and the corresponding derivatives as mixed second derivatives $\partial^2\omega / \partial z\,\partial a_{ij}$. Higher-order derivatives and mixed products follow recursively by this formalism, reducing computation from an _∞∞ problem to a sequence of _11 and _∞1 operations suitable for reverse-mode AD.

### Algorithmic Realization

The procedure, expressed in PyTorch pseudocode for two spatial dimensions $(x, y)$, is as follows:

```python
# Step 1: Create shift scalars and dummy weights
z_x = torch.tensor(0., requires_grad=True)
z_y = torch.tensor(0., requires_grad=True)
a = torch.ones((M, N), requires_grad=True)

# Step 2: Shift coordinates
X = x + z_x
Y = y + z_y

# Step 3: Network forward pass
u = f_theta(p, (X, Y))   # shape M x N

# Step 4: Collapse outputs into a scalar root
omega = (a * u).sum()

# Step 5: First derivatives w.r.t. z
q_x = autograd.grad(omega, z_x, create_graph=True)[0]
q_y = autograd.grad(omega, z_y, create_graph=True)[0]

# Step 6: Second derivatives (Laplacian terms)
s_xx = autograd.grad(q_x, z_x)[0]
s_yy = autograd.grad(q_y, z_y)[0]

# Step 7: Recover M x N field
g = autograd.grad((s_xx + s_yy), a)[0]
return g
```

This method requires only three scalar-to-scalar (_11) and one vector-to-scalar (_∞1) AD operations, eliminating the need for loops or coordinate duplication.

### Implementation Considerations

The ZCS methodology is implemented by subclassing DeepXDE’s operator modules, overriding coordinate handling to use shift scalars and compressing the AD graph. The approach treats $M$ as a batch dimension, with coordinate tensors never duplicated across input functions. The bulk memory and compute thus remain comparable to $M=1$ PINN instances, with the only overhead coming from a small number of scalars and possibly the $a_{ij}$ tensor.

### Empirical Performance

Extensive benchmarking on representative PDE problems (reaction-diffusion, Burgers, Kirchhoff–Love plates, Stokes) demonstrates that ZCS reduces computational graph memory consumption and wall time by an order of magnitude compared to previous approaches. For example, in a Burgers equation with $(M, N, P) = (50, 12,800, 2)$, ZCS achieves 0.20 GB graph memory, 0.36 GB peak memory, and 15 s per 1000 batches, versus 7.8 GB, 7.9 GB, and 316 s for function loop methods [2311.00860]. Relative errors are essentially unchanged between methods.

A summary table is provided:

| Problem (M,N,P)         | Method    | Graph Mem (GB) | Peak Mem (GB) | Time/1000 batches (s) | Rel. Error   |
|-------------------------|-----------|----------------|---------------|-----------------------|--------------|
| Reaction–Diff (50,1000,2) | FuncLoop  | 0.96           | 0.98          | 181                   | 8.3% ±2.0    |
|                         | DataVect  | 0.97           | 1.46          | 24                    | 9.5% ±2.5    |
|                         | ZCS       | **0.02**       | **0.05**      | **10**                | 8.2% ±2.0    |
| Burgers (50,12800,2)    | FuncLoop  | 7.84           | 7.91          | 316                   | 7.5%         |
|                         | DataVect  | 7.73           | 11.40         | 221                   | 7.2%         |
|                         | ZCS       | **0.20**       | **0.36**      | **15**                | 7.1% ±0.5    |
| Kirchhoff–Love (36,10000,4) | FuncLoop  | 77.6           | 77.6          | 4081                  | 27.3%        |
|                         | DataVect  | –              | –             | –                     | –            |
|                         | ZCS       | **2.36**       | **3.30**      | **144**               | 26.9% ±0.5   |
| Stokes (50,5000,2)      | FuncLoop  | 74.4           | 78.6          | 4253                  | 10.3%        |
|                         | DataVect  | –              | –             | –                     | –            |
|                         | ZCS       | **1.99**       | **3.30**      | **147**               | 10.4% ±0.6   |

A plausible implication is that ZCS facilitates the training of operator learning models at previously infeasible scales.

### Scope and Limitations

ZCS is agnostic to the choice of network architecture, PDE type, or data sampling, and imposes no loss in training accuracy. It does not yield advantages for architectures leveraging explicit grid structure (e.g., Fourier Neural Operators or CNN-based PINNs) with optimized finite difference or FFT differentiators, where $O(1)$ stencils remain optimal. ZCS is formulated for reverse-mode AD, and its application to forward-mode scenarios is an area for potential extension [2311.00860].

## 2. ZCS and Zero Correlation Zone (ZCZ) Sequences 

### Definition and Theoretical Foundations

In communication theory, Zero Coordinate Shift (ZCS) is closely related to the zero correlation zone (ZCZ) property of sequences used for multi-carrier code division multiple access (MC-CDMA) and related systems. For two complex-valued sequences $a = (a_0, ..., a_{L-1})$ and $b = (b_0, ..., b_{L-1})$, the aperiodic cross-correlation is defined as
$$
R_{a, b}(\tau) = \sum_{n=0}^{L-1-\tau} a_n b^*_{n+\tau}, \quad \tau = 0, 1, ..., L-1.
$$
ZCZ of width $Z$ means $R_{a,b}(\tau) = 0$ for all nonzero $|\tau| \leq Z$. For a collection of sequences, when auto- and cross-correlations satisfy this property, the set is a Z-complementary code set (ZCCS) [2105.10147].

A $(M,N,L,Z)$–ZCCS is a family $\mathcal{A} = \{A_0, ..., A_{M-1}\}$ of $M$ sequence sets, each $A_i$ consisting of $N$ sequences of length $L$. The ZCCS definition requires
1. $R_{A_i, A_i}(0) = NL,\;\; R_{A_i, A_i}(\tau) = 0$ for $0 < |\tau| < Z$ (sum auto-correlation vanishes outside zero shift).
2. $R_{A_i, A_j}(\tau) = 0$ for all $i \neq j$ and $|\tau| < Z$ (mutual orthogonality over the zone).

The Feng–Fan–Zhou bound asserts $M \leq N\left\lfloor L/Z \right\rfloor$; equality means optimality.

### Construction via Extended Boolean Functions

Optimal ZCCS constructions are achieved using extended Boolean functions (EBF) $f : \mathbb{Z}_q^m \rightarrow \mathbb{Z}_q$:
$$
f(\vec{x}) = a \cdot \sum_{k=1}^{m-v-1} x_{T(k)} x_{T(k+1)} + \sum_{\ell=1}^m \sum_{t=1}^{q-1} c_{\ell, t} x_\ell^t \pmod{q},
$$
with $q \geq 2$, $m \geq 2$, $0 \leq v \leq m$, $T$ a permutation of $\{1, ..., m-v\}$, $a \in \mathbb{Z}_q^*$, $c_{\ell, t} \in \mathbb{Z}_q$.

Define
$$
s_{p, n}(\vec{x}) = f(\vec{x}) + n x_{T(1)} + \sum_{i=1}^{v+1} p_i x_{m-v+i} \pmod{q},
$$
for $p \in \{0, 1, ..., q^{v+1}-1\}$ (in base $q$), $n \in \mathbb{Z}_q$.

Partitioning the $q^{v+1}q$ sequences into $M = q^{v+1}$ sets $S_p = \{s_{p,0}, ..., s_{p,q-1}\}$ yields an optimal $(q^{v+1}, q, q^m, q^{m-v})$–ZCCS: correlation vanishes up to shifts $|τ| < q^{m-v}$ (proof by generalized Golay argument and monomial analysis) [2105.10147].

A concrete illustration: for $q=2$, $m=3$, $v=1$, $f(\vec{x})=x_1x_2$, one recovers a $(4,2,8,4)$–ZCCS by explicit computation.

### Practical Implications

ZCZ sequences constructed this way are critical for asynchronous MC-CDMA, enabling multiple users to transmit signals robustly with delay offsets up to $Z$ without multi-user or multi-path interference. The practical trade-off—$M \leq N\left\lfloor L/Z \right\rfloor$—permits system design flexibility by adjusting sequence length $L$, flock size $N$, and zone width $Z$. The $q$-ary construction generalizes beyond binary codes, supporting large, optimal ZCCS for high-capacity networks [2105.10147].

## 3. Applications in Deep Learning and Communications

In physics-informed deep learning, ZCS enables efficient training of neural operators for PDEs, with empirical evidence indicating scalability to previously unreachable sizes for operator dimensions $M$ and collocation points $N$, without increased training errors [2311.00860].

In wireless communications, ZCS in the ZCCS/ZCZ sense is foundational to MC-CDMA and analogous systems, supporting robust, interference-free multi-user transmission and flexible protocol design [2105.10147].

## 4. Comparative Summary and Terminological Distinctions

ZCS, as developed independently in operator learning and communications code design, addresses orthogonal concerns: backpropagation efficiency in the former, and interference-robust sequence structure in the latter. In operator learning, ZCS is a differentiable shift variable trick; in sequence design, it is a condition on correlation functions. Each is grounded in rigorous mathematical analysis and has led to substantial advances in its respective field.

## 5. Limitations and Prospective Extensions

In physics-informed operator learning, ZCS is not advantageous for grid-based models amenable to $O(1)$ stencil or FFT-based differentiation, and is primarily realized for reverse-mode AD systems. Anticipated extensions include porting to forward-mode AD in frameworks such as JAX or Julia for further performance improvement at high derivative orders.

In code design, ZCS/ZCZ constructions hinge on the algebraic constraints of the employed EBFs, and while highly flexible in $q$-ary settings, remain subject to combinatorial existence bounds (Feng–Fan–Zhou).

## 6. References

- "Zero Coordinate Shift: Whetted Automatic Differentiation for Physics-informed Operator Learning" [2311.00860]
- "New Construction of Z-Complementary Code Sets and Mutually Orthogonal Complementary Sequence Sets" [2105.10147]

Source: https://www.emergentmind.com/topics/zero-coordinate-shift-zcs