import numpy as np
import galsim