An upper semi-continuous function $u : \Omega \to \mathbb{R}$, $\Omega \subseteq \mathbb{C}^n$ is said to be *subharmonic* if it satisfies the submean inequality $u(a) \leq \mu\_S(u;a,r)$, where $\mu\_S(u;a,r)$ is the mean value of $u$ on the sphere $S(a,r)$ with center $a$ and radius $r$ (for any $a$,$r$ such that the closed ball $\overline{B(a,r)}$ is contained in $\Omega$). 


The function is said to be *plurisubharmonic* if it is subharmonic restricted to any complex line. It is straight-forward to see that plurisubharmonic implies subharmonic, see for example Demailly, "Complex Differential and Algebraic Geometry", §I.5.A. It is also not so difficult to see that plurisubharmonicity is invariant under holomorphic changes of coordinates, see for example Demailly, Theorem I.5.11.


The entry [Encyclopedia of Mathematics: Pluripotential theory](https://www.encyclopediaofmath.org/index.php/Pluripotential_theory) states that "Plurisubharmonic functions are precisely the subharmonic functions invariant under a holomorphic change of coordinates". One inclusion is clear from above, but does anyone have an idea of how to prove the other inclusion, i.e.: if $u$ is subharmonic with respect to any local holomorphic coordinates, then $u$ is plurisubharmonic?



Subharmonic in any holomorphic coordinates = Plurisubharmonic?
It is actually enough to assume that your function $u$ defined on a neighborhood of $0\in \mathbb C^n$ remains subharmonic after composing with any linear transformation. 


Indeed, let $0\neq \lambda \in \mathbb C^n$ and set, for $\sigma \in \mathbb C^\*$,


$$A(\sigma) = \begin{pmatrix}
\lambda\_1 & \sigma a\_{12} & \cdots & \sigma a\_{1n}\\\
\vdots & \vdots & & \vdots \\\
\vdots & \vdots & &\vdots\\\
\lambda\_n & \sigma a\_{n2} & \cdots & \sigma a\_{nn}
\end{pmatrix}$$


with $a\_{ij} $ chosen such that the matrix is invertible. Then one finds, at the origin


$$\begin{eqnarray}
\Delta (u \circ A(\sigma))&=&\sum\_{k=1}^n \frac{\partial^2(u \circ A(\sigma))}{\partial z\_k \partial \bar{z\_k}}(0)\\\
&=&\sum\_{k=1}^n \sum\_{1\leq i,j\leq n}A(\sigma)\_{ik} \overline{A(\sigma)\_{jk}} \frac{\partial^2 u}{\partial z\_i \partial \bar {z\_j}}(0)\\\
&=&Hu\_{0}(\lambda) + |\sigma|^2 \sum\_{\ell=2}^n Hu\_{0}((a\_{j \ell})\_{j})
\end{eqnarray}$$


where $H$ denotes the complex hessian. This quantity is non-negative by assumption, and one gets the result by making $\sigma$ approach $0$.


