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())
Exemple #2
0
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)