How can I change a FOTRAN code into a MATLAB code withour re-writing the code all over?

I will add the following to my earlier answer, just in case the MEX file approach isn't clear. MEX is short for MATLAB Executable. A MEX acts as a wrapper (or glue or gateway or interface) between MATLAB and existing FORTRAN (or C/C++) code. The MEX file approach basically consist of writing a gateway procedure that calls the existing FORTRAN subroutine or function.

MATLAB sees a MEX as a MATLAB function. You pass it arguments in MATLAB and it returns arguments to MATLAB, just as any MATLAB function does. The MEX gateway function receives the arguments as MATLAB objects, which it processes and passes the required information to the FORTRAN subroutine or function. The return values produced by the FORTRAN subroutine or function are processed by the gateway subroutine into MATLAB objects, which are then passed to MATLAB.

MEX files may be a bit cumbersome to write. But, in my opinion, they are a much better way to make existing FORTRAN code work with MATLAB. Much better than porting the FORTRAN code to MATLAB code using third party software, with results that are not guaranteed to work as they should. Especially if the code to be ported is long and complex and requires changes before porting.

Here's further information:

http://en.wikipedia.org/wiki/MEX_file

http://www.mathworks.se/help/matlab/create-mex-files.html

Vineet Kumar Chaurasiya

please convert this code into matlab

SUBROUTINE 50LV5 COMMON IINPT11 WW(SO),ALFAO,ALFA1 COMMON IBLCOI NP,NPT,NX,NXT,NTR,IT COMMON IGROI X(SO),ETA(S1),DETA(S1),A(S1) COMMON IBLC11 F(S1,2),U(S1,2),V(S1,2),B(S1,2),G(S1,2),P(S1,2), C(S1,2),D(S1,2),E(S1,2),RMU(S1),BC(S1) COMMON IBLC31 DELF(S1),DELU(S1),DELV(S1),DELG(S1),DELP(S1) COMMON/BLCSI 51(S1),52(S1),53(S1),54(S1),55(S1),5S(S1),57(S1), 58(S1),B1(S1),B2(S1),B3(S1),B4(S1),B5(S1),BS(S1), B7(S1),B8(S1),B9(S1),B10(S1),R(5,S1) DIMENSION A11(S1),A12(S1),A13(S1),A14(61),A15(61),A21(61), A22(61),A23(61),A24(S1),A25(61),A31(61),A32(61), A33(61),A34(61),A35(61),G11(61),G12(61),G13(61), G14(61),G15(61),G21(S1),G22(61),G23(61),G24(61), G25(61),G31(61),G32(61),G33(61),G34(61),G35(S1), W1 (61) ,W2 (61) ,W3 (61) ,W4 (61), WS (61) - - - - - - - - - - - - - - - - - -C - C - - - - - ELEMENTS OF TRIANGLE MATRIX,

A11 (1) =1.0 A12(1)=0.0 A13 (1) =0.0 A14 (1) =0. 0 A15(1)=0.0 A21(1)=0.0 A22 (1) =1. 0 A23(1)=0.0 A24 (1) =0. 0 A25(1)=0.0 A31 (1) =0.0 A32(1)=0.0 A33(1)=0.0 A34(1)=ALFAO A35(1)=ALFA1 C ELEMENTS 0 W-VECTOR, ACCORDING TO EO. (13.32A) W1 (1)=R(1,1) W2 (1) =R (2, 1) 414 Finite-Difference Solution of Boundary-Layer Equations W3 (1)=R(3.1) W4 (1) =R (4.1) WS (1) =R (S.l) C - FORWARD SWEEP C DEFINITIONS DO 30 J=2.NP AA1=A(J)*A24(J-l)-A2S(J-l) AA2=A(J)*A34(J-l)-A3S(J-l) AA3=A(J)*A12(J-l)-A13(J-l) AA4=A(J)*A22(J-l)-A23(J-l) AAS=A(J)*A32(J-l)-A33(J-l) AAS=A(J)*A14(J-l)-A1S(J-l) AA7=A(J)ISS(J)-S2(J) AA8=S8 (J) *A (J) AAS=A(J)*BS(J)-Bl0(J) AA10=A(J)*B8(J)-B2(J) C ELEMENTS OF TRIANGLE MATRIX. ACCORDING TO EO. (13.31B) DET=Al1(J-l)*(AA4*AA2-AA1*AAS)-A21 (J-l)* (AA3*AA2-AAS*AAS)+ 1 A31(J-l)*(AA3*AA1-AA4*AAS) Gll(J) =(-(AA4*AA2-AAS*AA1) +A(J)**2* (A21 (J-l)*AA2-A31 (J -1)*AA1))1 1 DET G12(J)=((AA3*AA2-AAS1AAS)-A(J)**2*(All(J-l)*AA2-A3l(J-l)*AAS))/DETG13(J) =(-(AA3*AA1-AA4*AAS) +A(J)**2* (All (J-l)*AA1-A21 (J-l)*AAS))1 1 DET G14(J)=Gl1(J)*A12(J-l)+G12(J)*A22(J-l)+G13(J)*A32(J-l)+A(J) G1S(J)=Gll(J)*A14(J-l)+G12(J)*A24(J-l)+G13(J)*A34(J-l) G21(J)=(S4(J)*(AA2*AA4-AA1*AAS)+A3l(J-l)*(AA1*AA7-AA4*AA8)+ 1 A21(J-l)*(AAS*AA8-AA7*AA2))/DET G22(J)=(All(J-l)*(AA2*AA7-AAS*AA8)+A3l(J-l)*(AA3*AA8-AAS*AA7)+ 1 S4(J)*(AAS*AAS-AA2*AA3))/DET G23(J)=(All (J-l)*(AA4*AA8-AA1*AA7) +S4(J)* (AA3*AA1-AA4* AAS) +1 A2l(J-l)*(AA7*AAS-AA3*AA8))/DET G24(J)=G2l (J)*A12(J-l) +G22(J)*A22(J-l) +G23(J)*A32(J-l) -SS(J) G2S(J)=G2l(J)*A14(J-l)+G22(J)*A24(J-l)+G23(J)*A34(J-l)-S8(J) G3l(J)=(B4(J)*(AA4*AA2-AAS*AA1)-AAS*(A2l(J-l)*AA2-A3l(J-l)*AA1)+ 1 AA10* (A2l (J-l)*AAS-A3l (J-l)*AA4))/DET G32(J)=(-B4(J)*(AA3*AA2-AA5*AAS)+AAS*(All(J-l)*AA2-A3l(J-l)aAAS)-1 AA1D* (All (J-l)aAA5-A3l (J-l)*AA3))/DET G33(J)=(B4(J)*(AA3*AA1-AA4*AAS)-AAS*(Al1(J-l)*AA1-A2l(J-l)*AAS)+ 1 AA10* (All (J-l)*AA4-A2l (J-l)*AA3))/DET G34(J)=G3l(J)*A12(J-l)+G32(J)*A22(J-l)+G33(J)*A32(J-l)-BS(J) G35(J)=G3l(J)*A14(J-l)+G32(J)*A24(J-l)+G33(J)*A34(J-l)-B8(J) C ELEMENTS OF TRIANGLE MATRIX. ACCORDING TO EO. (13.31C) A11 (J) =1.0 A12 (J) =-A (J) -G14 (J) A13(J)=A(J)*G14(J) A14 (J) =-G1S (J) A1S(J)=A(J)*G1S(J) A2l (J) =S3 (J) A22(J)=SS(J)-G24(J) A23(J)=Sl(J)+A(J)*G24(J) A24(J)=-G2S(J)+S7(J) A2S(J)=A(J)*G25(J) A3l (J) =B3 (J) A32(J)=B5(J)-G34(J) A33(J)=B9(J)+A(J)*G34(J) A34(J)=B7(J)-G35(J) A3S (J) =Bl (J) + A(J) *G35 (J) C ELEMENTS OF W-VECTOR. ACCORDING TO EO. (13. 32B) Wl(J) =R(1.J)-Gl1(J)*Wl(J-l)-G12(J)*W2(J-l)-G13(J)*W3(J-l)- 1 G14(J)*W4(J-l)-G1S(J)*W5(J-l) W2(J) =R(2.J)-G21(J)*Wl(J-l)-G22(J)*W2(J-l)-G23(J)*W3(J-l)- 1 G24(J)*W4(J-l)-G2S(J)*WS(J-l) W3(J) =R(3.J)-G31(J)*Wl(J-l)-G32(J)*W2(J-l)-G33(J)*W3(J-l)-1 G34(J)*W4(J-l)-G35(J)*WS(J-l) W4 (J) =R (4. J) WS (J) =R (S. J) -3D CONTINUE C - BACKWARD SWEEP J =NPC DEFINITIONSDP -(A31(J)*(A13(J)*W2(J)-Wl(J)*A23(J))-A32(J)*(All(J)*1 W2 (J) -Wl (J) *A21 (J)) + W3 (J) * (All (J) *A23 (J) -A13 (J) *A2l (J)))DV -(A31(J)*(Wl(J)*A25(J)-W2(J)*A1S(J))-W3(J)*(All(J)*A25(J)1 -A15(J)*A2l (J)) +A35(J)* (All (J)*W2 (J) -Wl (J) *A2l (J)))OF -(W3(J)*(A13(J)*A25(J)-A23(J)*A15(J))-A33(J)*(Wl(J)*A25(J) Fortran Program for Coupled Boundary-Layer Flows 415 -A15 (J) '(W2 (J)) + A35 (J) if (Wl (J) *A23 (J) -A13 (J) *W2 (J))) - (A31 (J) * (A13 (J) *A25 (J) -A23 (J) *A15 (J)) -A33 (J) '( (All (J) '( A25 (J) -A21 (J) '(A15 (J)) +A35 (J) * (All (J) *A23 (J) -A21 (J) *A13 (J)))Dl 1 C ELEMENTS OF DELTA-VECTOR FOR J=NP, ACCORDING TO EO.

DELP (J) DEL V(J) DELF (J) DELG(J) DELU (J) DP ID1 DV ID1 DF ID1 D.O 0.040 J = J-1 C DEFINITIONS BB1=DELU (J+1) -A (J+1) *DELV (J+1) -W4 (J) BB2=DELG(J+1)-A(J+1)*DELP(J+1)-W5(J) CC1=W1(J)-A12(J)'(BB1-A14(J)'(BB2 CC2=W2(J)-A22(J)*BB1-A24(J)*BB2 CC3=W3(J)-A32(J)'(BB1-A34(J)ifBB2 DD1=A13(J)-A12(J)'(A(J+l) DD2=A23(J)-A22(J)*A(J+l) DD3=A33(J)-A32(J)*A(J+l) EE1=A15(J)-A14(J)*A(J+l) EE2=A25(J)-A24(J)*A(J+1) EE3=A35(J)-A34(J)'(A(J+l) DETT=A11 (J)*DD2*EE3+A21 (J)*DD3'(EE1+A31 (J)'(DD1*EE2 1 -A3l (J)*DD2*EE1-A21 (J)*DD1*EE3-A11 (J)*DD3*EE2 C ELEMENTS OF DELTA-VECTOR, ACCORDING TO EO. (13.33B)DELF(J)=(CC1*DD2*EE3+CC2*DD3*EE1+CC3*DD1*EE2-CC3*DD2ifEEl -CC2*DD1*EE3-CC1*DD3*EE2)/DETT DELV(J)=(A11(J)*CC2*EE3+A21(J)*CC3*EE1+A31(J)*CC1*EE2- A31 (J)*CC2*EE1-A21 (J)*CC1 I EE3-Al1 (J)*CC3*EE2)/DETT DELP(J)=(Al~(J)*CC3*DD2+A21(J)*CC1*DD3+A3l(J)*CC2*DD1-1 1 1 A3l(J)*CC1*DD2-A21(J)*CC3*DD1-Al1(J)*CC2*DD3)/DETT DELU(J)=BB1-A(J+l).DELV(J) DELG(J)= BB2-A(J+l)*DELP(J) IF(J .GT. 1) GO TO 40 C NEW VALUES OF F,U,V,G,P, ACCORDING TO EO. (13.61) DO 50 J=l,NP F(J,2)=F(J,2)+DELF(J) U(J, 2) =U (J, 2) +DELU (J) V(J, 2) =V (J, 2) +DELV (J) G(J, 2) =G (J, 2) +DELG (J) P (J, 2) =P (J, 2) +DELP (J) CONTINUE U (1, 2) = O. 0RETURN END

Didier Pinchon

You may be interested by the following url with reference to program f2matlab:

http://engineering.dartmouth.edu/~d30574x/consulting/consultingIndex.html

that I obtained by asking to Google "Fortran to Matlab translator".

Sayantan Mukherjee

There is a program that can convert fortan code to matlab code.You can download from the following link..

http://www.mathworks.com/matlabcentral/fileexchange/5260-f2matlab

It would be a much better idea to write a MEX wrapper for your FORTRAN code. This way there is no need to port anything to MATLAB code. Also, porting FORTRAN code to MATLAB code will most likely make execution much slower.

Pankaj K Mishra

Look at a similar discussion here

http://stackoverflow.com/questions/12116916/converting-fortran-code-to-matlab-references

Olaniyi Samuel Iyiola

Thanks to all

Gbenga Ojo

Please check the site, hope it will be helpful

http://math.hawaii.edu/~dale/190/06l.pdf

Afaq Ahmad

During conversion I have following errors.

Error using fscanf

Invalid file identifier. Use fopen to generate a valid file identifier.

Error in addExtraFiles (line 20)

temp1=fscanf(fid,'%c'); fclose(fid);

Error in f2matlab (line 5201)

addExtraFiles

Please guide me.

Regards

Ben Earl Barrowes

Afaq, I am just reading this, but I thought I would write this in case you are still working on this translation. Something seems to have gone wrong with file permissions when f2matlab is trying to find the file writeFmt.m. Please contact me so I can help you with this one.

Ghouse, Contact me for fortran conversion help: [email protected]

Oluwaseyi Aliu

Ben Earl Barrowes kindly assist with a conversion from fortran to matlab

Luis Eduardo Olivar Robayo

Hi, does anyone have the matlab codes for data mollification ?

Please what is the best method of extraction for health supplement from botanical ?

Are there relationships between the temperature, concrete strength development and the air permeability of hardened concrete?

How can I get "Carbon nanotube-based chemical sensing materials"?

Good day all, please what concentration of extract is best to be used for testing of susceptibility of Salmonella species?

Which source material will be appropriate?

Is there a better (efficient) way to code the coefficient matrix of Riesz space fractional reaction diffusion equation both in 1D & 2D in MATLAB?

Climate change modelling

COVID-19 and Aquatic resources (fish species)

How to compare two groups with only two measurements?

Enhanced Yellow Fluorescent Protein (Aequorea victoria) DNA sequence?

What kind of fluid could I use for a pitfall trap that also does not invalidate molecular testing?

What's the best way to measure growth rates in House sparrow chicks from day 2 to day 10?

Do we need to write protocol for a second study of a systematic review that has been registered and published?

How to increase a model accuracy ?

Do I need to do linear regression before running mediation analysis using PROCESS MACRO by Hayes ?

Sample dilution for quantitative analysis using ELISA?

Is There Any Feasible Method To Test Fluorescent Compounds Other Than Fluorescence Spectrometers ??

How to calculate average trend in maximum annual NDVI from MODIS?