The algorithm of the original Reflectivity Method has been vectorized and implemented on a CDC CYBER 205 computer. Calculation times are shortened by a factor of 20 to 30 compared with a general purpose computer with a capacity of several million floating point operations per second (MFLOP). The rapid calculation of synthetic seismograms for complex models, high frequency sources and all offset ranges is a provision for modeling not only particular phases but the whole observed wavefield. As an example we model refraction data of the Black Forest, Southwest Germany and are able to derive rather tight constraints on the physical properties of the lower crust.