Превод с фортрана на СИ

Модераторы: Hawk, Romeo, Absurd, DeeJayC, WinMain

grigzema
Сообщения: 10
Зарегистрирован: 01 ноя 2009, 21:17

void main()
{
double e[50][1000], f[50][1000], t[50][50][1000],b1,b2,b,b3,b4,b5,sumtau,a1,a2,t1,l122,t2,l111,l112,l121,l211,l222,l221,l212,c,tsr,tw,rrx,rry,rr,tau,r,dy,dx,df,deltax,deltaf1,deltay,deltaf2;
double x1,y1,tl,tlfx,tlfy,nx,mx,ny,my,nnx,nny,ny1,nx1,kfs,j,k,i,bn,rf,mws,kj,mks,ts;


//ввод данных


//выбор шага сетки по Х
a1=4.95E-6;
a2=5.78E-6;
deltax=x1/1000;
deltaf1=tlfx/1000;
dx=sqrt(6*a1*tau);
df=sqrt(4*a2*tau);
nx=int((deltax/dx)+0.5)+1;
mx=int((deltaf1/df)+0.5)+nx+1;
nnx=mx-nx;
printf ("Число узлов в металле по Х=",nx);
printf ("Число узлов в форме по Х=",nnx);

//выбор шага сетки по У
deltaf2=tlfy/1000;
deltay=y1/1000;
dy=sqrt(4*a1*tau);
ny=int((deltay/dy)+0.5)+1;
my=int((deltaf2/df)+0.5)+ny+1;
nny=my-ny;
printf ("Число узлов в металле по Y=",ny);
printf ("Число узлов в форме по Y=",nny);

//задание начальных условий в отливке
k=1;
do
{ j=1;
do
{ i=1;

t[j][k]=t1;
}
while(i=nx);
}
while(i=ny);
grigzema
Сообщения: 10
Зарегистрирован: 01 ноя 2009, 21:17

Люди ну помогите ктонить я уже готов заплатить
atavin-ta
Сообщения: 585
Зарегистрирован: 30 янв 2009, 06:38

А фортранный исходник можешь растолковать? Здесь не аксакалы собрались, а сионисты и в товоём коде мало кто разберётся, а быстро вообще никто. Кстати, с чего ты взял, что это фортарн и его надо переводить? По-моему это уже c.
Вопрос: "Почему вы все сионисты? Нельзя ли писать на чём то другом?".
Ответ: "Писать можно на чём угодно. Но зачем же так себя ограничивать? Пиши на С!".
grigzema
Сообщения: 10
Зарегистрирован: 01 ноя 2009, 21:17

все правильно это уже сишный код, изначально он был написан на фортране. Программа используя неявные конечно разностные уравнения считает коэффициенты уравнений, после чего используя обратную прогонку вычисляются температуры в узлах сетки. массив t отвечает за описание просчитываемой детали, пошагово. Выкладываю весь код если пригодится.
Вложения

[Расширение txt было запрещено, вложение больше недоступно.]

[Расширение txt было запрещено, вложение больше недоступно.]

[Расширение txt было запрещено, вложение больше недоступно.]

[Расширение txt было запрещено, вложение больше недоступно.]

atavin-ta
Сообщения: 585
Зарегистрирован: 30 янв 2009, 06:38

Так ты температуру моделишь? Я на дипломе, на диссертации и на гранте моделил температуру же. Так что если есть конкретные вопросы, давай их сюда.
Вопрос: "Почему вы все сионисты? Нельзя ли писать на чём то другом?".
Ответ: "Писать можно на чём угодно. Но зачем же так себя ограничивать? Пиши на С!".
grigzema
Сообщения: 10
Зарегистрирован: 01 ноя 2009, 21:17

я температуру не моделирую я ее считаю, из уравнения фурье я выделяю температуру и считаю ее согласно ранее посчитаным коэффициентам уравнения. http://www.uchites.ru/files/nummethod_b ... 5-1234.pdf вот тут есть про эту схему.
atavin-ta
Сообщения: 585
Зарегистрирован: 30 янв 2009, 06:38

Вычисление температуры по уравнениям теплопроводности и конвективного теплопереноса - это и есть математическое моедлирования тепловых процессов (или кратнко - моделирование температуры).
Вопрос: "Почему вы все сионисты? Нельзя ли писать на чём то другом?".
Ответ: "Писать можно на чём угодно. Но зачем же так себя ограничивать? Пиши на С!".
grigzema
Сообщения: 10
Зарегистрирован: 01 ноя 2009, 21:17

уговорили)
Аватара пользователя
Romeo
Сообщения: 3126
Зарегистрирован: 02 мар 2004, 17:25
Откуда: Крым, Севастополь
Контактная информация:

atavin-ta, ты вместо того, чтобы спорить о терминологии лучше бы помог человеку, если разбираешься в этом :)
Entites should not be multiplied beyond necessity @ William Occam
---
Для выделения С++ кода используйте конструкцию [ code=cpp ] Код [ /code ] (без пробелов)
---
Сообщение "Спасибо" малоинформативно. Благодарность правильнее высказать, воспользовавшись кнопкой "Reputation" в виде звёздочки, расположенной в левом нижнем углу рамки сообщения.
Ответить