专栏首页数值分析与有限元编程Fortran如何实现矩阵与向量的乘法运算

Fortran如何实现矩阵与向量的乘法运算

矩阵是二维数组,而向量是一维数组,内置函数matmul不能实现矩阵与向量的乘法运算。在这一点Fortran不如matlab灵活。

Fortran如何实现矩阵与向量的乘法运算,现有以下三种方法供参考。

一)将一维数组看作二维数组的退化形式,比如a(3)可以看作a(3,1)或者a(1,3),这样就可以用matmul函数计算了。

二)用spread函数将一维数组扩展成二维数组,同样可用matmul函数计算。

来看过程。

数组c的第一列就是需要的计算结果。

spread(B,2,2)就是按列扩展,成为二维数组

三)利用dot_product函数。dot_product函数是向量点积运算函数,可将二维数组的每一行抽取出来,和一维数组作dot_product运算。

程序员为什么会重复造轮子?现在的软件发展趋势,越来越多的基础服务能够“开箱即用”、“拿来用就好”,越来越多的新软件可以通过组合已有类库、服务以搭积木的方式完成。这是趋势,将来不懂开发语言的人都可以通过利用现有软件组件快速构建出能解决实际问题的软件产品。对程序员来讲,在一开始的学习成长阶段,造轮子则具有特殊的学习意义,学习别人怎么造,了解内部机理,自己造造看,这是非常好的锻炼。每次学习新技术都可以用这种方式来练习。

本文分享自微信公众号 - 数值分析与有限元编程(program_fem)

原文出处及转载信息见文内详细说明,如有侵权,请联系 yunjia_community@tencent.com 删除。

原始发表时间:2019-06-14

本文参与腾讯云自媒体分享计划,欢迎正在阅读的你也加入,一起分享。

我来说两句

0 条评论
登录 后参与评论

相关文章

  • C++STL常用算法

    在iterator对标识元素范围内,查找一对相邻重复元素,找到则返回指向这对元素的第一个元素的迭代器。否则返回past-the-end。

    羊羽shine
  • C++基础字符串

    string类提供了向前和向后遍历的迭代器iterator,迭代器提供了访问各个字符的语法,类似于指针操作,迭代器不检查范围。用string::iterator...

    羊羽shine
  • C++基础构造函数(constructor)

    C++中的类可以定义与类名相同的特殊成员函数,这种与类名相同的成员函数叫做构造函数,构造函数在定义时可以有参数;但是是没有任何返回类型的声明。

    羊羽shine
  • Go基础——Map

    map声明语法为map[K]V,其中K和V分别对应key和value。map中所有的key都有相同的类型,所有的value也有着相同的类型,但是key和valu...

    羊羽shine
  • Golang包——sort

    sort 包 在内部实现了四种基本的排序算法:插入排序(insertionSort)、归并排序(symMerge)、堆排序(heapSort)和快速排序(qui...

    羊羽shine
  • Go基础——常量

    关键字 const 被用于表示常量,常量用于存储不会改变的数据,常量中的数据类型只可以是布尔型、数字型(整数型、浮点型和复数)和字符串型。

    羊羽shine
  • C++基础 模板(Templates)

      泛型即是指具有在多种数据类型上皆可操作的含意。 泛型编程的代表作品:STL 是一种高效、泛型、可交互操作的软件组件。泛型编程最初诞生于 C++中, 目的是为...

    羊羽shine
  • Go基础——reflect反射

    reflect包实现了运行时反射,允许程序操作任意类型的对象。典型用法是用静态类型interface{}保存一个值,通过调用TypeOf获取其动态类型信息,该函...

    羊羽shine
  • C++基础变量

    变量的基本形式是:首先是类型说明符,随后紧跟有一个或者多个的变量名组成的列表,其中含有多个变量名时用逗号隔开,最后以分号结束。

    羊羽shine
  • Java基础——多线程

    进程(process)是计算机中已运行程序的实体。在面向线程设计的系统中,进程本身不是基本运行单位,而是线程的容器。程序本身只是指令、数据及其组织形式的描述,进...

    羊羽shine

扫码关注云+社区

领取腾讯云代金券