I have a map $\mathcal{N}$ from the space of two-qubit subnormalised density matrices $\mathcal{S}(\mathcal{H}_2 \otimes \mathcal{H}_2)$ to itself (positive operators with trace between 0 and 1). Notice that it is not necessarily a physical one (e.g. it may be nonlinear, that's perfectly fine). I know $\mathcal{N}$ is
linear for all product states $\rho_A \otimes \rho_B$. By this I mean $\mathcal{N} \left( \sum_{ijkl} c_{ik}d_{jl} |i\rangle \langle k| \otimes |j \rangle \langle l | \right)= \sum_{ijkl} c_{ik}d_{jl} \mathcal{N} \left( |i\rangle \langle k| \otimes |j \rangle \langle l | \right)$ for every $\rho_A = \sum_{ik}c_{ik}|i\rangle \langle k|, \rho_B = \sum_{jl}d_{jl}|j \rangle \langle l| \in \mathcal{S}(\mathcal{H}_2)$
continuous on the whole space of two-qubit density matrices $\mathcal{S}(\mathcal{H}_2 \otimes \mathcal{H}_2)$.
Does it automatically follow that $\mathcal{N}$ has to be linear over all $\mathcal{S}(\mathcal{H}_2 \otimes \mathcal{H}_2)$?