order plane-wave incidence is proposed as an efficient calculation approach to simulate
three-dimensional photonic crystal devices. As an example, the dispersion relations and
quality factors are calculated for resonant cavity arrays embedded in a woodpile photonic
crystal. An interesting ultraslow negative group velocity is observed in this structure.