Вопросы программирования на FORTRAN (ФОРТРАН) - [328] :: Прикладное программирование :: Компьютерный форум Ru.Board

Перейти из форума на сайт.

Новости • Файловые архивы
Поиск • Активные темы • Топ лист
Правила • Кто в on-line?

Вход • Забыли пароль? • Первый раз на этом сайте? • Регистрация

Компьютерный форум Ru.Board » Компьютеры » Прикладное программирование » Вопросы программирования на FORTRAN (ФОРТРАН)

Модерирует : ShIvADeSt
Версия для печати • Подписаться • Добавить в закладки
Страницы: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330

Открыть новую тему

Написать ответ в эту тему

akaGM

Platinum Member

Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору

Обсуждаются все вопросы, связанные с программированием на ФОРТРАН, как общего так и конкретного характера.
Постарайтесь дать как можно больше информации о возникшей проблеме -- это в конце концов в ваших же интересах чтобы вам помогли...

прежде чем просить помощи в задании
платное решение задач

ресурсы этого топика
ссылка на подборку ресурсов, собранных посетителями этого форума

то, чем мы решили поделиться
ссылка на страничку программ etc собственного изготовления, которыми любезно делятся наши форумчане

если вам вдруг не отвечают или ответ вас не устраивает
и вообще полезно прочитать всем спрашивающим

просьба к пишущим и отвечающим все большие листинги оформлять тегом more
и отключать графические смайлики при размещении фортран-кода

Всего записей: 24120 | Зарегистр. 06-12-2002 | Отправлено: 18:11 14-01-2007 | Исправлено: akaGM, 09:47 01-03-2020

akaGM

Platinum Member

Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору

херовый костыль...
проще/быстрее из Ф С-функу звать...
жаль в Ф асм-вставок нет...

Добавлено:

Цитата:

Преобразовать unsigned

зачем тогда было DWORD / ULONG заводить...

Добавлено:

код-то копеечный, 10 строк...

Код:

unsigned long WR512(void) {
unsigned long a, b, c, d;
a = state[index];
c = state[(index+13)&15];
b = a^c^(a<<16)^(c<<15);
c = state[(index+9)&15];
c ^= (c>>11);
a = state[index] = b^c;
d = a^((a<<5)&0xDA442D24UL);
index = (index + 15)&15;
a = state[index];
state[index] = a^b^d^(a<<2)^(b<<18)^(c<<28);
return state[index];
}

мля, сбой только в одной красной строке!

Добавлено:

c = IXOR(c, SHIFTR(c, 11))

Всего записей: 24120 | Зарегистр. 06-12-2002 | Отправлено: 17:50 18-11-2023 | Исправлено: akaGM, 18:00 18-11-2023

MBK2

Silver Member

Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору

OMFG, это что рандомный полином, зачем такое страшное кодить на фортране?
Если вас это утешит - в JS, например, вообще нет unsigned int32, приходится какие то безумные костыли делать. Еще можно преобразовать число в массив байт (или даже бит) и оперировать с этим массивом, потом преобразовать обратно.

Всего записей: 4576 | Зарегистр. 18-09-2018 | Отправлено: 18:30 18-11-2023 | Исправлено: MBK2, 18:32 18-11-2023

akaGM

Platinum Member

Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору

Цитата:

зачем такое страшное кодить на фортране?

это не для себя,
меня попросили, я делаю...

Добавлено:

/Zl -- рулит!

Всего записей: 24120 | Зарегистр. 06-12-2002 | Отправлено: 18:39 18-11-2023

akaGM

Platinum Member

Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору

Цитата:

Если вас это утешит - в JS, например, вообще нет unsigned int32

так в том-то и дело, что есть!
два типа -- DWORD и ULONG, тока в них знаковый бит всё равно защищён в побитовых операциях...

Всего записей: 24120 | Зарегистр. 06-12-2002 | Отправлено: 16:38 20-11-2023

MBK2 Silver Member	Редактировать \| Профиль \| Сообщение \| Цитировать \| Сообщить модератору Я не пойму у вас получилось сишный код подлинковать или нет?
	Всего записей: 4576 \| Зарегистр. 18-09-2018 \| Отправлено: 18:28 20-11-2023

akaGM

Platinum Member

Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору

Цитата:

Я не пойму у вас получилось сишный код подлинковать или нет?

нет, блин, не получилось, третий день линкую :)
чему здесь получаться, это последнее что бы мне хотелось иметь...

в итоге вбил костыль в самом фортране -- обошёлся одной проверкой на отрицательность в теле функи...

Всего записей: 24120 | Зарегистр. 06-12-2002 | Отправлено: 18:41 20-11-2023 | Исправлено: akaGM, 18:49 20-11-2023

MBK2 Silver Member	Редактировать \| Профиль \| Сообщение \| Цитировать \| Сообщить модератору Да, пожалуй, это самый простой костыль.
	Всего записей: 4576 \| Зарегистр. 18-09-2018 \| Отправлено: 19:21 20-11-2023

ne_viens Advanced Member	Редактировать \| Профиль \| Сообщение \| Цитировать \| Сообщить модератору INTEGER*8 c !?
	Всего записей: 1530 \| Зарегистр. 01-11-2004 \| Отправлено: 10:11 27-11-2023

akaGM

Platinum Member

Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору

ne_viens

1) как я уже говорил, это не себе, и как только клиент сказал Большое Спасибо и что он уже вовсю юзает, я об этом просто забыл...

2) оригинальный алгоритм, как я писал выше:

Код:

// return 32 bit random number
WR512(void)

неплохой, кстати, алгоритм

3) зачем чинить несломанное? :)

Добавлено:

сам мой код, вдруг кто решит воспользоваться...
WELL RNG 512bit

Всего записей: 24120 | Зарегистр. 06-12-2002 | Отправлено: 16:03 27-11-2023 | Исправлено: akaGM, 17:08 27-11-2023

MBK2 Silver Member	Редактировать \| Профиль \| Сообщение \| Цитировать \| Сообщить модератору Это неполный код, где инициализация state? Вот как надо
	Всего записей: 4576 \| Зарегистр. 18-09-2018 \| Отправлено: 22:31 27-11-2023

akaGM

Platinum Member

Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору

не знаю как надо, а С-оригинал выглядел вот так:

Код:

subroutine InitWELLRNG512(init)
implicit none

integer(ULONG) DIM_SIZE
parameter (DIM_SIZE = 16)

integer(ULONG) init(0:DIM_SIZE-1)
integer(ULONG) j

state_i = 0 ! global var

do j = 0, DIM_SIZE-1
state(j) = init(j)
enddo

end subroutine InitWELLRNG512

Всего записей: 24120 | Зарегистр. 06-12-2002 | Отправлено: 22:48 27-11-2023

MBK2

Silver Member

Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору

akaGM
Я не пойму, вы стегаетесь или вообще в алгоритм не вникали?
Без этой инициализации рандомный генератор работать не будет, а где она у вас?

Код:

void well512_init( Well512* ws, uint32_t seed )
{
int i;
uint32_t prev;
uint32_t* wstate = ws->wstate;

ws->wi = 0;
wstate[0] = seed;
for( i = 1; i < WELL512_STATE_SIZE; ++i )
{
prev = wstate[i-1];
wstate[i] = (1812433253 * (prev ^ (prev >> 30)) + i);
}
}

Всего записей: 4576 | Зарегистр. 18-09-2018 | Отправлено: 12:51 28-11-2023

akaGM

Platinum Member

Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору

MBK2

Цитата:

Я не пойму, вы стегаетесь

нет, я не садо-мазо...

без какой _этой_ инициализации?

вот мой пример использоваеия:

Код:

integer(ULONG) init(0:DIM_SIZE-1)/1, 634511, 17, 286, 8, 56789,
& 10001, 33888, 1000000, 4000000, 55000000,
& 1111, 99999933, 10, 14, 666/
integer(ULONG) i, j, N/10/
real*8 r

call InitWELLRNG512(init)

do i = 1, N
r = WELLRNG512()
write(*,'(F10.6)') r
enddo

END

как юзверь инициализирует init(), хоть штатным рандомом, хоть как, хоть кверху каком -- не моё собачье дело,
я по ТЗ переписал две функи, всё...
и да, я
Цитата:

в алгоритм не вникали

Всего записей: 24120 | Зарегистр. 06-12-2002 | Отправлено: 14:24 28-11-2023 | Исправлено: akaGM, 14:25 28-11-2023

MBK2

Silver Member

Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору

Ну, оно, то, конечно, можно и так инициализировать, но в итоге получаетса не

Цитата:

неплохой, кстати, алгоритм

,
а какая то колхозная самодеятельность по его мотивам.
Я и так не особо понимаю, за каким хреном его к фортрану прикручивать, там что, своего штатного рандома нет?

Всего записей: 4576 | Зарегистр. 18-09-2018 | Отправлено: 15:53 28-11-2023

akaGM

Platinum Member

Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору

вот я поднёс старушке сумку до подъезда, а ты мне:

1) чё, у старушки внуков нет?
2) а чё ты ей за херню в сумке поднёс, она домой заказать не могла?
етц

оч похоже...
если моя деятельность так уж тебя раздражает, то нафига ко мне цепляться?

Добавлено:

для тех, кто тоже считает алгоритм колхозным:

доки по WELL RNG (1.2 MB)

https://disk.yandex.ru/d/SHdo8RBAdAaleQ

Всего записей: 24120 | Зарегистр. 06-12-2002 | Отправлено: 17:31 28-11-2023

MBK2

Silver Member

Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору

Я ничего против алгоритма не имею, просто инициализация массива state вызывает нарекания.
Вот, к примеру реализация на js, так тут они, хотя бы randomizом инитят. Кстати, обратите внимание, как они выкрутились с отсутствием uint на js.

Всего записей: 4576 | Зарегистр. 18-09-2018 | Отправлено: 17:57 28-11-2023

akaGM

Platinum Member

Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору

Цитата:

инициализация массива state вызывает нарекания

Цитата:

как юзверь инициализирует init(), хоть штатным рандомом, хоть как, хоть кверху каком -- не моё собачье дело

выяснил -- штатным рандомом...
полегчало?

Всего записей: 24120 | Зарегистр. 06-12-2002 | Отправлено: 18:11 28-11-2023

MBK2

Silver Member

Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору

Цитата:

1, 634511, 17, 286, 8, 56789,
& 10001, 33888, 1000000, 4000000, 55000000,
& 1111, 99999933, 10, 14, 666/

Это не штатный рандом, а какая то колхозная фигня. Так и передайте

Всего записей: 4576 | Зарегистр. 18-09-2018 | Отправлено: 18:45 28-11-2023 | Исправлено: MBK2, 18:45 28-11-2023

akaGM

Platinum Member

Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору

передал...

это не штатный рандом, а мой тестовый пример использования, точно такую же последовательность я использовал в Си-версии, надо же с чем-то сравнивать...
и полное совпадение -- критериий того что я нигде не накосячил по дороге...
сложно понять?
конечно сложно, проще считать меня колхозным дебилом...

Код:

F C Delphi
0.070572 0.070572 0.070572
0.593556 0.593556 0.593556
0.966379 0.966379 0.966379
0.230200 0.230200 0.230200
0.506667 0.506667 0.506667
0.742368 0.742368 0.742368
0.736956 0.736956 0.736956
0.493384 0.493384 0.493384
0.141041 0.141041 0.141041
0.029140 0.029140 0.029140

я заодно и на дельфях написал -- пусть будет...

Всего записей: 24120 | Зарегистр. 06-12-2002 | Отправлено: 18:55 28-11-2023

MBK2

Silver Member

Редактировать | Профиль | Сообщение | Цитировать | Сообщить модератору

Цитата:

проще считать меня колхозным дебилом...

Нет, но я все время говорил, что вы интересная личность
1) Вы приводите свой код как финал-апофеоз титанической работы и пример кошерного использования
2) Я вам указываю на его неполноту и привожу действительно кошерный пример данного кода
3) Вы упираетесь и шлете бесполезный фрагмент
4) Я указываю на его бесполезность и шлю еще полезный фрагмент
5) Вы наконец шлете (как оказалось) слабый тестовый недостающий фрагмент кода
6) Я указываю на его слабину и шлю еще один пример хорошего годного кода
7) Вы совсем обижаетесь и в отместку зачем то шлете мне тестовую таблицу результатов своего кода
.......
!!!!!

Всего записей: 4576 | Зарегистр. 18-09-2018 | Отправлено: 19:12 28-11-2023

Открыть новую тему

Написать ответ в эту тему

Страницы: 1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60 61 62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78 79 80 81 82 83 84 85 86 87 88 89 90 91 92 93 94 95 96 97 98 99 100 101 102 103 104 105 106 107 108 109 110 111 112 113 114 115 116 117 118 119 120 121 122 123 124 125 126 127 128 129 130 131 132 133 134 135 136 137 138 139 140 141 142 143 144 145 146 147 148 149 150 151 152 153 154 155 156 157 158 159 160 161 162 163 164 165 166 167 168 169 170 171 172 173 174 175 176 177 178 179 180 181 182 183 184 185 186 187 188 189 190 191 192 193 194 195 196 197 198 199 200 201 202 203 204 205 206 207 208 209 210 211 212 213 214 215 216 217 218 219 220 221 222 223 224 225 226 227 228 229 230 231 232 233 234 235 236 237 238 239 240 241 242 243 244 245 246 247 248 249 250 251 252 253 254 255 256 257 258 259 260 261 262 263 264 265 266 267 268 269 270 271 272 273 274 275 276 277 278 279 280 281 282 283 284 285 286 287 288 289 290 291 292 293 294 295 296 297 298 299 300 301 302 303 304 305 306 307 308 309 310 311 312 313 314 315 316 317 318 319 320 321 322 323 324 325 326 327 328 329 330

Компьютерный форум Ru.Board » Компьютеры » Прикладное программирование » Вопросы программирования на FORTRAN (ФОРТРАН)

Реклама на форуме Ru.Board.

Powered by Ikonboard "v2.1.7b" © 2000 Ikonboard.com
Modified by Ru.B0ard
© Ru.B0ard 2000-2024

BitCoin: 1NGG1chHtUvrtEqjeerQCKDMUi6S6CG4iC