matrix相关内容

如何将一个巨大的矩阵逐行写入文件(fortran 90)

我想将包含大量数据的矩阵逐行写入文件.例如,我有一个 100*100 的矩阵,我想在文件中以 100*100 的形式拥有它.但是,它不起作用.以下是我的代码和一些描述.N 和 M 是大约数百个整数.RECL 是我设置文件的预期长度,但这里似乎这个命令不起作用.当 N 设置为 99 时,输出为 198 行M 设置为 200. Vec 是双精度复数矩阵.如何输出 Vec 的值并保持其原始格式 N*M? ..
发布时间:2022-01-14 10:22:28 其他开发

数组内的 Fortran 数组

我正在尝试在 Fortran 中创建一个类似于 MATLAB 中的单元格的数组. 基本上(例如)我正在尝试创建一个数组 X(10) 其中元素 X(1) 是一个维度为 (20,2) 的数组, X(2) 是一个维度为 (25,2) 等的数组 我该怎么做? 解决方案 使用包含单个组件的派生类型来实现特定情况的等价物.元胞数组对应于该派生类型的数组,位于元胞数组每个元素内的数组就是每个 ..
发布时间:2022-01-14 10:12:19 其他开发

如何构造一个包含 9 个较小矩阵的矩阵

我有九个矩阵,其维度为 (N x N)A1(i,j),A2(i,j),A3(i,j),A4(i,j),A5(i,j),A6(i,j),A7(i,j),A8(i,j),A9(i,j) 然后我想构造一个更大的矩阵(3N x 3N),包括这九个矩阵: A = [A1 A2 A3A4 A5 A6[A7 A8 A9] 在fortran中,我可以使用命令行吗 做 i=1,FN做 j=1,FMLA( ..
发布时间:2022-01-14 10:11:33 其他开发

在 Fortran 中使用 ZGETRI 的逆矩阵错误

我正在尝试使用 ZGETRI 计算复杂矩阵的逆矩阵,但是即使它执行没有错误(info = 0),它没有给我正确的逆矩阵,我绝对有不知道错误来自哪里. 程序求解线性隐式无整数 :: i,j,info,lwork整数,维度(3)::ipiv复杂(16),维度(3,3)::C,Cinv,M,LU复杂(16),可分配 :: 工作(:)信息=0lwork=100分配(工作(lwork))ipiv=0工作= ..
发布时间:2022-01-14 10:08:49 其他开发

R,按名称访问矩阵的列向量

在 R 中,我可以通过以下方式访问列矩阵的列向量中的数据: mat2[,1] mat2 的每一列都有一个名称.如何使用 name 属性而不是 [,1] 从第一列检索数据? 例如,假设我的第一列名为“saturn".我想要类似的东西 mat2[,1] == mat2[saturn] 解决方案 以下应该做到: mat2[,'saturn'] 例如: >x ..
发布时间:2022-01-13 18:27:08 其他开发

rbind/data.frame 转换类型层次结构 R

在使用 rbind 和 data.frame 时,谁能系统地向我解释字符/数字/因子之间类型转换的层次结构? 在我的理解中,rbind 放在一个矩阵中,矩阵只能有一种类型.因此,如果存在类型冲突,将转换为什么类型?其他类型的矩阵创建函数(例如 cbind、matrix)是否以同样的方式工作?示例: >sapply(rbind(“a",“b"),类)一个“性格"“性格">sapply(rbi ..
发布时间:2022-01-13 08:17:36 其他开发

LibGDX 移动旋转相机

我在使用 LibGDX 构建的应用程序中遇到了一种情况,我需要相机能够旋转、用户可以移动以及能够缩放.我的用于移动和缩放相机的相机控件可以很好地工作,但是,当相机旋转时,相机会根据该旋转移动.例如,如果相机向左旋转 45 度,并且用户向右拖动以移动相机,它将向右上角移动,因为这是相机右侧现在指向的位置.相反,我需要相机始终相对于屏幕/用户输入移动,无论其旋转如何. 我确实有一个解决这个问题的 ..
发布时间:2022-01-12 16:57:34 其他开发

绘制在 libgdx 中旋转的 BitmapFont

我似乎无法弄清楚如何正确旋转位图字体.我认为您修改了 SpriteBatch 的转换矩阵.但是,尝试旋转会使文本围绕某个点旋转,我不知道如何相对于文本本身旋转它. 解决方案 你可以试试下面的代码: Matrix4 mx4Font = new Matrix4();BitmapFont 字体;SpriteBatch spriteFont;font = new BitmapFont(Gdx.f ..
发布时间:2022-01-12 16:54:20 Java开发

Prolog:寻找所有解决方案

标题可能看起来像一毛钱,但事实并非如此.这个方案的目的就是拿这些类(需求) 需要([[ece2090,1,m,13,16],[ece3520,1,tu,11,14],[ece4420,1,w,13,16]]. 并将他们与具有教授课程证书并且在此期间也有空的大学助教配对(资源;对于 TA 的一天、开始和停止的价值意味着他目前不可用.) resources([[joel, [ece2090,ec ..
发布时间:2022-01-12 10:55:24 其他开发

如何在prolog中转置矩阵

如何将 [[1,2,3][4,5,6][6,7,8]] 之类的列表转置为 [[1,4,6],[2,7,8],[3,6,9]]? 描绘它:我想将矩阵向左翻转 90 度.我该怎么做? 解决方案 不确定你的例子是否正确,但我明白了. 如果使用 SWI-PROLOG,您可以使用 CLPFD 模块,例如所以: :- use_module(library(clpfd)). 允许 ..
发布时间:2022-01-12 10:03:25 其他开发

从 R 中的单个数据帧运行多个线性回归

我有一个包含 21 列的单个国家/地区的出口贸易数据集.第一列表示年份(1962-2014),而其他 20 个是贸易伙伴.我正在尝试对年列和其他列进行线性回归.我已经尝试过这里推荐的方法:运行多个简单的线性回归来自 R 中需要使用的数据框 combn(names(DF), 2, function(x){lm(DF[, x])}, 简化 = FALSE) 然而,这只会产生每对的截距,这对我来说不 ..
发布时间:2022-01-11 09:34:12 其他开发

在numpy中将n个子矩阵编译成NxN矩阵

处理矩阵结构分析中的一个问题.我正在用 Python(使用 Anaconda 3)编写一个程序来分析桁架.每个单独的桁架成员生成一个 4x4 矩阵,总共有 n 个 4x4 矩阵.然后,这些 4x4 矩阵被编译成一个 NxN 矩阵,排列如下,对于矩阵 A、B、C: 如您所见,每个连续的子矩阵都被放置在前一个子矩阵的上一行和下一行.此外,由于桁架的大小和桁架节点(节点)的数量由用户指定,因此 N ..
发布时间:2022-01-10 21:21:16 Python

使用 CUDA 进行动态矩阵乘法

我一直在尝试编写的简单程序的想法是从用户那里获取输入,以查看要乘以多大的矩阵. 我希望将输入 x 乘以 x,目前我不希望将两个不同的大小相乘. 你们会建议我如何完成这项工作? 很抱歉我的问题不够清楚,我想修改这个内核,以便它可以处理任何大小的矩阵(其中 x 和 y 是等价的,以保持简单).而不是 16 的倍数. 我不确定您是否需要我当前的代码,但这里是内核代码: //CU ..
发布时间:2022-01-10 15:58:41 其他开发

CPU 和 GPU 中的 SVD 速度

我在 Matlab R2014a 中测试 svd 似乎没有 CPU vs GPU 加速.我正在使用 GTX 460 卡和 Core 2 duo E8500. 这是我的代码: %test SVDn=10000;%主持人Mh=兰特(n,1000);抽动%[Uh,Sh,Vh]= svd(Mh);svd(兆赫);目录%设备Md = gpuArray.rand(n,1000);抽动%[Ud,Sd,V ..
发布时间:2022-01-10 15:55:26 其他开发

如何使用 numba 在 GPU 上泛化快速矩阵乘法

最近我一直在尝试使用 Numba 库在 Python 中进行 GPU 编程.我一直在使用那里的教程在他们的网站上阅读它,目前我坚持使用他们的示例,可以在这里找到:https://numba.pydata.org/numba-doc/latest/cuda/examples.html.我试图概括一下快速矩阵乘法的示例(其形式为 A*B=C).在测试时,我注意到维度不能完全被每块线程数 (TPB) 整 ..
发布时间:2022-01-10 15:27:01 Python