Massive stars end their lives through core collapse. Some explode as core-collapse supernovae (CCSN), leaving behind neutron stars (NS), while others collapse to stellar-mass black holes (BH). Determining the conditions under which a massive star will explode remains a significant challenge in CCSN theory. I will present recent advancements in analytic theory and numerical simulations focusing on predicting the outcomes of massive star evolution. The fate of massive stars depends on an (im)balance of various forces, including the effects of neutrino heating and the ram pressure of accreting material. We have developed an analytic force explosion condition (FEC) that provides a criterion for successful explosions based on four dimensionless parameters: net neutrino heating in the gain region, neutrino opacity, buoyant driving, and the radial component of Reynolds stress. I will highlight our progress in validating the FEC through multi-dimensional simulations. Additionally, I will discuss our efforts in developing next-generation CCSN simulations designed to improve computational speed and efficiency. These advancements will enhance our understanding of CCSN explosion mechanism and the resulting outcomes such as the mass distributions of neutron stars and black holes. By integrating analytical and numerical approaches, we aim to refine predictions regarding the fate of massive stars and the properties of the compact objects they leave behind.