This paper explains the methodology used to develop a high-resolution, multi-dimensional Euler solver that is capable of handling non-ideal equation of state and stiff chemical source terms. We have developed a pointwise implementation that has computational advantages for our intended applications, as opposed to a finite volume implementation. Our solver allows for the placement of internal reflective boundaries and the standard inflow and outflow and reflective boundaries at the edge of the domain. We discuss the spatial discretization and the temporal integration schemes, upwinding and flux splitting and the combined use of the Lax - Friedrichs and Roe schemes to solve for the required fluxes. A complete description of the pointwise internal boundary method is given. An overall summary of a representative code structure is given. We provide details on the verification of our integrated set of algorithms that resulted in an application code. We demonstrate the order of convergence for test problems. Two example applications from measurement of detonation shock dynamics and deflagration to detonation transition in porous energetic materials are presented.