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

时间:2013-10-17 15:06:45

标签: c++ linux fortran

好的,所以我想让我的代码工作。这是一个带有“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,我会非常感激!

谢谢,

Josh Derrick

2 个答案:

答案 0 :(得分:0)

你必须意识到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()的{​​{1}} 错误。

答案 1 :(得分:0)

现在我已经第二次看了你的Fortran,我可以告诉你你的语言知识相当有限。

  1. END指定程序的结束,而不是循环/ if语句/模块的结尾(这些是ENDDO / END DOENDIF / {{分别为1}}和END IF
  2. END MODULE表示声明所有变量 必须 ,并且您从未声明IMPLICIT NONE
  3. i不是必需的,因为代码将在完成后返回
  4. 虽然没有必要,但您可能应该将RETURN声明为a并将INTENT(INOUT)声明为b(但是,不可否认,我不确定这是否会导致C- Fortran互操作性)。
  5. 您跳过了INTENT(IN)(看似END,给出了缩进)
  6. 您的Fortran代码应该更像

    IF