如何最好地实现与Matlab相交功能等效的Fortran

[poster]

令我惊讶的是，这里没有内在的例程，但是我需要一个代码片段，该片段应采用两个预先排序的1D整数数组，并以有效的方式返回两个数组中都存在的元素。

integer, allocatable :: A (:), B(:), C(:) A=[1,2,3,6,10] B=[3,10,5] C=intersect(A,B) !< returns [3,10] 在Fortran中实现此功能的最佳方法是什么，或者我可以看一个好的开源实现？我也无法在数字食谱或Rosetta代码中找到它。



更多&回答...