Архив: Перевести программу с Fortran на Си
Бюджет
10$
/ 600
руб
в месяц
Создан: 5 лет назад
Закрыт
- Описание
- Перевод небольшой программы на си
integer lijp1
real mu/1E-3/
dimension pij(0:100,0:100),vxij(0:100,0:100),vyij(0:100,0:100),vx0ij(0:100,0:100),vy0ij(0:100,0:100), &
aPij(0:100,0:100),aEij(0:100,0:100),aWij(0:100,0:100),aNij(0:100,0:100),aSij(0:100,0:100), &
bPij(0:100,0:100),bEij(0:100,0:100),bWij(0:100,0:100),bNij(0:100,0:100),bSij(0:100,0:100)
data Nx/50/,Ny/50/,dx/0.1/,dy/0.1/,ro/1./,eps/1E-5/
hx=dx;hy=dy;omega=0.5
lijp1=hx/ro;dip1j=hy/ro;delt=1E20
open(1,file='dan2.txt')
vx0ij=0.;vy0ij=0.
vxij=5.;vyij=5.
do while (delt>eps)
delt=0.
do i=0,Nx
do j=0,Ny
pij(i,j)=300.+1.*float(i)+2.*float(j)
end do
end do
!write(1,*) ' pij'
!write(1,1001) ((i,j,pij(i,j),j=1,Ny),i=1,Nx)
1001 format(2i5,g15.7)
do i=0,Nx-1
do j=1,Ny
aPij(i,j)=-2.*mu/hx**2*vxij(i+1,j);aEij(i,j)=mu/hx**2*vxij(i+2,j);aWij(i,j)=(mu/hx**2+mu/hy**2)*vxij (i,j);aNij(i,j)=-2.*mu/hy**2*vxij(i,j+1);aSij(i,j)=mu/hy**2*vxij(i,j+2)
aPij(i,j)=aPij(i,j)*hx*hy;aEij(i,j)=aEij(i,j)*hx*hy;aWij(i,j)=aWij(i,j)*hx*hy;aNij(i,j)=aNij(i,j)*hx *hy;aSij(i,j)=aSij(i,j)*hx*hy
vxij(i,j)=(omega*(-aPij(i,j)-aEij(i,j)-aNij(i,j)-aSij(i,j)-dip1j*(pij(i,j)-pij(i+1,j)))+(1.-omega)*a Wij(i,j))/(mu/hx**2+mu/hy**2)
delt=delt+(vxij(i,j)-vx0ij(i,j))**2
end do
end do
do i=1,Nx
do j=0,Ny-1
bPij(i,j)=-2.*mu/hy**2*vyij(i,j+1);bEij(i,j)=mu/hy**2*vyij(i,j+2);bWij(i,j)=(mu/hx**2+mu/hy**2)*vyij (i,j);bNij(i,j)=-2.*mu/hx**2*vyij(i+1,j);bSij(i,j)=mu/hx**2*vyij(i+2,j)
bPij(i,j)=bPij(i,j)*hy*hx;bEij(i,j)=bEij(i,j)*hy*hx;bWij(i,j)=bWij(i,j)*hy*hx;bNij(i,j)=bNij(i,j)*hy *hx;bSij(i,j)=bSij(i,j)*hy*hx
vyij(i,j)=(omega*(-bPij(i,j)-bEij(i,j)-bNij(i,j)-bSij(i,j)-lijp1*(pij(i,j)-pij(i+1,j)))+(1.-omega)*b Wij(i,j))/(mu/hx**2+mu/hy**2)
delt=delt+(vyij(i,j)-vy0ij(i,j))**2
end do
end do
vx0ij=vxij;vy0ij=vyij
end do !while
write(1,*) 'vx'
write(1,1001) ((i,j,vxij(i,j),j=1,Ny),i=1,Nx)
write(1,*) 'vy'
write(1,1001) ((i,j,vyij(i,j),j=1,Ny),i=1,Nx)
do i=0,Nx-1
do j=1,Ny
aPij(i,j)=-2.*mu/hx**2*vxij(i+1,j);aEij(i,j)=mu/hx**2*vxij(i+2,j);aWij(i,j)=(mu/hx**2+mu/hy**2)*vxij (i,j);aNij(i,j)=-2.*mu/hy**2*vxij(i,j+1);aSij(i,j)=mu/hy**2*vxij(i,j+2)
aPij(i,j)=aPij(i,j)*hx*hy;aEij(i,j)=aEij(i,j)*hx*hy;aWij(i,j)=aWij(i,j)*hx*hy;aNij(i,j)=aNij(i,j)*hx *hy;aSij(i,j)=aSij(i,j)*hx*hy
vxij(i,j)=(omega*(-aPij(i,j)-aEij(i,j)-aNij(i,j)-aSij(i,j)-dip1j*(pij(i,j)-pij(i+1,j)))+(1.-omega)*a Wij(i,j))/(mu/hx**2+mu/hy**2)
end do
end do
do i=1,Nx
do j=0,Ny-1
bPij(i,j)=-2.*mu/hy**2*vyij(i,j+1);bEij(i,j)=mu/hy**2*vyij(i,j+2);bWij(i,j)=(mu/hx**2+mu/hy**2)*vyij (i,j);bNij(i,j)=-2.*mu/hx**2*vyij(i+1,j);bSij(i,j)=mu/hx**2*vyij(i+2,j)
bPij(i,j)=bPij(i,j)*hy*hx;bEij(i,j)=bEij(i,j)*hy*hx;bWij(i,j)=bWij(i,j)*hy*hx;bNij(i,j)=bNij(i,j)*hy *hx;bSij(i,j)=bSij(i,j)*hy*hx
vyij(i,j)=(omega*(-bPij(i,j)-bEij(i,j)-bNij(i,j)-bSij(i,j)-lijp1*(pij(i,j)-pij(i+1,j)))+(1.-omega)*b Wij(i,j))/(mu/hx**2+mu/hy**2)
end do
end do
do ii=1,13
write(1,*) 'Sum_a'
do i=0,Nx-1
do j=1,Ny
aPij(i,j)=-2.*mu/hx**2*vxij(i+1,j);aEij(i,j)=mu/hx**2*vxij(i+2,j);aWij(i,j)=(mu/hx**2+mu/hy**2)*vxij (i,j);aNij(i,j)=-2.*mu/hy**2*vxij(i,j+1);aSij(i,j)=mu/hy**2*vxij(i,j+2)
aPij(i,j)=aPij(i,j)*hx*hy;aEij(i,j)=aEij(i,j)*hx*hy;aWij(i,j)=aWij(i,j)*hx*hy;aNij(i,j)=aNij(i,j)*hx *hy;aSij(i,j)=aSij(i,j)*hx*hy
vxij(i,j)=(omega*(-aPij(i,j)-aEij(i,j)-aNij(i,j)-aSij(i,j)-dip1j*(pij(i,j)-pij(i+1,j)))+(1.-omega)*a Wij(i,j))/(mu/hx**2+mu/hy**2)
write(1,1001) i,j,aPij(i+1,j)+aEij(i+2,j)+aWij(i,j)+aNij(i+1,j+1)+aSij(i+1,j-1)+dip1j*(pij(i,j)-pij(i+1,j))
end do
end do
write(1,*) 'Sum_b'
do i=1,Nx
do j=0,Ny-1
bPij(i,j)=-2.*mu/hy**2*vyij(i,j+1);bEij(i,j)=mu/hy**2*vyij(i,j+2);bWij(i,j)=(mu/hx**2+mu/hy**2)*vyij (i,j);bNij(i,j)=-2.*mu/hx**2*vyij(i+1,j);bSij(i,j)=mu/hx**2*vyij(i+2,j)
bPij(i,j)=bPij(i,j)*hy*hx;bEij(i,j)=bEij(i,j)*hy*hx;bWij(i,j)=bWij(i,j)*hy*hx;bNij(i,j)=bNij(i,j)*hy *hx;bSij(i,j)=bSij(i,j)*hy*hx
vyij(i,j)=(omega*(-bPij(i,j)-bEij(i,j)-bNij(i,j)-bSij(i,j)-lijp1*(pij(i,j)-pij(i+1,j)))+(1.-omega)*b Wij(i,j))/(mu/hx**2+mu/hy**2)
write(1,1001) i,j,bPij(i,j+1)+bEij(i+1,j+1)+bWij(i-1,j+1)+bNij(i,j+2)+bSij(i,j)+lijp1*(pij(i,j)-pij(i,j+1))
end do
end do
end do !ii
close(1)
end
- Категория
Похожие проекты
0
3
Необходимо разработать систему, которая: Подключается к API Diagnocat и iDent (через токены доступа). Получает данные о пациенте, результатах диагностики, снимках и назначениях. Автоматически формирует структурированный отчет по шаблону (включая врачей, этапы лечения, стоимость и описание процедур).......
2 месяца назад
0
0
Необходимо связать БД 1C:ERP с БД Outlook, MS Access, Битрикс24. Например при вводе нового сотрудника в 1С:ERP он должен автоматически появляться в вышеперечисленных БД....
3 месяца назад
0
1
Требуется реализация отчёта на кастомизированной 1С:УПП 1. 2, ТЗ и контрольный пример во вложении, 80% кода должно быть покрыто модульными тестами....
3 месяца назад
$7
0
2
Мне нужно копировать мой текст и вставлять в одну группу ВК с открытой стеной для публикации рекламы. Задача скрипта/сервиса/программы - примерно каждые 10 секунд брать мой текст и вставлять в эту группу. Я пробовал через одну программу, которая записывает действия мышки и потом повторяет, но она как-то......
5 месяцев назад
0
0
Написать прошивку для A94B114, используя 2 подобные рабочие прошивки. Все файлы есть....
5 месяцев назад
$21
0
2
Нужно разархивировать файл архив. ujems, есть распаковщик, но выдает ошибку. Прислать пруфы(скрины)....
7 месяцев назад
0
1
Входные: Имеется скрипт по обработке email писем и добавление их в систему Redmine Проблема: Система обработки этих писем написана некорректная и негибкая Задача: Переписать логику обработки входящих писем, чтобы исключить их потерю при незначительных отклонениях от ожидаемой структуры....
7 месяцев назад
$37
0
1
- В обязанности сотрудника входит поддержание мобильного приложения системы автоматизации бизнеса. - Работа сдельная: задачи выдаются на оценку, далее после согласования передаются на реализацию. - Заключаем договор и NDA (соглашение о неразглашении)....
7 месяцев назад
$71
0
0
Информация о проекте: Проект включает в себя React Native приложение с Expo внутри которого находится webview с сайтом На устройствах с ОС IOS проблема с воспроизведением видео....
8 месяцев назад
0
0
Разрабатываю десктопные программы на языке Java, а так же Java + SQL....
8 месяцев назад
0
0
0
6
1. Чтобы можно было писать описание события (без заголовка) - превью в списке берется N количества букв с начала. (в списке писать дату события и сколько осталось до него) . Выставление даты и времени напоминания и ИНТЕРВАЛА ПОВТОРА ОТ 1 МИНУТЫ! Выставлять количество повторов, далее напоминания прекращаются.......
9 месяцев назад
0
4
Сайт arg2028t. beget. tech/. Нужно сделать вот такой калькулятор. airprint. by/outdoor/wide-format. и еще некоторые доработки....
9 месяцев назад
$214
0
1
0
7
Краткое описание: Мы ищем исполнителя для тестирования функционала мобильного приложения на платформе Android. Задача включает проверку основных функций приложения, таких как работа комнаты, передача видео, звука и чата, а также создание подробной отчетности о найденных ошибках. Требования к исполнителю:......
10 месяцев назад
0
12
$143
0
0
Oпpocы (1008pyб/чaс) Заходитe на сaйт: gonsù....
11 месяцев назад
0
4
Требуется создать скрипт для premiere pro. Данная суть скрипта: в выделеном бине нужно что бы автоматически проставлялся in out ровно по серидине файла так что бы выделенная область была длиной 3 секунды. И требуется кнопка что бы работало через окно expressions. Проблема текущего в том что он не работает......
11 месяцев назад
0
2
Нужно создать приложения под гемблинг/беттинг, на котлин. Приложение должно работать Firebase. Детали вышлю при обсуждение....
11 месяцев назад
$21
0
4
Нужно разархивировать файл archive. ujems, есть распаковщик, но выдает ошибку. Прислать пруфы ....
1 год назад
0
3
Разработать впн сервис ( Приложения на ios и android + сайт для оплаты подписки). Настройка Серверов. Оплата по договоренности. Работаем черед безопасную сделку....
1 год назад
Похожая удалённая работа
- Удаленная работа для веб-программиста
- Удаленная работа для разработчика баз данных
- Удаленная работа для прикладного программиста
- Удаленная работа для системного программиста
- Удаленная работа для разработчика игр
- Удаленная работа для разработчика мобильных приложений
- Удаленная работа для тестировщика ПО
- Удаленная работа для 1С программиста
- Удаленная работа для разработчика встраиваемых систем
- Удаленная работа для разработчика CRM и ERP
Перетащите