用 Fortran 编写矩阵

2024-01-11

我需要用 FORTRAN 输出一个矩阵。我有一个计算值的工作代码,但我得到的不是矩阵,而是一列。矩阵很大,ixj = ~2000x2000。

这是我的示例代码:

  open(19, file="results1.txt", status="old", position="rewind", 
 & action="write") 

  do j=0,p
  do i=0,o
  write(19,*) mat_user_yield_surface(d, eps(i), deps(j), 200.0d0)
  end do
  end do

  close(19)

使用隐含的 do 循环:

do j=0,p
   write(19,'(2000g22.14)') (mat_user_yield_surface(d, eps(i), deps(j),200.0d0),i=0,o)
end do

我建议不要使用“o”作为变量名,因为它很容易与零混淆。

本文内容由网友自发贡献,版权归原作者所有,本站不承担相应法律责任。如您发现有涉嫌抄袭侵权的内容,请联系:hwhale#tublm.com(使用前将#替换为@)

用 Fortran 编写矩阵 的相关文章

  • 如何使用 Julia 查找矩阵中的连通分量

    假设我有以下矩阵 此处用 Julia 语言定义 mat 1 1 0 0 0 1 1 0 0 0 0 0 0 0 1 0 0 0 1 1 将一组值为 1 的相邻元素视为一个 分量 如何识别该矩阵有 2 个分量以及每个分量由哪些顶点组成 对于矩
  • 如何找到在matlab中重复的矩阵的每一行的索引?

    我想找到矩阵中所有有重复项的行的索引 例如 A 1 2 3 4 1 2 3 4 2 3 4 5 1 2 3 4 6 5 4 3 要返回的向量将是 1 2 4 很多类似的问题建议使用unique函数 我已经尝试过 但我能得到的最接近我想要的功
  • 使用矩阵代数来操作字符串:可行吗?

    我正在尝试使用矩阵代数来操作字符串 这意味着能够使用字符串或字符串数 组的串联和粘贴来实现多个类似矩阵的结构 我之前尝试在 R 上实现这个东西 但这是不可能的 因为矩阵只能有一维条目 我希望足够的与语言无关和抽象 但为了清楚起见 我将使用类
  • R矩阵右除

    在 R 中进行左右矩阵除法的最简洁 最快 数值最稳定 最 R 惯用的方法是什么 我理解左除法inv A B通常是用solve a b 但是怎么样B inv A 确实是最好的计算方法t solve t A t B It is B solve
  • Matlab没有优化以下内容吗?

    我有一个很长的向量 1xrv 和一个很长的向量w1xs 和一个矩阵Arxs 它是稀疏的 但维度非常大 我期望 Matlab 对以下内容进行优化 这样我就不会遇到内存问题 A v w 但看起来 Matlab 实际上是在尝试生成完整的v w矩阵
  • R 中的自定义对比:对比系数矩阵或对比矩阵/编码方案?以及如何到达那里?

    自定义对比在分析中应用非常广泛 例如 这个三水平因子的第 1 级和第 3 级的 DV 值是否有显着差异 直观上 这种对比用单元均值表示为 c 1 0 1 这些对比中的一个或多个 以列的形式绑定 形成对比系数矩阵 例如 mat matrix
  • 更快的四元数向量乘法不起作用

    我的数学库需要一个更快的四元数向量乘法例程 现在我正在使用规范v qv q 1 它产生的结果与向量乘以由四元数组成的矩阵相同 所以我对它的正确性充满信心 到目前为止 我已经实现了 3 种替代 更快 的方法 1 我不知道我从哪里得到这个 v
  • R::bigmemory - 如何创建角色big.matrix?

    我尝试使用bigmemory封装在R我一开始就陷入困境 我愿意 temp lt matrix paste a 1 10 5 2 并得到一个字符矩阵 没关系 但后来我尝试 x lt as big matrix temp type char 我
  • 如何使用 MPI_Scatterv 将矩阵的行发送到所有进程?

    我正在使用 MPI 接口 我想分割一个矩阵 按行 并将各个部分分配给每个进程 例如 我有这个7x7的方阵M M 0 00 1 00 2 00 3 00 4 00 5 00 6 00 7 00 8 00 9 00 10 00 11 00 12
  • 如何将 numpy.matrix 提高到非整数幂?

    The 运算符为numpy matrix不支持非整数幂 gt gt gt m matrix 1 0 0 5 0 5 gt gt gt m 2 5 TypeError exponent must be an integer 我想要的是 oct
  • 在重载算术运算符中调用 C++ 析构函数

    我有一个用于神经网络程序和重载算术运算符的定制矩阵库 这是类声明 class Matrix public int m int n double mat Matrix int int Matrix int Matrix const Matri
  • 单个对象的 Monogame XNA 变换矩阵?

    我读过一些解释 XNA Monogame 变换矩阵的教程 问题是这些矩阵应用于 SpriteBatch Begin matrix 这意味着所有 Draw 代码都将被转换 如何将变换矩阵应用于单个可绘制对象 就我而言 我想转换滚动背景 使其自
  • 在 Visual Studio 2010 中从 Fortran 调用 C++ 函数

    我想从 Fortran 调用 C 函数 为此 我在 Visual Studio 2010 中创建了一个 FORTRAN 项目 之后 我将一个 Cpp 项目添加到该 FORTRAN 项目中 当我要构建程序时出现以下错误 Error 1 unr
  • Android获取地球坐标系上的加速度计

    我想在 Android 上获取加速度计并将它们放在地球坐标系上 就像这个主题一样从设备坐标系到绝对坐标系的加速度 https stackoverflow com questions 11578636 acceleration from de
  • 将矩阵 csv 数据导入 R - 如何转换为数据框

    我有一组以矩阵格式保存的 csv 数据附图是矩阵的示例 https i stack imgur com kTIN6 png我想将数据加载到 R 中并将其存储为包含 x Year x Death x ASMR 的数据框 我怎样才能做到这一点
  • 如何使用 BLAS 执行向量矩阵乘法?

    BLAS 定义了 GEMV 矩阵向量乘法 2 级运算 如何使用 BLAS 库执行向量矩阵乘法 这可能很明显 但我不知道如何使用 BLAS 运算进行此乘法 我原本期望进行 GEVM 手术 M x N 矩阵与 N x 1 向量的矩阵向量乘法将得
  • 向具有 TARGET 属性的虚拟参数的过程提供不具有 TARGET 属性的参数

    在 Fortran 语言中 向具有 TARGET 属性的虚拟参数的过程提供不具有 TARGET 属性的参数应该会导致无效代码 但是 当使用 gfortran 5 1 0 或 ifort 14 0 0 编译以下代码时 不会检测到错误 并且程序
  • Julia:按第 2 列然后按第 3 列对矩阵进行排序

    我想对我的矩阵进行排序A按第 2 列 然后按第 3 列 A round randn 100 4 也许是这样的 sort A 0 2 3 100x4 Array Float64 2 0 0 2 0 2 0 1 0 1 0 2 0 1 0 1
  • 如何在Python中实现部分旋转的LU分解?

    我想实现我自己的 LU 分解 P L U my lu A 以便给定矩阵 A 通过部分旋转计算 LU 分解 但我只知道如何在不旋转的情况下做到这一点 有人可以帮忙做部分旋转吗 def lu A import numpy as np Retur
  • & 位于第 5 列和行尾

    我今天遇到一些代码 如下所示 subroutine sub hello world this routine takes a crazy large number of arguments so many that it gets spli

随机推荐

  • 使用 Dagger 模块而不使用“injects”指令

    我试图让 Dagger 在 Module 注释中没有 injects 指令的情况下工作 我的测试项目基于Android 简单 Dagger 示例 https github com square dagger tree master exam
  • 如何迭代结构体的元素?

    我正在编写一个用于加密网络通信的小型客户端 服务器程序 并具有以下结构以允许端点协商功能 struct KeyExchangePacket kexinit u8 replay cookie u8 32 kex algorithms Stri
  • Java - 日期构造函数接受日期字符串,但已弃用。尝试过替代方案,但没有运气

    String temp date 07 28 2011 11 06 37 AM Date date new Date temp date Depricated SimpleDateFormat sdf new SimpleDateForma
  • 在自定义 UIControl 对象中定义自定义触摸区域

    我正在创建一个自定义 UIControl 对象 详细信息here http www thinkandbuild it how to build a custom control in ios 除了触摸区域外 一切都运行良好 我想找到一种方法
  • 如何检查 CGPoint 是否在 UIImageView 内部?

    In touchesBegan CGPoint touch point touches anyObject locationInView self view 有数十个UIImageView周围 存储在NSMutableArray image
  • 如何在打印前隐藏按钮并在打印过程完成后显示它?

    我有一个打印按钮id print req 我编写了一些用于打印页面的Javascript代码 这是通过单击此按钮触发的 我还想在打印之前隐藏此按钮 并在整个打印过程后显示它 我的意思是不在我的打印文档中打印按钮 这是我的代码 documen
  • 如何在 Doctrine 映射中描述多列外键

    我们有一个数据库模式 以简化的 稍微做作的 形式 看起来像 从用户到域的外键设置在列 domainId groupId 上 以保证引用完整性 这种结构可以很好地达到预期目的 但是 对于与同一数据库通信的新应用程序 我现在需要为映射上述结构的
  • 将不同的参数传递给每个映射器

    我有一份使用多个映射器和一个减速器的工作 映射器几乎相同 除了 a 的值不同之外String他们用来产生结果 目前我有几个类 每个类对应一个值String我提到过 感觉应该有更好的方法 不需要那么多的代码重复 有没有办法通过这些String
  • 如何通过matplotlib删除plot_wireframe中的隐藏线?

    是否matplotlib在 python 2 x 中支持隐藏线删除吗 我自己该如何实施呢 这是此处解决的另一篇文章的进一步问题 如何通过Python获取3D彩色表面 https stackoverflow com questions 249
  • React Hook 表单错误 - 类型“UseFormRegister”不可分配给类型“LegacyRef”|不明确的'

    我正在阅读一本 React 教程书 有一次 作者向我展示了如何构建一个简单的表单 效果很好 然后 他又将其升级为 React Hook Form 一切都按照说明运行良好 直到最后一行 他说了以下内容 将 ref 属性添加到 Header 组
  • UIView 越界时隐藏子视图

    我有一个视图 比如说 100x100 它有一组 uiviews 作为子视图 假设 30x30 如果子视图具有左上角坐标 90 90 我希望只能看到该子视图的一部分 但我仍然可以看到整个区域 即使它超出了父母的范围 问题是如何使 uiview
  • 显示 Silverlight 的组合框组标题

    我想展示一个ComboBox with OPTGROUPSilverlight 中的样式标头分组 我找到的每个网站 包括有关 SO 的问题 都解决了此链接到过时链接的问题 并且轻松地没有显示可供我使用的代码片段 E g 那么我该怎么做呢 看
  • InstallReferrerClient.startConnection() 抛出错误

    我的一些用户在使用 InstallReferrerClient startConnection 时收到以下错误 我正在跟进android 开发者教程 https android developers googleblog com 2017
  • 如何使用实体框架仅更新一个字段?

    这是表格 Users UserId UserName Password EmailAddress 和代码 public void ChangePassword int userId string password code to updat
  • SQL Server 数据库字段处理韩文和中文字符

    SQL Server中是否有可能有一个字段可以存储中文 韩文和欧洲字符 我的汉字就变成了 数据类型是NVARCHAR以及 NVARCHAR是正确的类型 它将所有内容存储在 2 字节 Unicode 中 与您合作时需要注意的是NVARCHAR
  • Oracle SQL 语句动态模式变量

    我有数百行代码来创建表 但模式名称是硬编码的 如所有 DDL 语句中的 SCHEMA NAME1 TABLE NAME 我怎样才能将其作为变量并在所有地方使用 以便我们可以轻松地在一个地方更改 SCHEMA NAME 要求的多个位置 请给出
  • jquery 模式对话框 onclick?

    我要这个 http jqueryui com demos dialog modal message http jqueryui com demos dialog modal message 当您单击 ClickMe 时发生 这个怎么做 di
  • Firebase 数据库和应用本地化

    我正在为 iOS 和 Android 平台构建一个问答游戏 我希望能够处理本地化 我正在使用 Firebase 的实时数据库解决方案来解决游戏中存在的所有问题 我将问题硬编码到 Firebase 的数据库中 每个问题对象都有 2 个参数 I
  • 基于多个切片器值 Power BI 显示用户计数

    我正在研究一个Viewership table它告诉客户手表资产的数量 有两个资产过滤器 一种用于观看 另一种用于未观看 根据资产过滤器 我需要显示观看节目的客户数量和未观看节目的客户数量 下面是我用过的方法 首先 创建一个新表 Asset
  • 用 Fortran 编写矩阵

    我需要用 FORTRAN 输出一个矩阵 我有一个计算值的工作代码 但我得到的不是矩阵 而是一列 矩阵很大 ixj 2000x2000 这是我的示例代码 open 19 file results1 txt status old positio