Tensor method for optimal control problems constrained by fractional three-dimensional elliptic operator with variable coefficients