many fortran C mixed programing examples added

Author: Foadsf

.gitignore

@@ -155,3 +155,15 @@ lapacke_mangling_with_flags.h
 lapacke_mangling.h
 lapack_utils.h
 lapacke.h
+
+
+libblas.a
+librefblas.a
+libcblas.a
+libf2c.a
+libgfortran.dylib
+liblapack.a
+liblapacke.a
+liblapacke.dylib
+librefblas.a
+libtmglib.a

A_fortran/other/cfortran/cfortran-20110621.tar.bz2

@@ -20,7 +20,21 @@ sources to be looked at:
 8. https://docs.oracle.com/cd/E19059-01/stud.9/817-6694/11_cfort.html
 9. https://www.math.utah.edu/software/c-with-fortran.html
 10. https://software.intel.com/en-us/forums/intel-visual-fortran-compiler-for-windows/topic/290051
+11. http://luthaf.github.io/calling-c++-from-fortran.html
+12. http://alignment.hep.brandeis.edu/Software/Mixing/Mixing_Manual.html
+13. http://www.fortran.bcs.org/2002/interop.htm
+14. http://www.pgroup.com/userforum/viewtopic.php?p=7398&sid=2a05571a922045c6c4ae26d65bbf9fba
+15. http://www.unidata.ucar.edu/software/netcdf/examples/programs/
+16. http://people.sc.fsu.edu/~jburkardt/c_src/mixed/mixed.html
+17. https://docs.oracle.com/cd/E19422-01/819-3685/11_cfort.html
+
 
 points:
 
 1. name of the routines/subroutines/functions and structs/commons in some versions or dialects of fortran must be in small letters and and an underscore _ must be added to the C function
+
+
+issues:
+1. test22 result not correct
+2. test19 does not compile
+3. test24 not compiling

A_fortran/other/f90ppr/f90ppr.f90.gz

@@ -1,9 +1,11 @@
-this file shows how we can call a C function from fortran
+fortran code written if FORTRAN 77?
+
+this file shows how we can call a C function from fortran and how to pass a string properly from fortran to C
 source: http://www.yolinux.com/TUTORIALS/LinuxTutorialMixingFortranAndC.html
 
 
 points:
-1. when you pass a string from fortran to C you actually pass two variables including an array of single characters and also the length of the array.
+1. when you pass a string from fortran to C you actually pass two variables including a pointer to an array of single characters (string) and also the length of the array.
 2. when you pass a string to C you need to first be sure that you allocate at least one letter more than the characters you need and in C first of se "\0" or NULL as the end of the string
 3. common is the equivalent of struct in C
   in C:
@@ -18,7 +20,9 @@ points:
     ...
     common/mystruct/ i, c, ...
 
-4. when fortran passes a variable to C function, the C function receives a pointer
+4. when fortran passes a variable to C function, the C function receives a pointer?
+5. when you want to call a C function in fortran you don't need to declare a dummy in the fortran code, unlike C ,fortran doesn't care if the function is not declared and then on run time it will look for the function (how about variables and structs?)
+6. when passing variables, functions or srtucts from/to C/fortran it is Case-insensitive. but you should be careful that inside each world it is case sensitive.
 
 
 
@@ -28,3 +32,11 @@ questions:
 I think we needed a
   external doubleijk
   here, but we don't!
+
+
+
+
+fortran points:
+1. in FORTRAN 77 (?) one can print to terminal with something like:
+        write(6,<line number for format e.g. 10>) var_1,var_2,...
+  10    format('some text', typeof_var_N,'some other text',...)

A_fortran/other/floppy7/floppy7.tgz

@@ -1,7 +1,8 @@
       program test
 
       integer ii, jj, kk
-      common/ijk/ ii, jj, kk
+      common/IjK/ ii, jj, kk
+
       real*8  ff
       character*25 cc
       character*7 dd

A_fortran/other/forcheckIIDE/forcheckIDE-14.6-Qt5-source.tar.gz

@@ -0,0 +1,8 @@
+all:
+	gfortran -c ranvec.f90
+	gcc -c main.c
+	gfortran -o result.out ranvec.o main.o
+	rm -rf *.o
+
+clean :
+	rm -rf *.out *~ *.bak *.o

A_fortran/readme.txt

@@ -0,0 +1,14 @@
+#include <stdio.h>
+
+void ranvec_(double *, int *);
+
+int main() {
+  int n = 5;
+  double a[5];
+
+  ranvec_(a, &n);
+
+  printf("a = %f, %f, %f, %f, %f\n", a[0], a[1], a[2], a[3], a[4]);
+
+  return 0;
+}

A_fortran/test1/readme.txt

@@ -0,0 +1,8 @@
+SUBROUTINE ranvec(v, n)
+  INTEGER                        :: n
+  DOUBLE PRECISION, DIMENSION(n) :: v
+
+  CALL RANDOM_SEED()
+  CALL RANDOM_NUMBER(v)
+
+END SUBROUTINE ranvec

A_fortran/test1/test1F.f

@@ -0,0 +1 @@
+source: http://www.math.chalmers.se/Math/Grundutb/CTH/tma881/0809/Assignments/Mixing_C_Fortran.html

A_fortran/test10/Makefile

@@ -0,0 +1,8 @@
+all:
+	gfortran -c add.f
+	gcc -c test11.c
+	gfortran -o result.out add.o test11.o
+	rm -rf *.o
+
+clean :
+	rm -rf *.out *~ *.bak *.o

A_fortran/test10/main.c

@@ -0,0 +1,7 @@
+      REAL*8 FUNCTION ADD (A, B, C, D)
+      REAL*8 B,D
+      INTEGER*4 A,C
+      DIMENSION B(4), D(4)
+      ADD = B(A) + D(C)
+      RETURN
+      END

A_fortran/test10/ranvec.f90

@@ -0,0 +1 @@
+source: http://www.ibm.com/support/knowledgecenter/SSGH2K_13.1.2/com.ibm.xlc131.aix.doc/proguide/tucsampl.html

A_fortran/test10/readme.txt

@@ -0,0 +1,21 @@
+#include <stdio.h>
+
+double add_(int *, double[], int *, double[]);
+
+double ar1[4] = {1.0, 2.0, 3.0, 4.0};
+double ar2[4] = {5.0, 6.0, 7.0, 8.0};
+
+int main() {
+  int x, y;
+  double z;
+
+  x = 3;
+  y = 3;
+
+  z = add_(&x, ar1, &y, ar2); /* Call Fortran add routine */
+  /* Note: Fortran indexes arrays 1..n */
+  /* C indexes arrays 0..(n-1) */
+
+  printf("The sum of %1.0f and %1.0f is %2.0f \n", ar1[x - 1], ar2[y - 1], z);
+  return 0;
+}

A_fortran/test11/Makefile

@@ -0,0 +1,56 @@
+# CC = gcc
+# FC = gfortran
+# # Make sure the CC and FC compilers were built with the same glibc library
+# # g77 will work fine on most systems
+# CFLAGS = -Wall
+#
+# OBJS = testdot.o testret.o testmat.o
+#
+# all: driver lib
+#
+# driver : driver.o $(OBJS)
+# 	$(FC) -o driver driver.o $(OBJS) -lstdc++
+#
+# testret.o : testret.cpp
+# 	$(CC) $(CFLAGS) -c testret.cpp
+#
+# testdot.o : testdot.cpp
+# 	$(CC) $(CFLAGS) -c testdot.cpp
+#
+# testmat.o : testmat.cpp
+# 	$(CC) $(CFLAGS) -c testmat.cpp
+#
+# driver.o : driver.f
+# 	$(FC) -c driver.f
+#
+# # Default Targets for Cleaning up the Environment
+# clean :
+# 	rm *.o
+# 	rm *.a
+#
+# pristine :
+# 	rm *.o
+# 	rm *.a
+# 	touch *.cpp
+#
+# ctags :
+# 	ctags *.cpp
+#
+# # Target for making the library
+#
+# lib: $(OBJS)
+# 	ar -rc libbblas.a $(OBJS)
+# 	ranlib libbblas.a
+
+
+all:
+	gfortran -c driver.f
+	g++ -c testdot.cpp
+	g++ -c testmat.cpp
+	g++ -c testret.cpp
+	g++ -o result.out testret.o testmat.o testdot.o driver.o -L /usr/local/Cellar/gcc/6.2.0/lib/gcc/6/ -lgfortran
+	rm -rf *.o
+
+
+clean :
+	rm -rf *.out *~ *.bak *.o

A_fortran/test11/add.f

@@ -0,0 +1,56 @@
+      program driver
+
+      integer veclen matdim
+      real mat, dot
+      real vec1(2), vec2(2)
+      real matrix1(3,3), matrix2(3,3), matrix3(3,3)
+
+      veclen = 2
+      matdim = 3
+
+      vec1(1) = 1.0
+      vec1(2) = 2.0
+      vec2(1) = 3.0
+      vec2(2) = 4.0
+
+      matrix1(1,1) = 1.0
+      matrix1(2,1) = 0.0
+      matrix1(3,1) = 0.0
+      matrix1(1,2) = 0.0
+      matrix1(2,2) = 1.0
+      matrix1(3,2) = 0.0
+      matrix1(1,3) = 0.0
+      matrix1(2,3) = 0.0
+      matrix1(3,3) = 1.0
+
+      matrix2(1,1) = 1.0
+      matrix2(2,1) = 0.0
+      matrix2(3,1) = 0.0
+      matrix2(1,2) = 0.0
+      matrix2(2,2) = 2.0
+      matrix2(3,2) = 0.0
+      matrix2(1,3) = 0.0
+      matrix2(2,3) = 0.0
+      matrix2(3,3) = 3.0
+
+      mat = 0.0;
+
+* mat should get incremented by 1 in testret
+      do 10 i = 1, 10
+         call testret(mat)
+10    continue
+
+* testdot should return the dot product of vec1 and vec2 as a real
+      dot = testdot(veclen, vec1, vec2)
+
+* testmat uses matrix1 and matrix2 to calculate matrix3
+      call testmat(matdim, matrix1, matrix2, matrix3)
+
+      print *, mat, '  ', dot
+      print *,
+      do 20 i = 1, matdim
+        write (6,*) (matrix3(j,i), j=1,matdim)
+20    continue
+
+
+      end

A_fortran/test11/redme.txt

@@ -0,0 +1 @@
+source: http://theochem.mercer.edu/interlanguage/

A_fortran/test11/test11.c

@@ -0,0 +1,19 @@
+
+extern "C" {
+    float testdot_( int *N, float *vec1, float *vec2);
+    }
+
+
+float testdot_(  int *N, float *vec1, float *vec2 ) 
+{
+    int i;
+    float dot;
+
+    dot = 0.0;
+
+    for (i=0; i<*N; i++) { 
+      dot = dot + *(vec1+i) * *(vec2+i);
+      }
+    return dot;
+}
+

A_fortran/test12/Makefile

@@ -0,0 +1,24 @@
+
+extern "C" {
+    void testmat_( int *N, float *mat1, float *mat2, float *mat3);
+    }
+
+
+void testmat_(  int *N, float *mat1, float *mat2, float *mat3 ) 
+{
+    int i,j,k;
+    float dot;
+
+    // Assign stuff to mat3 from mat1 and mat2
+    for (i=0; i<*N; i++) { 
+      for (j=0; j<*N; j++) {
+         dot = 0.0;
+         for (k=0; k<*N; k++ ) {
+           dot = dot + *(mat1+(*N*j+k)) * *(mat2+(*N*i+j));
+         }
+         *(mat3+(*N*j+i)) = dot;
+       }
+    }
+
+}
+

A_fortran/test12/driver.f

@@ -0,0 +1,11 @@
+
+extern "C" {
+    void testret_( float *mat);
+    }
+
+
+void testret_( float *mat ) 
+{
+    *mat = *mat + 1;
+}
+

A_fortran/test12/readme.txt

@@ -0,0 +1,8 @@
+all:
+	gfortran -c test13F.f
+	gcc -c test13C.c
+	gfortran -o result.out test13F.o test13C.o
+	rm -rf *.o
+
+clean :
+	rm -rf *.out *~ *.bak *.o

A_fortran/test12/testdot.cpp

@@ -0,0 +1 @@
+source: http://striky.ece.jhu.edu/~sasha/SOFTWARE/FORTRAN/007-3696-005/sgi_html/ch07.html

A_fortran/test12/testmat.cpp

@@ -0,0 +1,15 @@
+#include <stdio.h>
+#include <string.h>
+
+typedef char fstr_16[16];
+
+void prt_(fstr_16, float *, fstr_16, int, int);
+
+int main() {
+  float val = 2.1828e0;
+  fstr_16 bef, aft;
+  strncpy(bef, "Before..........", sizeof(bef));
+  strncpy(aft, "...........After", sizeof(aft));
+  prt_(bef, &val, aft, sizeof(bef), sizeof(aft));
+  return 0;
+}

A_fortran/test12/testret.cpp

@@ -0,0 +1,6 @@
+      SUBROUTINE PRT(BEF, VAL, AFT)
+      CHARACTER*(*) :: BEF, AFT
+      REAL :: VAL
+      PRINT *, BEF, VAL, AFT
+      RETURN
+      END SUBROUTINE PRT

A_fortran/test13/Makefile

@@ -0,0 +1,8 @@
+all:
+	gfortran -c test13F.f
+	gcc -c test13C.c
+	gfortran -o result.out test13F.o test13C.o
+	rm -rf *.o
+
+clean :
+	rm -rf *.out *~ *.bak *.o

A_fortran/test13/readme.txt

@@ -0,0 +1,16 @@
+source: http://striky.ece.jhu.edu/~sasha/SOFTWARE/FORTRAN/007-3696-005/sgi_html/ch07.html
+
+
+on the fortran side if you do not specify a the type of the variable it assumes it's real!
+
+when transferring a string or an array from C to fortran we must use strlen or size of function to pass the size as well
+
+  void func_(int *, int, char *, int);
+
+  char * mychar = "some text";
+  int * myint = {1,2,3,4};
+
+  func_(myint, sizeof(myint), mychar , strlen(mychar));
+
+
+what surprises me is that somehow miraculously even if we change the order of the variables fortran understands to assign the output of strlen and sizeof functions to the right array/string