如何从Fortran90矩阵中删除随机行

huangapple
go评论85阅读模式
英文:

How to remove random rows from a matrix in Fortran90

问题

好的，以下是翻译好的内容：

问题如下。在我的作业中，我应该创建一个(3,1000)的矩阵，然后创建它后，我应该随机删除271行矩阵。最后，我应该获得一个(3,729)的矩阵，其中包含所有初始值，除了我删除的与随机数对应的那个。

为了做到这一点，我编写了这个脚本，其中
Mc=1000
N=729

  1. h=1
  2. do a=1, (Mc-N)
  4.    call random_number(rand_num)
  5.    rand_num=int(rand_num*Mc)-1
  6.    !print*, rand_num
  8.    do b=1, Mc
  10.       if (b /= rand_num) then
  12.       NMatrix(:,h) = Matrix(:,b)
  14.       h=h+1
  16.       else
  18.       h=h
  20.       endif
  21.     enddo
  22. enddo

但是当我运行它时，它报告了这个错误：程序接收到信号SIGABRT：进程中止信号。
我应该怎么办？请帮我，这非常重要。

我知道do循环和数组的内存存在一些问题，但我无法弄清楚具体问题所在。你能向我提问吗？

英文:

OK so here is the problem. In my assignment I should create a (3,1000) matrix and after I created it I should delete RANDOMLY 271 rows of the matrix. At the end I should obtain a (3,729) Matrix with all the initial values except the one I deleted corresponding to the random numbers.

To do so I wrote this script where
Mc=1000
N=729

  1. h=1
  2. do a=1, (Mc-N)
  4.    call random_number(rand_num)
  5.    rand_num=int(rand_num*Mc)-1
  6.    !print*, rand_num
  8.    do b=1, Mc
  10.       if (b /= rand_num) then
  12.       NMatrix(:,h) = Matrix(:,b)
  14.       h=h+1
  16.       else
  18.       h=h
  20.       endif
  21.     enddo
  22. enddo

But when I run it it reports me this error: Program receiver signal SIGABRT : Process abort signal.
What should I do? Please it's really important

I know there is some problem with the do loops and the memory of the arrays, but I cannot get to the point. Can you please ask me?

答案1

得分: 0

这是您要求的翻译部分：

  1. 这是您要求的翻译部分：
  2. ```fortran
  3. 这是您要求的翻译部分：
  4. ```fortran
  5. 这是您要求的翻译部分：
  6. ```fortran
  7. 这是您要求的翻译部分：
  8. ```fortran
  9. 这是您要求的翻译部分：
  10. ```fortran
  11. 这是您要求的翻译部分：
  12. ```fortran
  13. 这是您要求的翻译部分：
  14. ```fortran
  15. 这是您要求的翻译部分：
  17. ```fortran
  18. 这是您要求的翻译部分：
  20. ```fortran
  21. 这是您要求的翻译部分：
  22. ```fortran
  23. 这是您要求的翻译部分：
  24. ```fortran
  25. 这是您要求的翻译部分：
  27. ```fortran
  28. 这是您要求的翻译部分：
  30. ```fortran
  31. 这是您要求的翻译部分：
  33. ```fortran
  34. 这是您要求的翻译部分：
  36. ```fortran
  37. 这是您要求的翻译部分：
  39. ```fortran
  40. 这是您要求的翻译部分：
  42. ```fortran
  43. 这是您要求的翻译部分：
  45. ```fortran
  46. 这是您要求的翻译部分：
  48. ```fortran
  49. 这是您要求的翻译部分：
  51. ```fortran
  52. 这是您要求的翻译部分：
  54. ```fortran
  55. 这是您要求的翻译部分：
  57. ```fortran
  58. 这是您要求的翻译部分：
  60. ```fortran
  61. 这是您要求的翻译部分：
  63. ```fortran
  64. 这是您要求的翻译部分：
  66. ```fortran
  67. 这是您要求的翻译部分：
  69. ```fortran
  70. 这是您要求的翻译部分：
  72. ```fortran
  73. 这是您要求的翻译部分：
  75. ```fortran
  76. 这是您要求的翻译部分：
  78. ```fortran
  79. 这是您要求的翻译部分：
  81. ```fortran
  82. 这是您要求的翻译部分：
  84. ```fortran
  85. 这是您要求的翻译部分：
  87. ```fortran
  88. 这是您要求的翻译部分：
  90. ```fortran
  91. 这是您要求的翻译部分：
  93. ```fortran
  94. 这是您要求的翻译部分：
  96. ```fortran
  97. 这是您要求的翻译部分：
  99. ```fortran
  100. 这是您要求的翻译部分：
  102. ```fortran
  103. 这是您要求的翻译部分：
  105. ```fortran
  106. 这是您要求的翻译部分：
  108. ```fortran
  109. 这是您要求的翻译部分：
  111. ```fortran
  112. 这是您要求的翻译部分：
  114. ```fortran
  115. 这是您要求的翻译部分：
  117. ```fortran
  118. 这是您要求的翻译部分：
  120. ```fortran
  121. 这是您要求的翻译部分：
  123. ```fortran
  124. 这是您要求的翻译部分：
  126. ```fortran
  127. 这是您要求的翻译部分：
  129. ```fortran
  130. 这是您要求的翻译部分：
  132. ```fortran
  133. 这是您要求的翻译部分：
  135. ```fortran
  136. 这是您要求的翻译部分：
  138. ```fortran
  139. 这是您要求的翻译部分：
  141. ```fortran
  142. 这是您要求的翻译部分：
  144. ```fortran
  145. 这是您要求的翻译部分：
  147. ```fortran
  148. 这是您要求的翻译部分：
  150. ```fortran
  151. 这是您要求的翻译部分：
  153. ```fortran
  154. 这是您要求的翻译部分：
  156. ```fortran
  157. 这是您要求的翻译部分：
  159. ```fortran
  160. 这是您要求的翻译部分：
  162. ```fortran
  163. 这是您要求的翻译部分：
  165. ```fortran
  166. 这是您要求的翻译部分：
  168. ```fortran
  169. 这是您要求的翻译部分：
  171. ```fortran
  172. 这是您要求的翻译部分：
  174. ```fortran
  175. 这是您要求的翻译部分：
  177. ```fortran
  178. 这是您要求的翻译部分：
  180. ```fortran
  181. 这是您要求的翻译部分：
  183. ```fortran
  184. 这是您要求的翻译部分：
  186. ```fortran
  187. 这是您要求的翻译部分：
  189. ```fortran
  190. 这是您要求的翻译部分：
  192. ```fortran
  193. 这是您要求的翻译部分：
  195. ```fortran
  196. 这是您要求的翻译部分：
  198. ```fortran
  199. 这是您要求的翻译部分：
  201. ```fortran
  202. 这是您要求的翻译部分：
  204. ```fortran
  205. 这是您要求的翻译部分：
  207. ```fortran
  208. 这是您要求的翻译部分：
  210. ```fortran
  211. 这是您要求的翻译部分：
  213. ```fortran
  214. 这是您要求的翻译部分：
  216. ```fortran
  217. 这是您要求的翻译部分：
  219. ```fortran
  220. 这是您要求的翻译部分：
  222. ```fortran
  223. 这是您要求的翻译部分：
  225. ```fortran
  226. 这是您要求的翻译部分：
  228. ```fortran
  229. 这是您要求的翻译部分：
  231. ```fortran
  232. 这是您要求的翻译部分：
  234. ```fortran
  235. 这是您要求的翻译部分：
  237. ```fortran
  238. 这是您要求的翻译部分：
  240. ```fortran
  241. 这是您要求的翻译部分：
  243. ```fortran
  244. 这是您要求的
  246. <details>
  247. <summary>英文:</summary>
  249. Here an example ouput of what you are asking for (I think).
  251. ```lang-text
  252. rows=3, cols=1000
  253. 1.000 2.000 3.000 4.000 5.000 6.000 7.000 8.000 9.000 10.00 11.00 12.00 ...
  254. 1001. 1002. 1003. 1004. 1005. 1006. 1007. 1008. 1009. 1010. 1011. 1012. ...
  255. 2001. 2002. 2003. 2004. 2005. 2006. 2007. 2008. 2009. 2010. 2011. 2012. ...
  257. rows=3, cols=729
  258. 1.000 3.000 4.000 9.000 10.00 11.00 12.00 13.00 14.00 16.00 17.00 19.00 ...
  259. 1001. 1003. 1004. 1009. 1010. 1011. 1012. 1013. 1014. 1016. 1017. 1019. ...
  260. 2001. 2003. 2004. 2009. 2010. 2011. 2012. 2013. 2014. 2016. 2017. 2019. ...

Note that by convention, the 1st index is the rows, and the 2nd index is the columns. So you are asking to remove 271 columns from the matrix.

I decided to use an array of index values indx=[1,2,3,4 .. n], which I shuffle randomly, then slice off the last 271 values (keep n entries) and then sort back to increasing order.

Finally, the result is essentially NMatrix(:,:) = Matrix(:, indx) or in my case I call them a and b in my function.

  1. function RemoveRandomColumns(a,n_col) result(b)
  2. real, intent(in) :: a(:,:)
  3. integer, intent(in) :: n_col
  4. real, allocatable :: b(:,:)
  5. integer :: i,j,n,m,t
  6. integer, allocatable :: indx(:)
  7. real :: u !used in shuffle part
  8. m = size(a, 1)
  9. n = size(a, 2)
  10. ! must remove m-n_col columns from a
  11. allocate(indx(n))
  12. indx = [ (i,i=1,n) ]
  13. ! indx = [1,2,3,4,5,6,7,8,..,n]
  14. ! suffle the index
  15. call RANDOM_SEED()
  16. do i=1, n
  17. call RANDOM_NUMBER(u)
  18. j = 1 + floor(u*(n-1))
  19. t = indx(i)
  20. indx(i) = indx(j)
  21. indx(j) = t
  22. end do
  23. ! indx = [18,10,3,16,21,...]
  24. ! take only n_col items
  25. indx = indx(1:n_col)
  26. ! indx = [18,10,3,...]
  27. ! sort the index
  28. do i=1, n_col
  29. do j=1, i-1
  30. if(indx(i)&lt;indx(j)) then
  31. t = indx(i)
  32. indx(i) = indx(j)
  33. indx(j) = t
  34. end if
  35. end do
  36. end do
  37. ! indx = [2,3,4,6,7,9,..,n]
  38. ! now slice the array to get the resukt
  39. allocate(b(m,n_col))
  40. do j=1, m
  41. b(j,:) = a(j,indx)
  42. end do
  43. end function

and the sample code to test the above is

  1. program FortranProgram1
  2. use, intrinsic :: iso_fortran_env
  3. implicit none
  4. integer, parameter :: Mc = 1000, n = 729
  5. ! Variables
  6. real :: matrix(3, Mc)
  7. real :: reduced(3, n)
  8. integer :: i,j
  9. do j=1, 3
  10. do i=1, Mc
  11. matrix(j,i) = real(Mc*(j-1)+i)
  12. end do
  13. end do
  14. print &#39;(&quot;rows=&quot;,g0,&quot;, cols=&quot;,g0)&#39;, size(matrix,1), size(matrix,2)
  15. print &#39;(*(g0.4,&quot; &quot;))&#39;, matrix(1,1:12), &quot;...&quot;
  16. print &#39;(*(g0.4,&quot; &quot;))&#39;, matrix(2,1:12), &quot;...&quot;
  17. print &#39;(*(g0.4,&quot; &quot;))&#39;, matrix(3,1:12), &quot;...&quot;
  18. print *, &quot;&quot;
  19. reduced = RemoveRandomColumns(matrix, n)
  20. print &#39;(&quot;rows=&quot;,g0,&quot;, cols=&quot;,g0)&#39;, size(reduced,1), size(reduced,2)
  21. print &#39;(*(g0.4,&quot; &quot;))&#39;, reduced(1,1:12), &quot;...&quot;
  22. print &#39;(*(g0.4,&quot; &quot;))&#39;, reduced(2,1:12), &quot;...&quot;
  23. print &#39;(*(g0.4,&quot; &quot;))&#39;, reduced(3,1:12), &quot;...&quot;
  24. contains
  25. function RemoveRandomColumns(a,n_col) result(b)
  26. ...
  27. end function
  28. end program

huangapple
  • 本文由 发表于 2023年2月16日 18:22:59
  • 转载请务必保留本文链接：https://go.coder-hub.com/75470847.html
  • arrays
  • fortran
  • memory
  • sigabrt
匿名

发表评论

匿名网友

:?: :razz: :sad: :evil: :!: :smile: :oops: :grin: :eek: :shock: :???: :cool: :lol: :mad: :twisted: :roll: :wink: :idea: :arrow: :neutral: :cry: :mrgreen:

确定