← All posts tagged fortran

PineappleZombie

LAPACK, что ты делаешь? Прекрати!

The following LAPACK subroutine base names are supported for single precision (s) and double precision (d+), in both the high-level and middle-level interfaces:

bdsdc bdsqr disna gbbrd gbcon gbequ gbequb gbrfs gbrfsx gbsv gbsvx gbsvxx gbtrf gbtrs gebak gebal gebrd gecon geequ geequb gees geesx geev geevx gehrd gejsv gelqf gels gelsd gelss gelsy geqlf geqp3 geqpf geqrf geqrfp gerfs gerfsx gerqf gesdd gesv gesvd gesvj gesvx gesvxx getrf getri getrs ggbak ggbal gges ggesx ggev ggevx ggglm gghrd gglse ggqrf ggrqf ggsvd ggsvp gtcon gtrfs gtsv gtsvx gttrf gttrs hgeqz hsein hseqr opgtr opmtr orgbr orghr orglq orgql orgqr orgrq orgtr ormbr ormhr ormlq ormql ormqr ormrq ormrz ormtr pbcon pbequ pbrfs pbstf pbsv pbsvx pbtrf pbtrs pftrf pftri pftrs pocon poequ poequb porfs porfsx posv posvx posvxx potrf potri potrs ppcon ppequ pprfs ppsv ppsvx pptrf pptri pptrs pstrf ptcon pteqr ptrfs ptsv ptsvx pttrf pttrs sbev sbevd sbevx sbgst sbgv sbgvd sbgvx sbtrd sfrk spcon spev spevd spevx spgst spgv spgvd spgvx sprfs spsv spsvx sptrd sptrf sptri sptrs stebz stedc stegr stein stemr steqr sterf stev stevd stevr stevx sycon syequb syev syevd syevr syevx sygst sygv sygvd sygvx syrfs syrfsx sysv sysvx sysvxx sytrd sytrf sytri sytrs tbcon tbrfs tbtrs tfsm tftri tfttp tfttr tgevc tgexc tgsen tgsja tgsna tgsyl tpcon tprfs tptri tptrs tpttf tpttr trcon trevc trexc trrfs trsen trsna trsyl trtri trtrs trttf trttp tzrzf

PineappleZombie

Коллега портирует систему обработки, написанную на фортране, на 64 бита. В процессе обнаружилось, что массив в 5000 элементов пишется 18000 элементов, ну и там ещё по мелочи. Спасибо gfortran за -fbounds-check

Особенно меня радует тот факт, что с помощью этой системы сделано, не скажу сразу сколько, но много публикаций

PineappleZombie

Нет, я хуею с этих физиков! И так уже 50-й коммит:

real axcal(7,7,9),aycal(7,7,9)
integer nudum
c Implicits
— REAL asava,aplv
+ REAL aplv
REAL ameanv,amav,asaha
REAL aplh,ameanh,wqq,awt,aaa,axmean,assa,amah,axma
REAL assah,assav,assath,assatv,axmat,flsah,flsav
@@ -635,7 +635,6 @@
imuxa=1
endif
if(imuxa.eq.1.and.assa.lt.250.0)then
— asava=assa
nmuver=nmuver+1
lmuver(nmuver)=is
endif