Materials with strong electronic correlations often exhibit a superconducting phase in close competition with other insulating phases, which is outstandingly difficult to resolve, e.g., for a group of quasi-two-dimensional (Q2D) materials such as the cuprates, even for the simplified minimal model of these materials, the doped 2D Hubbard model. The present work shows how quasi-one-dimensional (Q1D) systems, 2D and three-dimensional (3D) arrays of weakly coupled 1D correlated electrons, are much more amenable to resolving such competition, treating both instabilities on equal footing. Using the recently established matrix product state plus mean field (MPS+MF) approach for fermions [Bollmark et al., Phys. Rev. X 13, 011039 (2023)], we demonstrate that large systems can be reached readily in these systems, which opens up the thermodynamic regime via extrapolation. Focusing on basic model systems, 3D arrays of negative-U Hubbard chains with additional nearest-neighbor interaction V, we show that despite the MF component of the MPS+MF technique, we can reproduce the expected coexistence of the superconductivity and charge density wave at V = 0 for density n = 1. We then show how we can tune away from coexistence by both tuning V and doping the system. This work thus paves the way to deploy two-channel MPS+MF theory on some highly demanding high-Tc superconducting systems, such as 3D arrays of repulsive-U doped Hubbard ladders; we recently characterized the properties of such arrays in single-channel MPS+MF calculations [Bollmark et al., Phys. Rev. X 13, 011039 (2023)].