Computation of correlation matrices for tetraquark candidates with Jp = 0+ and flavor structure q1q2q3q3