A method of inducing a periodic variation of nonlinearity value in a sample of ferroelectric material comprises arranging a pair of electrodes on opposite faces of the sample, one electrode defining a desired pattern of nonlinearity variation, applying a pre-bias voltage across the sample for a predetermined time using the electrodes, the pre-bias voltage being less than the coercive field of the ferroelectric material; and after the predetermined time, applying a current-controlled poling voltage across the sample using the electrodes, to produce domain inversion in the sample according to the desired pattern of nonlinearity variation. The pre-bias voltage may be 75% of the coercive field or more, and applied for a pre-determined time between 1 and 100 seconds.