mddat = numpy.array([[3.4, 4.5], [5.6, 6.7], [7.8, 8.9]], dtype=numpy.float64) mcfdat = numpy.array(numpy.vectorize(complex)(mfdat, 8 + mfdat), dtype=numpy.complex64) mcddat = numpy.array(numpy.vectorize(complex)(mddat, 16 + mddat), dtype=numpy.complex128) mf = lal.gsl_matrix_float(mfdat.shape[0], mfdat.shape[1]) mf.data = mfdat mfout = lal.gsl_matrix_float(mfdat.shape[0], mfdat.shape[1]) mfout.data = numpy.zeros(numpy.shape(mfdat), dtype=mfdat.dtype) assert (lal.swig_lal_test_viewin_gsl_matrix_float(mfout, mf)) assert ((mfout.data == mf.data).all()) mfout.data = numpy.zeros(numpy.shape(mfdat), dtype=mfdat.dtype) assert (lal.swig_lal_test_viewin_gsl_matrix_float(mfout, mfdat)) assert ((mfout.data == mfdat).all()) mfout.data = numpy.zeros(numpy.shape(mfdat), dtype=mfdat.dtype) assert (lal.swig_lal_test_viewinout_gsl_matrix_float(mfout, mf)) assert ((2 * mfout.data == mf.data).all()) mfout.data = numpy.zeros(numpy.shape(mfdat), dtype=mfdat.dtype) assert (lal.swig_lal_test_viewinout_gsl_matrix_float(mfout, mfdat)) assert ((2 * mfout.data == mfdat).all()) mf.data = mfdat assert (lal.swig_lal_test_copyinout_gsl_matrix_float(mf)) assert ((mf.data == 3 * mfdat).all()) mf.data = mfdat retn, mf = lal.swig_lal_test_copyinout_gsl_matrix_float(mf) assert (retn) assert ((mf.data == 3 * mfdat).all()) mf = mfdat retn, mf = lal.swig_lal_test_copyinout_gsl_matrix_float(mf) assert (retn) assert ((mf == 3 * mfdat).all())
lal.CheckMemoryLeaks() mfdat = numpy.array([[1.2, 2.3, 3.4], [4.5, 5.6, 6.7]], dtype=numpy.float32) mddat = numpy.array([[3.4, 4.5], [5.6, 6.7], [7.8, 8.9]], dtype=numpy.float64) mcfdat = numpy.array(numpy.vectorize(complex)(mfdat, 8 + mfdat), dtype=numpy.complex64) mcddat = numpy.array(numpy.vectorize(complex)(mddat, 16 + mddat), dtype=numpy.complex128) mf = lal.gsl_matrix_float(mfdat.shape[0], mfdat.shape[1]) mf.data = mfdat mfout = lal.gsl_matrix_float(mfdat.shape[0], mfdat.shape[1]) mfout.data = numpy.zeros(numpy.shape(mfdat), dtype=mfdat.dtype) assert(lal.swig_lal_test_viewin_gsl_matrix_float(mfout, mf)) assert((mfout.data == mf.data).all()) mfout.data = numpy.zeros(numpy.shape(mfdat), dtype=mfdat.dtype) assert(lal.swig_lal_test_viewin_gsl_matrix_float(mfout, mfdat)) assert((mfout.data == mfdat).all()) mfout.data = numpy.zeros(numpy.shape(mfdat), dtype=mfdat.dtype) assert(lal.swig_lal_test_viewinout_gsl_matrix_float(mfout, mf)) assert((2 * mfout.data == mf.data).all()) mfout.data = numpy.zeros(numpy.shape(mfdat), dtype=mfdat.dtype) assert(lal.swig_lal_test_viewinout_gsl_matrix_float(mfout, mfdat)) assert((2 * mfout.data == mfdat).all()) del mf del mfout del mfdat lal.CheckMemoryLeaks() md = lal.gsl_matrix(mddat.shape[0], mddat.shape[1]) md.data = mddat mdout = lal.gsl_matrix(mddat.shape[0], mddat.shape[1]) mdout.data = numpy.zeros(numpy.shape(mddat), dtype=mddat.dtype) assert(lal.swig_lal_test_viewin_gsl_matrix(mdout, md)) assert((mdout.data == md.data).all()) mdout.data = numpy.zeros(numpy.shape(mddat), dtype=mddat.dtype)