📄 matrix.c.htm
字号:
<?xml version="1.0" encoding="iso-8859-1"?><!DOCTYPE html PUBLIC "-//W3C//DTD XHTML 1.1//EN" "http://www.w3.org/TR/xhtml11/DTD/xhtml11.dtd"><html xmlns="http://www.w3.org/1999/xhtml"><head> <title>matrix.c</title> <link rel="stylesheet" type="text/css" href="highlight.css"/></head><body><pre><span class="line"> 1 </span><span class="dir">#include <stdlib.h></span><span class="line"> 2 </span><span class="dir"></span><span class="dir">#include <stdio.h></span><span class="line"> 3 </span><span class="dir"></span><span class="line"> 4 </span><span class="dir">#include</span> <span class="dstr">"matrix.h"</span><span class="dir"></span><span class="line"> 5 </span><span class="dir"></span><span class="line"> 6 </span><span class="kwb">void</span> matrix_mult_f<span class="sym">(</span><span class="kwb">float</span> <span class="sym">*</span>m1<span class="sym">,</span> <span class="kwb">unsigned int</span> m1_rows<span class="sym">,</span> <span class="kwb">unsigned int</span> m1_cols<span class="sym">,</span><span class="line"> 7 </span><span class="sym"></span> <span class="kwb">float</span> <span class="sym">*</span>m2<span class="sym">,</span> <span class="kwb">unsigned int</span> m2_rows<span class="sym">,</span> <span class="kwb">unsigned int</span> m2_cols<span class="sym">,</span><span class="line"> 8 </span><span class="sym"></span> <span class="kwb">float</span> <span class="sym">*</span>ans<span class="sym">)</span><span class="line"> 9 </span><span class="sym">{</span><span class="line"> 10 </span><span class="sym"></span> <span class="slc">// position markers in the output matrix</span><span class="line"> 11 </span><span class="slc"></span> <span class="kwb">long int</span> ans_row<span class="sym">,</span> ans_col<span class="sym">;</span><span class="line"> 12 </span><span class="sym"></span><span class="line"> 13 </span><span class="sym"></span><span class="line"> 14 </span><span class="sym"></span> <span class="slc">// check the matrix dimensions match</span><span class="line"> 15 </span><span class="slc"></span> <span class="kwa">if</span> <span class="sym">(</span>m1_cols <span class="sym">!=</span> m2_rows<span class="sym">)</span><span class="line"> 16 </span><span class="sym"> {</span><span class="line"> 17 </span><span class="sym"></span> printf<span class="sym">(</span><span class="str">"ERROR:</span><span class="esc">\n </span> <span class="str">matrix_mult_f: Matrix dimensions do not match, m1_cols and m2_rows must be the same. %d != %d</span><span class="esc">\n</span><span class="str">"</span><span class="sym">,</span><span class="line"> 18 </span><span class="sym"></span> m1_cols<span class="sym">,</span> m2_rows<span class="sym">);</span><span class="line"> 19 </span><span class="sym"></span><span class="line"> 20 </span><span class="sym"></span> exit<span class="sym">(-</span><span class="num">1</span><span class="sym">);</span><span class="line"> 21 </span><span class="sym"> }</span><span class="line"> 22 </span><span class="sym"></span><span class="line"> 23 </span><span class="sym"></span><span class="line"> 24 </span><span class="sym"></span> <span class="kwa">for</span> <span class="sym">(</span> ans_row<span class="sym">=</span><span class="num">0</span><span class="sym">;</span> ans_row <span class="sym"><</span> m1_rows<span class="sym">;</span> ans_row<span class="sym">++ )</span><span class="line"> 25 </span><span class="sym"> {</span><span class="line"> 26 </span><span class="sym"></span> <span class="kwa">for</span> <span class="sym">(</span> ans_col<span class="sym">=</span><span class="num">0</span><span class="sym">;</span> ans_col <span class="sym"><</span> m2_cols<span class="sym">;</span> ans_col<span class="sym">++)</span><span class="line"> 27 </span><span class="sym"> {</span><span class="line"> 28 </span><span class="sym"></span> <span class="slc">// loop variable for next 'for' loop</span><span class="line"> 29 </span><span class="slc"></span> <span class="kwb">int</span> i<span class="sym">;</span><span class="line"> 30 </span><span class="sym"></span> <span class="kwb">float</span> <span class="sym">*</span>ans_addr <span class="sym">=</span> ans <span class="sym">+ (</span>ans_row <span class="sym">*</span> m2_cols<span class="sym">) +</span> ans_col<span class="sym">;</span><span class="line"> 31 </span><span class="sym"></span><span class="line"> 32 </span><span class="sym"> *(</span>ans_addr<span class="sym">) =</span> <span class="num">0</span><span class="sym">;</span> <span class="slc">//initialise the answer value</span><span class="line"> 33 </span><span class="slc"></span><span class="line"> 34 </span> <span class="kwa">for</span><span class="sym">(</span> i <span class="sym">=</span> <span class="num">0</span><span class="sym">;</span> i <span class="sym"><</span> m1_cols<span class="sym">;</span> i<span class="sym">++)</span><span class="line"> 35 </span><span class="sym"> {</span><span class="line"> 36 </span><span class="sym"></span> <span class="kwb">float</span> m1_val<span class="sym">,</span>m2_val<span class="sym">;</span><span class="line"> 37 </span><span class="sym"></span><span class="line"> 38 </span><span class="sym"></span> m1_val <span class="sym">= *(</span>m1 <span class="sym">+ (</span>ans_row <span class="sym">*</span> m1_cols <span class="sym">+</span> i<span class="sym">));</span><span class="line"> 39 </span><span class="sym"></span> m2_val <span class="sym">= *(</span>m2 <span class="sym">+ (</span>i <span class="sym">*</span> m2_cols <span class="sym">+</span> ans_col<span class="sym">));</span><span class="line"> 40 </span><span class="sym"> *(</span>ans_addr<span class="sym">) +=</span> m1_val <span class="sym">*</span> m2_val<span class="sym">;</span><span class="line"> 41 </span><span class="sym"> }</span><span class="line"> 42 </span><span class="sym"></span><span class="line"> 43 </span><span class="sym"> }</span><span class="line"> 44 </span><span class="sym"> }</span><span class="line"> 45 </span><span class="sym"></span> <span class="kwa">return</span><span class="sym">;</span><span class="line"> 46 </span><span class="sym">}</span><span class="line"> 47 </span><span class="sym"></span><span class="line"> 48 </span><span class="sym"></span><span class="line"> 49 </span><span class="sym"></span><span class="kwb">void</span> comp_matrix_mult_f<span class="sym">(</span><span class="kwb">float</span> <span class="sym">*</span>m1_re<span class="sym">,</span> <span class="kwb">float</span> <span class="sym">*</span>m1_im<span class="sym">,</span><span class="line"> 50 </span><span class="sym"></span> <span class="kwb">unsigned int</span> m1_rows<span class="sym">,</span> <span class="kwb">unsigned int</span> m1_cols<span class="sym">,</span><span class="line"> 51 </span><span class="sym"></span> <span class="kwb">float</span> <span class="sym">*</span>m2_re<span class="sym">,</span> <span class="kwb">float</span> <span class="sym">*</span>m2_im<span class="sym">,</span><span class="line"> 52 </span><span class="sym"></span> <span class="kwb">unsigned int</span> m2_rows<span class="sym">,</span> <span class="kwb">unsigned int</span> m2_cols<span class="sym">,⌨️ 快捷键说明
复制代码
Ctrl + C
搜索代码
Ctrl + F
全屏模式
F11
切换主题
Ctrl + Shift + D
显示快捷键
?
增大字号
Ctrl + =
减小字号
Ctrl + -