C++ 类函数调用 fortran 子例程

Question

好的，所以我正在尝试让我的代码工作。它是一个带有类“CArray”的简单 C++ 程序。这个类有 2 个属性，数组大小和值。我希望 C++ 主程序创建类 CArray 的两个实例。在类 CArray 中，我有一个名为“AddArray( CArray )”的函数，它将另一个数组添加到当前数组。我遇到的问题是，我希望函数“AddArray”在 Fortran 中添加两个数组。我知道，要复杂得多，但这就是我需要的。我在将两者链接到类代码中时遇到问题。

#include <iostream>
using namespace std;

class CArray
{
    public:
    CArray();
    ~CArray();
    int Size;
    int*    Val;

    void    SetSize( int );
    void    SetValues();
    void    GetArray();


    extern "C" 
    {
    void    Add( int*, int*, int*, int*);
    void    Subtract( int*, int*, int*, int*);
    void    Muliply( int*, int*, int *, int* );
    }
    void    AddArray( CArray );
    void    SubtractArray( CArray );
    void    MultiplyArray( CArray );
};

这里还有 CArray 函数文件。

#include "Array.h"

#include <iostream>
using namespace std;

CArray::CArray()
{
}


CArray::~CArray()
{
}


void CArray::SetSize( int s )
{
    Size = s;
    for ( int i=0; i<s; i++ )
    {
    Val = new int[Size];
    }
}


void CArray::SetValues()
{
    for ( int i=0; i<Size; i++ )
    {
    cout << "Element " << i+1 << ": ";
    cin >> Val[i];
    }
}


void CArray::GetArray()
{
    for ( int i=0; i<Size; i++ )
    {
    cout << Val[i] << " ";
    }
}


void CArray::AddArray( CArray a )
{
    if ( Size == a.Size )
    {
    Add(&Val, &a.Val);
    }
    else
    {
    cout << "Array dimensions do not agree!" << endl;
    }
}


void CArray::SubtractArray( CArray a )
{
    Subtract( &Val, &a, &Size, &a.Size);
    GetArray();
}

这是我的 Fortran 代码。

module SubtractArrays
    use ico_c_binding
    implicit none

    contains

    subroutine Subtract(a,b,s1,s2) bind(c,name='Subtract')
    integer s1,s2
    integer a(s1),b(s2)

    if ( s1.eq.s2 )
        do i=1,s1
        a(i) = a(i) - b(i)
        end
    return 
    end
end

如果有人可以帮助我设置将整数数组从 C++ 类发送到 Fortran，我将不胜感激!

谢谢，

乔希德里克

Answer 1

您必须意识到 Val 已经是一个指针。所以要调用减法例程，您只需这样做:

void CArray::SubtractArray( CArray a )
{
    Subtract( Val, a.val, &Size, &a.Size);
    GetArray();
}

请注意，我不懂 Fortran，所以我无法对代码的 Fortran 部分发表评论。

---

此外，您的代码中还有很多问题。

1. SetSize(int s) 完全是胡说八道。它分配 s 个数组，每个数组的大小为 s，并泄漏除最后一个以外的所有数组。你可能只是这个意思:

   void CArray::SetSize( int s )
   {
       Size = s;
       delete[] Val;
       Val = new int[Size];
   }
   

2. 你到处都是内存泄漏。您不应该手动管理内存，只需使用 std::vector:

   class CArray
   {
     public:
       CArray();
   
       void    SetSize( size_t newSize );
       void    SetValues();
       void    GetArray();
   
       void    AddArray( const CArray &a );
       void    SubtractArray( const CArray &a );
       void    MultiplyArray( const CArray &a );
   
     private:
       std::vector<int> Val;
       int Size;
   };
   
   extern "C" 
   {
   void    Add( int*, int*, int*, int*);
   void    Subtract( int*, int*, int*, int*);
   void    Muliply( int*, int*, int *, int* );
   }
   

   像这样实现:

   CArray::CArrary : Size(0)
   {}
   
   
   void CArray::SetSize( int newSize )
   {
     Val.resize(newSize);
     Size = newSize;
   }
   
   
   void CArray::SetValues()
   {
       for ( int i=0; i<Size; ++i )
       {
         cout << "Element " << i+1 << ": ";
         cin >> Val[i];
       }
   }
   
   
   void CArray::GetArray()
   {
       for ( int i=0; i<Size; ++)
       {
         cout << Val[i] << " ";
       }
   }
   
   
   void CArray::SubtractArray( const CArray &a )
   {
   
       Subtract( &Val[0], &a.Val[0], &Size, &a.Size);
       GetArray();
   }
   

3. GetArray() 最好重命名为 PrintArray()，理想情况下它应该将流作为参数。 SetValues() 最好命名为 ReadValuesFromInput() 或类似名称。

4. 您应该通过 const 引用而不是值来获取 CArray 类型的参数。我已经在上面的代码中这样做了。

5. 一项职能，一项责任。 SubtractArray() 调用 GetArray() 是错误的。