The first passage time (FPT) problem is an important problem with a wide range of applications in science, engineering, economics, and industry. Mathematically, such a problem can be reduced to estimating the probability of a stochastic process first to reach a boundary level. In most important applications in the financial industry, the FPT problem does not have an analytical solution and the development of efficient numerical methods becomes the only practical avenue for its solution. Most of our examples in this contribution are centered around the evaluation of default correlations in credit risk analysis, where we are concerned with the joint defaults of several correlated firms, the task that is reducible to a FPT problem. This task represents a great challenge for jump-diffusion processes (JDP). In this contribution, we develop further our previous fast Monte Carlo method in the case of multivariate (and correlated) JDP. This generalization allows us, among other things, to evaluate the default events of several correlated assets based on a set of empirical data. The developed technique is an efficient tool for a number of financial, economic, and business applications, such as credit analysis, barrier option pricing, macroeconomic dynamics, and the evaluation of risk, as well as for a number of other areas of applications in science and engineering, where the FPT problem arises.Many applications in science, engineering, finance, economics, and industry require the solution of the first passage problem. Examples range from population genetics and other applications in biology, to neurobiology and neurophysiology, to molecule kinetics and problems in chemistry and physics, to rare events modelling and quantum information, and to aerospace engineering applications (e.g. [1][2][3][4]). This problem also arises in financial applications. Although methodologies developed in this paper are applicable to a number of other areas, we exemplify our discussion here with problems from finance.In the financial world, individual companies are usually linked together via economic conditions, so default correlation, defined as the risk of multiple companies' default together, has been an important area of research in credit analysis with applications to joint defaults, credit derivatives, asset pricing, and risk management.Currently, there are two dominant groups of theoretical models used in default correlation. One is a reduced form model, such as in [5] that uses a Copula function to parameterize the default correlation. In this context, it is worthwhile noting that Chen and Sopranzetti [6] have translated the joint default probability into a bivariate normal probability function, making the analysis and applications of such models potentially more convenient.The second group of models for default correlation is a structural form model. Zhou [7] and Hull and White [8] were the first to incorporate default correlation into the Black-Cox first passage structural model. They obtained similar closed-form solutions for two ...