static void wrap() 
        {
            bp::class_< MultipleImageHelper > pyMultipleImageHelper(
                "MultipleImageHelper", bp::init<const MultipleImageHelper &>()
            );
            wrapTemplates_Multi<float>(pyMultipleImageHelper);
            wrapTemplates_Multi<double>(pyMultipleImageHelper);
            wrapTemplates_Multi<int32_t>(pyMultipleImageHelper);
            wrapTemplates_Multi<int16_t>(pyMultipleImageHelper);

            bp::class_< SBInterpolatedImage, bp::bases<SBProfile> > pySBInterpolatedImage(
                "SBInterpolatedImage", bp::init<const SBInterpolatedImage &>()
            );
            pySBInterpolatedImage
                .def(bp::init<const MultipleImageHelper&, const std::vector<double>&,
                     boost::shared_ptr<InterpolantXY>,
                     boost::shared_ptr<InterpolantXY>,
                     boost::shared_ptr<GSParams> >(
                         (bp::args("multi","weights"),
                          bp::arg("xInterp"),
                          bp::arg("kInterp"),
                          bp::arg("gsparams")=bp::object())
                     )
                )
                .def("calculateStepK", &SBInterpolatedImage::calculateStepK,
                     bp::arg("max_stepk")=0.)
                .def("calculateMaxK", &SBInterpolatedImage::calculateMaxK,
                     bp::arg("max_maxk")=0.)
                ;
            wrapTemplates<float>(pySBInterpolatedImage);
            wrapTemplates<double>(pySBInterpolatedImage);
            wrapTemplates<int32_t>(pySBInterpolatedImage);
            wrapTemplates<int16_t>(pySBInterpolatedImage);
        }
    void pyExportSBInterpolatedImage(PY_MODULE& _galsim)
    {
        py::class_<SBInterpolatedImage, BP_BASES(SBProfile)> pySBInterpolatedImage(
            GALSIM_COMMA "SBInterpolatedImage" BP_NOINIT);
        pySBInterpolatedImage
            .def("calculateMaxK", &SBInterpolatedImage::calculateMaxK);
        WrapTemplates<float>(_galsim, pySBInterpolatedImage);
        WrapTemplates<double>(_galsim, pySBInterpolatedImage);

        py::class_<SBInterpolatedKImage, BP_BASES(SBProfile)> pySBInterpolatedKImage(
            GALSIM_COMMA "SBInterpolatedKImage" BP_NOINIT);
        pySBInterpolatedKImage
            .def(py::init<const BaseImage<std::complex<double> > &,
                 double, const Interpolant&, GSParams>());
    }