FortranCPP.hpp File Reference

Go to the source code of this file.

Macros
#define	BEGIN_EXTERN_FORTRAN
#define	END_EXTERN_FORTRAN
#define	EXTERN_FFUNCTION(x)
#define	CALL_FFUNCTION(x)
#define	EXTERN_FSUBROUTINE(x)
#define	CALL_FSUBROUTINE(x)
#define	EXTERN_FMODFUNCTION(x, y)
#define	CALL_FMODFUNCTION(x, y)
#define	EXTERN_FMODSUBROUTINE(x, y)
#define	CALL_FMODSUBROUTINE(x, y)
#define	FREAL
#define	FREAL8
#define	FREAL16
#define	FINT
#define	FINT2
#define	FINT8
#define	FLOGICAL
#define	FLOGICAL1
#define	FCHARACTER
#define	FSTRING
#define	FSTRINGLEN
#define	STD2FSTRING(str)
#define	STD2FSTRINGLEN(str)
#define	STD2FSTRINGFULL(str)

Detailed Description

Author

Jean-Christophe Fabre jean-.nosp@m.chri.nosp@m.stoph.nosp@m.e.fa.nosp@m.bre@s.nosp@m.upag.nosp@m.ro.in.nosp@m.ra.f.nosp@m.r

See Also

http://www-h.eng.cam.ac.uk/help/tpl/languages/mixinglanguages.html

http://www.neurophys.wisc.edu/comp/docs/notes/not017.html

http://www.yolinux.com/TUTORIALS/LinuxTutorialMixingFortranAndC.html

http://arnholm.org/software/cppf77/cppf77.htm

http://www.aei.mpg.de/~jthorn/c2f.html

Below are examples of C++ simulators calling Fortran subroutines.

fortran source code:

subroutine multrealvalue(value,mult,result)

implicit none
real*8 value
real*8 mult
real*8 result

result=value*mult

return

end


! ==============================================================================


subroutine multintvalue(value,mult,result)

implicit none
integer value
integer mult
integer result

result=value*mult

return
end


! ==============================================================================


subroutine multrealmatrix(matrix,dim1,dim2,mult,result)

implicit none
integer dim1
integer dim2
real*8 matrix(dim2,dim1)
integer mult
real*8 result(dim2,dim1)
integer i, j, k
subroutine multrealvalue(value,mult,result)

implicit none
real*8 value
real*8 mult
real*8 result

result=value*mult

return

end


! ==============================================================================


subroutine multintvalue(value,mult,result)

implicit none
integer value
integer mult
integer result

result=value*mult

return
end


! ==============================================================================


subroutine multrealmatrix(matrix,dim1,dim2,mult,result)

implicit none
integer dim1
integer dim2
real*8 matrix(dim2,dim1)
integer mult
real*8 result(dim2,dim1)
integer i, j, k

do j=1,dim1
  do i=1,dim2
    result(i,j) = matrix(i,j) * mult
  end do
end do


return
end

do j=1,dim1
  do i=1,dim2
    result(i,j) = matrix(i,j) * mult
  end do
end do


return
end

C++ code calling fortran subroutines:

BEGIN_EXTERN_FORTRAN
  EXTERN_FSUBROUTINE(multrealvalue)(FREAL8*,FREAL8*,FREAL8*);
  EXTERN_FSUBROUTINE(multintvalue)(FINT*,FINT*,FINT*);
  EXTERN_FSUBROUTINE(multrealmatrix)(FREAL8*,FINT*,FINT*,FINT*,FREAL8*);
END_EXTERN_FORTRAN

// =====================================================================
// =====================================================================

bool FortranFunction::runStep(const openfluid::base::SimulationStatus* SimStatus)
{
  int i;


  // ====== double ======

  double DValue, DMult, DResult;
  DValue = 1.5436;
  DMult = 2.5;
  DResult = 0.0;

  CALL_FSUBROUTINE(multrealvalue)(&DValue,&DMult,&DResult);

  if (abs(DResult - (DValue*DMult)) > m_Precision)
    OPENFLUID_RaiseError("tests.fortran","incorrect fortran call (multrealvalue)");


  // ====== int ======

  int IValue, IMult, IResult;
  IValue = 45;
  IMult = 18;
  IResult = 0;

  CALL_FSUBROUTINE(multintvalue)(&IValue,&IMult,&IResult);

  if (IResult != (IValue*IMult))
    OPENFLUID_RaiseError("tests.fortran","incorrect fortran call (multintvalue)");


  // ====== matrix of double ======

  int MMult, MDim1,MDim2;
  double *MValue;
  double *MResult;
  MMult = 3;
  MDim1 = 2;
  MDim2 = 3;


  MValue = new double[MDim1*MDim2];
  MResult = new double[MDim1*MDim2];

  for (i=0; i < MDim1*MDim2;i++) MValue[i] = 1.5;
  for (i=0; i < MDim1*MDim2;i++) MResult[i] = 0.0;

  CALL_FSUBROUTINE(multrealmatrix)(MValue,&MDim1,&MDim2,&MMult,MResult);

  for (i=0; i < MDim1*MDim2;i++)
  {
    if (abs(MResult[i] - (MValue[i] * MMult)) > m_Precision)
      OPENFLUID_RaiseError("tests.fortran","incorrect fortran call (multrealmatrix)");
  }

  return true;
}

Macro Definition Documentation

#define BEGIN_EXTERN_FORTRAN

#define CALL_FFUNCTION

(

x

)

Macro for calling an external fortran function

Parameters

[in]

x

the name of the function

#define CALL_FMODFUNCTION	(		x,
			y
	)

Macro for calling an external fortran90 function in a module

Parameters

[in]	x	the name of the module
[in]	y	the name of the function

#define CALL_FMODSUBROUTINE	(		x,
			y
	)

Macro for calling an external fortran90 subroutine in a module

Parameters

[in]	x	the name of the module
[in]	y	the name of the subroutine

#define CALL_FSUBROUTINE

(

x

)

Macro for calling an external fortran subroutine

Parameters

[in]

x

the name of the subroutine

#define END_EXTERN_FORTRAN

#define EXTERN_FFUNCTION

(

x

)

Macro for declaration of an external fortran function

Parameters

[in]

x

the name of the function

#define EXTERN_FMODFUNCTION	(		x,
			y
	)

Macro for declaration of an external fortran90 function in a module

Parameters

[in]	x	the name of the module
[in]	y	the name of the function

#define EXTERN_FMODSUBROUTINE	(		x,
			y
	)

Macro for declaration of an external fortran90 subroutine in a module

Parameters

[in]	x	the name of the module
[in]	y	the name of the subroutine

#define EXTERN_FSUBROUTINE

(

x

)

Macro for declaration of an external fortran subroutine

Parameters

[in]

x

the name of the subroutine

#define FCHARACTER

#define FINT

Macro for fortran INT type in C++ (int)

#define FINT2

Macro for fortran INT*2 type in C++ (short int)

#define FINT8

Macro for fortran INT*8 type in C++ (long int)

#define FLOGICAL

#define FLOGICAL1

#define FREAL

Macro for fortran REAL type in C++ (float)

#define FREAL16

#define FREAL8

Macro for fortran REAL type in C++ (double)

#define FSTRING

#define FSTRINGLEN

#define STD2FSTRING

(

str

)

#define STD2FSTRINGFULL

(

str

)

#define STD2FSTRINGLEN

(

str

)