Model order reduction of discrete-time parameter-varying systems via bilinearization and low-rank Gramian approximation | Synapse