Анульовано

Convert a working Matlab to Fortran

Few lines of Matlab code need to be converted to Fortran. As you will see the lines are simple statements

## Deliverables

Matlab code

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

%C! Boundary Conditions - B.C.

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

Er=0; % Radial Volt

sigmar=10; % Radial Stress

sigmart=0; % Shear Stress

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

%C! This is to read properly the knowns or inputs

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

b=cell(1,order);

A=cell(1,order);

TM=cell(1,order);

x=zeros(NoEq,NoEq);

% Applying concentrated-point (delta loads) @ theta0

theta0=pi/4;

%c!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

%c! Fourier Expansion

%c! Eq1, Eq2 or u-r, sigma-r are multiplied by cos

%C! Eq3, Eq4 or u-theta, sigma-r-theta are multiplied by sin

%C!

%c!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

for n=1:order

b{n}(1,1)=(1/pi)*(-e33)*Er*cos(n*theta0);%Eq1 u-r r=Radial

b{n}(2,1)=(1/pi)*(-e33)*Er*cos(n*theta0);%Eq2 sigma-r

b{n}(3,1)=0; %Eq3= is always zero u-theta

b{n}(4,1)=(1/pi)*e33s*Er*sin(n*theta0); %Eq4 sigma-r-theta

b{n}(5,1)=(1/pi)*sigmar*cos(n*theta0); % Radial Stress Boundary

b{n}(6,1)=(1/pi)*sigmart*sin(n*theta0); % Shear Stress Boundary

end

%c!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

%c! Fourier Expansion

%c! Eq1, Eq2 or u-r, sigma-r are multiplied by cos

%C! Eq3, Eq4 or u-theta, sigma-r-theta are multiplied by sin

%C!

%c!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

for n=1:order

%C! n can not be n=0 as t=pi/2*n

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

%C!

%C! Layer 0

%C!

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

%C Evaluating A unknown coeficients in layer 0

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

n0=0;n1=2;

t=0;% cosines coefficients in radial while sines=0

x(m+1,n0+1)=Eq1_0s(n,t,gm,dyna0,A0,prp0,prpp);

% Equation 2 in sigma-theta 0 layer

x(m+2,n0+1)=Eq2_0s(n,t,gm,dyna0,A0,prp0,prpp);

t=pi/(2*n);% cosines coefficients in radial while sines=0

% Equation 3 in sigma r-theta 0 layer

x(m+3,n0+1)=Eq3_0s(n,t,gm,dyna0,A0,prp0,prpp);

% Equation 4 in Dr 0 layer

x(m+4,n0+1)=Eq4_0s(n,t,gm,dyna0,A0,prp0,prpp);

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

%C Evaluating F coeficients in layer 0

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

% Equation 1 in sigma-r 0 layer

t=0;% cosines coefficients in radial while sines=0

x(m+1,n0+2)=Eq1_0s(n,t,gm,dyna0,F0,prp0,prpp);

% Equation 2 in sigma-theta 0 layer

x(m+2,n0+2)=Eq2_0s(n,t,gm,dyna0,F0,prp0,prpp);

% Equation 3 in sigma r-theta 0 layer

t=pi/(2*n);% cosines coefficients in radial while sines=0

x(m+3,n0+2)=Eq3_0s(n,t,gm,dyna0,F0,prp0,prpp);

% Equation 4 in Dr 0 layer

x(m+4,n0+2)=Eq4_0s(n,t,gm,dyna0,F0,prp0,prpp);

m=0;n0=0;n1=2;

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

%C! Layer 1

%C!

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

%C! Physical unknown A

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

% Equation 1 in sigma-r general layer

t=0;% cosines coefficients in radial while sines=0

x(m+1,n1+1)=Eq1_1s(n,t,gm,dyna1,A1,prp1,prpp);

% Equation 2 in sigma-theta general layer

x(m+2,n1+1)=Eq2_1s(n,t,gm,dyna1,A1,prp1,prpp);

% Equation 3 in sigma r-theta general layer

t=pi/(2*n);% cosines coefficients in radial while sines=0

x(m+3,n1+1)=Eq3_1s(n,t,gm,dyna1,A1,prp1,prpp);

% Equation 4 in Dr general layer

x(m+4,n1+1)=Eq4_1s(n,t,gm,dyna1,A1,prp1,prpp);

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

%C Boundary Conditions (known surface traction)

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

t=0;% cosines coefficients in radial while sines=0

x(m+5,n1+1)=sigmarn(n,t,r2,dyna1,A1,prp1);

t=pi/(2*n);% cosines coefficients in radial while sines=0

x(m+6,n1+1)=sigmartn(n,t,r2,dyna1,A1,prp1);

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

%C! Physical Equations Unknown C

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

% Equation 1 in sigma-r general layer

t=0;% cosines coefficients in radial while sines=0

x(m+1,n1+2)=Eq1_1s(n,t,gm,dyna1,C1,prp1,prpp);

% Equation 2 in sigma-theta general layer

x(m+2,n1+2)=Eq2_1s(n,t,gm,dyna1,C1,prp1,prpp);

% Equation 3 in sigma r-theta general layer

t=pi/(2*n);% cosines coefficients in radial while sines=0

x(m+3,n1+2)=Eq3_1s(n,t,gm,dyna1,C1,prp1,prpp);

% Equation 4 in Dr general layer

x(m+4,n1+2)=Eq4_1s(n,t,gm,dyna1,C1,prp1,prpp);

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

%C Boundary Conditions (known surface traction)

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

t=0;% cosines coefficients in radial while sines=0

x(m+5,n1+2)=sigmarn(n,t,r2,dyna1,C1,prp1);

t=pi/(2*n);% cosines coefficients in radial while sines=0

x(m+6,n1+2)=sigmartn(n,t,r2,dyna1,C1,prp1);

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

%C! Physical Equations Unknowns F

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

% Equation 1 in sigma-r general layer

t=0;% cosines coefficients in radial while sines=0

x(m+1,n1+3)=Eq1_1s(n,t,gm,dyna1,F1,prp1,prpp);

% Equation 2 in sigma-theta general layer

x(m+2,n1+3)=Eq2_1s(n,t,gm,dyna1,F1,prp1,prpp);

% Equation 3 in sigma r-theta general layer

t=pi/(2*n);% cosines coefficients in radial while sines=0

x(m+3,n1+3)=Eq3_1s(n,t,gm,dyna1,F1,prp1,prpp);

% Equation 4 in Dr general layer

x(m+4,n1+3)=Eq4_1s(n,t,gm,dyna1,F1,prp1,prpp);

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

%C Boundary Conditions (known surface traction)

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

t=0;% cosines coefficients in radial while sines=0

x(m+5,n1+3)=sigmarn(n,t,r2,dyna1,F1,prp1);

t=pi/(2*n);% cosines coefficients in radial while sines=0

x(m+6,n1+3)=sigmartn(n,t,r2,dyna1,F1,prp1);

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

%C! Physical Equations Unknown H

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

% Equation 1 in sigma-r general layer

t=0;% cosines coefficients in radial while sines=0

x(m+1,n1+4)=Eq1_1s(n,t,gm,dyna1,H1,prp1,prpp);

% Equation 2 in sigma-theta general layer

x(m+2,n1+4)=Eq2_1s(n,t,gm,dyna1,H1,prp1,prpp);

% Equation 3 in sigma r-theta general layer

t=pi/(2*n);% cosines coefficients in radial while sines=0

x(m+3,n1+4)=Eq3_1s(n,t,gm,dyna1,H1,prp1,prpp);

% Equation 4 in Dr general layer

x(m+4,n1+4)=Eq4_1s(n,t,gm,dyna1,H1,prp1,prpp);

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

%C Boundary Conditions (known surface traction)

%C!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!!

t=0;% cosines coefficients in radial while sines=0

x(m+5,n1+4)=sigmarn(n,t,r2,dyna1,H1,prp1);

t=pi/(2*n);% cosines coefficients in radial while sines=0

x(m+6,n1+4)=sigmartn(n,t,r2,dyna1,H1,prp1);

TM{n}=x;

A{n}=TM{n}\b{n};

end

* * *This broadcast message was sent to all bidders on Wednesday Sep 28, 2011 5:10:04 PM:

As some of you have noticed the program calls functions and subrountines. I already have these functions or subroutines in Fortran. The Matlab code is a development on previous code that I had it translated previously. This is an on going work. I code the programs in Matlab then translate them.

* * *This broadcast message was sent to all bidders on Monday Oct 3, 2011 3:45:49 PM:

I used to code my programs then convert them to Fortran. Unfortunately, the same vworker is busy and unavailable. This is an on going project. I would like to pick someone who is a potential future worker. This will ease things for both of us.

Навички: Техніка, Архітектура ПЗ, Тестування ПЗ

Деталі: matlab fortran, vworker end, r.f. engineering, r architecture, my working, m.r.i. software, matlab architecture, equations work, cos software, convert fortran to matlab, a. m. general, translate matlab fortran, matlab all, architecture programs, a0, VOlt, stress a, sin, sin cos, matlab program

Про роботодавця:
( 105 відгуки(-ів) ) Orlando, Canada

ID Проекту: #3603634

4 фрілансерів(-и) у середньому готові виконати цю роботу за $69

koles23

See private message.

$50.15 USD за 10 дні(-в)
(17 відгуків(и))
6.0
matfizvw

See private message.

$90 USD за 10 дні(-в)
(46 відгуків(и))
5.9
spiderdevelopers

See private message.

$100.3 USD за 10 дні(-в)
(16 відгуків(и))
4.3
hervebags

See private message.

$34 USD за 10 дні(-в)
(4 відгуків(и))
3.4