a ›¬óz#test_matrix_basic..)réé)é(ré<)éFéPéZyð?csˆˆSrrr)rÚA4rrr$rrr r r éöÿÿÿ)r#r#r#r#csˆdS)Nrrr)ÚA7rrr-riœÿÿÿ©rr)ÚmatrixÚrangeÚeyeÚ _matrix__dataÚlistÚpytestÚraisesÚ ValueErrorÚtolistÚevalÚreprÚfpÚrowsÚcolsÚlenÚsumÚonesÚzerosÚ randmatrixZnstr) ÚA1ÚiZA3ÚA5ÚlZA6ÚxZA8ZA9ZA10r)rr"r$rÚtest_matrix_basicsZ    & &  ""r>cCsNtjdkrt d¡tgd¢gd¢gƒ}tddgddgd d ggƒ}td ƒd S) z° Test the PEP465 "@" matrix multiplication syntax. To avoid syntax errors when importing this file in Python 3.4 and below, we have to use exec() - sorry for that. iðz9'@' (__matmul__) is only supported in Python 3.5 or newerrrr rrrrrzassert A4 @ A5 == A4 * A5N)ÚsysÚ hexversionr+Úskipr&Úexec)r"r;rrrÚ test_matmul4s   rCcstgd¢gd¢gd¢gƒ‰tgd¢ƒ}ˆdd…dd…fˆksBJ‚ˆdd…dftdgdgdggƒksjJ‚ˆddd…ftgd¢gƒksŒJ‚ˆdd …dd …ftdd gdd ggƒks¸J‚|dd …td d gƒksÔJ‚t t‡fd d„¡td ƒ}ˆ|dd…dd…f<|dd…dd…ftgd¢gd¢gd¢gƒks8J‚tgd¢gƒ|ddd…f<|tgd¢gd¢gd¢gƒksvJ‚tdgdgdggƒ|dd…df<|tgd¢gd¢gd¢gƒksºJ‚tddgddggƒ|dd…dd…f<|tgd¢gd¢gd¢gƒksJ‚d|dd …<|tgd¢ƒks$J‚t t¡,ˆdd…df|ddd…f<Wdƒn1sb0Yt t¡0ˆdd…dd…f|ddd…f<Wdƒn1s®0Yd|dd…df<|tgd ¢gd!¢gd"¢gƒksìJ‚d#|dd…dd…f<|D]}|d#ksJ‚qdS)$Nrrr )rrrrr rrr r rr rrcsˆdd…dd…fS)Nrr rr©ÚArrrMrz$test_matrix_slices..)ré rré éé)rrFrG)rr rH)r r rIéééé)rJrKrG)rLrMrHr)rrrrr r)rJrKr)rLrMr)r r rr)r&r+r,Ú IndexErrorr-)ÚVr9r=rrDrÚtest_matrix_slicesAs\þ (",þ þ  þ$ þ  < @ þrPcCsptddgddggƒ}|d||ks(J‚|d|||ks@J‚|dt|ƒksTJ‚|dt||ƒkslJ‚dS)Nrrrrréþÿÿÿ)r&ZinverserDrrrÚtest_matrix_powerss rRcCsÐtddgddgddggƒ}|j| ¡krDtgd¢gd¢gƒksJnJ‚t|ddƒ|tddgddgddggƒksxJ‚ddg}t|d dƒ|ddgksœJ‚ttdƒgd ¢ƒtgd ¢gd ¢gd ¢gƒksÌJ‚dS)Nrrrrr r )rrr )rrr rr)rrrr)rrrr)rrrr)r&ÚTZ transposeZswap_rowÚextendr()rEr<rrrÚtest_matrix_transformzs0 " rUcCs„tdtdgdtggƒ}| ¡ttddƒdgdtddƒggƒksDJ‚| ¡|jkrzttddƒdgdtddƒggƒks€nJ‚dS)Nrrrr)r&ÚjÚ conjugateZmpcZtranspose_conjÚHrDrrrÚtest_matrix_conjugate„s ,  ÿrYcCstgd¢ƒtgd¢gd¢gd¢gƒks*J‚tddƒ}|jdkrH|jdksLJ‚|D]}|dksPJ‚qPtddƒ}|jdkr€|jdks„J‚|D]}|dksˆJ‚qˆtd ƒtd ƒks®J‚tdƒ}tdƒt||d|dg|d|d|d g|d|d |d ggƒksJ‚dS) Nr)rrr)rrr)rrrrrrrrrr ) Zdiagr&r6r2r3r7r8ÚmpfZhilbert)r9ÚarZonerrrÚtest_matrix_creationŠs*  þr\cCsÄtddgddgddggƒ}t|dƒdks,J‚t|tƒdks>J‚t|dƒtd ƒksTJ‚tdƒd ksdJ‚gd ¢}t|dƒd ks~J‚tt|dƒd ƒdks–J‚tt|d ƒd ƒdks®J‚t|tƒdksÀJ‚dS)NrrQrrrr rÚFrr)rrQr rérgPÃþ…z$,@g %^Ì(@r)r&ZmnormÚinfÚsqrtZnormÚround)rEr=rrrÚ test_normsšsrbcCsÌtgd¢ƒ}|tdgdgdgdgdggƒks0J‚|ddks@J‚t|jƒdksRJ‚t|ƒttdƒƒksjJ‚d|d<d|d<|ddksŠJ‚t|ƒt|jƒkr¨dks®nJ‚|j|td ggƒksÈJ‚dS) N)rrrrrrrrrrr r#ér)r&r4r)r*r'rS)r=rrrÚ test_vector¨s $$rdcCsZtdƒ}| ¡}| }||ks"J‚||ks.J‚d|d<||ksBJ‚d|d<||ksVJ‚dS)Nr rr%é*)r6Úcopy)rEÚBÚCrrrÚtest_matrix_copy´s   ricCsZz ddl}Wnty YdS0ddgddgddgg}| |¡}t|ƒt|ƒksVJ‚dS)Nrrrrrr r )ÚnumpyÚ ImportErrorÚarrayr&)rjr<r[rrrÚtest_matrix_numpy¿s   rmc Cs¦tddƒ}||d}tdƒ}||d}d}t |¡}||d}t |¡}t d¡}||||||||fD]4} | |t | g¡ksˆJ‚|| t | g¡kslJ‚qldS)z/Multiplication of iv.matrix and any scalar typerry@regX9´Èv¾ó?y@N)ÚmpirZr1ÚconvertÚivr6r&) r[ÚbÚcÚdÚeÚfÚgÚhÚMr=rrrÚ test_interval_matrix_scalar_multÈs       rycCs|tdƒ}t d¡}t d¡}|||fD]P}||t |¡ks@J‚|||ksPJ‚||t |¡ksfJ‚|||ks&J‚q&dS)z2Multiplication of iv.matrix and other matrix typesrN)r6r1rpr&)rErgrxÚXrrrÚ test_interval_matrix_matrix_mult×s  r{cCshtdƒt d¡t d¡fD]H}t |¡}t d¡}t|dƒt|dƒksNJ‚| ¡| ¡ksJ‚qdS)Nrr%)r(r1rpr&Útyper.)Zother_type_eyerErgrrrÚtest_matrix_conversion_to_ivãs   r}cCs$tdƒ}t|ggƒ}t |¡}t |gg¡}||ks8J‚||}||}||ksTJ‚|djdksfJ‚|djdksxJ‚|djdksŠJ‚|djdksœJ‚t d¡|dvs²J‚t d¡|dvsÈJ‚t t d¡¡t d¡tddƒt tddƒtddƒgtddƒtddƒgg¡ks J‚dS) Nz1.00000000000001r%g¼‰Ø—²ÒœrCrPrRrUrYr\rbrdrirmryÚmarkZxfailr{r}r€rrrrÚs$/ 2