�
^Ǖi� � � � S r S rS r\S:X Ga* SSKrSr\" \\- 5 V s/ s H n \" U S- 5 PM sn r\" \\- 5 V s/ s H n \" U S- 5 PM sn r \R " 5 r
\" S 5 H
r \" \\ \5 M \R " 5 \
-
r \
" S
\ S \ S \ S
S35 \" S5 V s/ s H n \" U 5 PM
sn r\" S5 V s/ s H n \" U 5 PM
sn r\R " 5 r
\" S5 H r \" \\5 M \R " 5 \
-
r \
" S\ S
S35 ggs sn f s sn f s sn f s sn f )zv
Pure Python matrix multiply - target for rustify-ml acceleration.
Expected gain: 50-100x via ndarray matmul in Rust.
c �� � S/X"- - n[ U5 HI n[ U5 H7 nSn[ U5 H nX`XB- U- XU- U- - -
nM XcXB- U- ' M9 MK U$ )u? Multiply two n×n matrices stored as flat Vec<f64> (row-major).� )�range)�a�b�n�result�i�j�total�ks �=d:\WindsurfProjects\rustify\rustify-ml\examples\matrix_ops.py�matmulr st � ��U�a�e�_�F�
�1�X���q��A��E��1�X���1�5�1�9� ��a�%�!�)� �4�4�� � %�1�5�1�9� � � � �M� c �X � Sn[ [ U 5 5 H nX U X - -
nM U$ )z0Compute dot product of two equal-length vectors.r )r �len)r r r r s r
� dot_productr s0 � � �E�
�3�q�6�]��
�1����� �� � �Lr �__main__� N� � � �� zmatmul �xz: z.3fzs for 200 iterationsi� i' zdot_product len=1000: zs for 10000 iterations)�__doc__r r �__name__�timer r �floatr r �start�_�elapsed�print�va�vb)r s 0r
�<module>r$ sW ��� �� �z���
�A� $�Q��U�|�,�|�!��q�1�u� �|�,�A� $�Q��U�|�,�|�!��q�1�u� �|�,�A� �I�I�K�E�
�3�Z���q�!�Q�� ��i�i�k�E�!�G� �G�A�3�a��s�"�W�S�M�)=�
>�?� !�$�K� (�K�q�%��(�K� (�B� !�$�K� (�K�q�%��(�K� (�B� �I�I�K�E�
�5�\���B��� ��i�i�k�E�!�G� � "�7�3�-�/E�
F�G�% �� -��,��
)�� (s �D;�E �;E�E