CODEWITHSUNDEEP
c++fortran
crevell
crevell

Reputation: 460

Fortran subroutines produce "Undefined symbols" error when linking to c++ program

I am writing a short program to test calling the fortran Stripack library from c++. The c++ and fortran files each compile successfully, but errors occur when linking.

The c++ code is as follows:

#include <iostream>
#include <cmath>

#ifdef __cplusplus
extern"C" {
    #endif
    void trmesh_(int&,float[],float[],float[],int[],int[],int[],int&,int[],int[],float[],int&);
    void trlist2_(int&,int[],int[],int[],int&,int[][3],int&);
    #ifdef __cplusplus
}
#endif


int main(){

// Variables for distributing points on sphere.
  int polar = 16;
  int azimuth = 32;
  int n = polar*azimuth-azimuth+2;
  float radius=1.0;

// Define variables needed by Stripack
  float xs[n];
  float ys[n];
  float zs[n];
  int list[6*(n-2)];
  int lptr[6*(n-2)];
  int lend[6*(n-2)];
  int near[n];
  int next[n];
  float dist[n];
  int ltri[2*n-4][3];
  int lnew;
  int ier;
  int nt;

// Distribute n points on surface of unit sphere .
// xs, ys, zs store x, y, and z components pf each point position.
  zs[0] = 1;
  xs[0] = 0;
  ys[0] = 0;
  zs[n] = -1;
  xs[n] = 0;
  ys[n] = 0;
  for (int ii=1; ii<polar; ii++){
    for (int jj=0; jj<azimuth; jj++){
      zs[(ii-1)*azimuth+jj+1] = radius*cos(ii*M_PI/polar);
      xs[(ii-1)*azimuth+jj+1] = radius*sin(ii*M_PI/polar)*sin(jj*2*M_PI/azimuth);
      ys[(ii-1)*azimuth+jj+1] = radius*sin(ii*M_PI/polar)*cos(jj*2*M_PI/azimuth);
    }
  }

// Call stripack subroutines to obtain list of triangles ltri
  trmesh_(n,xs,ys,zs,list,lptr,lend,lnew,near,next,dist,ier);
  trlist2_(n,list,lptr,lend,nt,ltri,ier);

// Output list of triangles
  for (int ii =0; ii<n; ii++){
    std::cout << ltri[ii][0] << " " << ltri[ii][1] << " " << ltri[ii][2] << std::endl;
  }

}

I compile the files as follows:

ifort -c stripack.f90
clang++ -c -O0 -std=c++11 -c -o main.o main.cpp -g
clang++ -o main stripack.o main.o

The first two compilations work fine, but the last one produces the following results. It seems like the subroutines in the fortran file can't find standard fortran functions? I have tried with gfortran and the same problem occurs. Any suggestions as to what is going on would be greatly appreciated. 

Undefined symbols for architecture x86_64:
  "___libm_sse2_sincos", referenced from:
      _trplot_ in stripack.o
      _vrplot_ in stripack.o
  "___svml_sincos2", referenced from:
      _trans_ in stripack.o
  "_for_date_and_time", referenced from:
      _timestamp_ in stripack.o
  "_for_stop_core", referenced from:
      _trmesh_ in stripack.o
      _addnod_ in stripack.o
  "_for_trim", referenced from:
      _timestamp_ in stripack.o
  "_for_write_seq_fmt", referenced from:
      _delnod_ in stripack.o
      _edge_ in stripack.o
      _timestamp_ in stripack.o
      _trlprt_ in stripack.o
      _trmesh_ in stripack.o
      _addnod_ in stripack.o
      _trplot_ in stripack.o
      ...
  "_for_write_seq_fmt_xmit", referenced from:
      _delnod_ in stripack.o
      _edge_ in stripack.o
      _timestamp_ in stripack.o
      _trlprt_ in stripack.o
      _trmesh_ in stripack.o
      _addnod_ in stripack.o
      _trplot_ in stripack.o
      ...
ld: symbol(s) not found for architecture x86_64
clang: error: linker command failed with exit code 1 (use -v to see invocation)

Upvotes: 0

Views: 964

Answers (1)

KL-Yang
KL-Yang

Reputation: 401

I will demonstrate this is a link issue by example, you need do a little bit more research to solve the problem, as the information you provide is not complete. 

!fortran code, named as x.f90
subroutine testFDLL(str, n) bind(c, name='testFDLL_as_C')
    use ISO_C_BINDING
    integer(c_int), value :: n
    character(kind=c_char), intent(in) :: str(n)
    write(6,*)" Hello FORTRAN : let us do something ...",str
    return
end

The following C code is used for demonstration (you have got C++ mostly right already).

//c named as y.c
#include <stdio.h>
#include <string.h>

int main()
{
    void testFDLL_as_C(char *str, int n);
    char str[] = "Hello from C";
    testFDLL_as_C(str, strlen(str));
    return 0;
}

If you compile and link use the following

ifort -c x.f90
gcc y.c x.o -W -Wall

Depend on version of ifort and OS, should get error similar as the following

x.o: In function `testFDLL_as_C':
x.f90:(.text+0x42): undefined reference to `for_write_seq_lis'
x.f90:(.text+0x74): undefined reference to `for_write_seq_lis_xmit'
collect2: error: ld returned 1 exit status

You may noticed the undefined reference name pattern is similar with yours, if you link with 

gcc y.c x.o -W -Wall -L/path/to/your/ifort_lib -lifcore -ldl

The problem should be solved. Depend on the FORTRAN feature you used, you may need link some more ifort library. This part need you do some research and figure out.

Upvotes: 1

Related Questions