def METHOD_NAME(self, locs):
    """Compute stress field at locations.
    """
    (npts, dim) = locs.shape
    METHOD_NAME = numpy.zeros((numSteps, npts, self.TENSOR_SIZE), dtype=numpy.float64)
    METHOD_NAME[:,:, 0] = sxx.reshape(numSteps, 1)
    METHOD_NAME[:,:, 1] = syy.reshape(numSteps, 1)
    METHOD_NAME[:,:, 2] = szz.reshape(numSteps, 1)
    METHOD_NAME[:,:, 3] = sxy.reshape(numSteps, 1)
    return METHOD_NAME