A suite of programs for the numerical simulation of two wave mixing pr
ocesses in sillenite crystals is presented, The interacting waves are
assumed to possess arbitrary amplitude and/or phase profiles in two sp
atial dimensions as well as arbitrary polarization states, Several doc
umented program tests demonstrate the capabilities of the suite, inclu
ding the handling of waves slowly modulated in time.