Symmetry-unrestricted Hartree-Fock-Bogoliubov calculations with finite-range density-dependent Gogny force are presented with a triaxial harmonic oscillator basis. We include all the contributions to the paring and Hartree-Fock potentials arising from Gogny and Coulomb interaction as well as the center of mass correction. Ground-state properties are calculated for C, O and Ne isotopes, and compared with available experiment, Hartree-Fock plus BCS, shell model and Relativistic-Hartree-Bogoliubov calculations. We show that predicted drip lines are strongly model dependent due to their sensitivity to various theoretical details. Shape coexistence phenomenon is discussed. ©2001 American Institute of Physics.