I posted my question at MS but unfortunately it is still without a response, so let me ask it here.

We can think about a bounded operator $T\colon c_0\to c_0$ as a double-infinite matrix $[T_{mn}]_{m,n\geq 1}$ which acts on a sequence $a=[a_1, a_2, a_3, \ldots ]\in c_0$ in the same way as usual (finite) matrices act on vectors ($n$-tuples of scalars), i.e.

$$ Ta=  [T_{mn}][a_n] = \left[ \sum_{n=1}^\infty T_{mn}a_n\right] $$

Suppose $a=[a_1, a_2, a_3, \ldots ]\in \ell^\infty = (c_0)^{**}$. Does the following formula still hold:

$$ T^{**}a=  [T_{mn}][a_n] = \left[ \sum_{n=1}^\infty T_{mn}a_n\right] $$

Are there any methods of computing the values of second conjugate operators acting on biduals $X^{**}$ of classical (non-reflexive) spaces $X$ at elements from $X$?