คำถามในหัวข้อ 'fortran'
โปรแกรมความเที่ยงตรงสองเท่าของ Fortran พร้อมรูทีน MKL BLAS อย่างง่าย
ในการพยายามผสมความแม่นยำในโปรแกรมง่ายๆ โดยใช้ทั้งของจริงและสองเท่า และใช้รูทีน ddot จาก BLAS ฉันพบผลลัพธ์ที่ไม่ถูกต้องสำหรับชิ้นส่วนที่มีความแม่นยำสองเท่า นี่คือรหัส:
program test
!! adding this statement narrowed the issue down to ddot being...
3115 มุมมอง
schedule
02.11.2023
การแบ่งส่วนความผิด C และ fortran
------ main.c ---------
#include <stdio.h>
#include <stdlib.h>
#include <dlfcn.h>
#include <string.h>
int main()
{
char* lib_name = "./a.out";
int array[5] = {1,2,3,4,5};
int size_a =...
758 มุมมอง
schedule
02.01.2024
การอ่านไฟล์ข้อมูลภายนอกใน gnuplot
นี่คือความต่อเนื่องของคำถามก่อนหน้าของฉัน อย่างที่ฉันบอกไปว่าไฟล์ข้อมูลถูกสร้างขึ้นจากโค้ด Fortran ไฟล์ข้อมูลทั้งหมดประกอบด้วยข้อมูลสองคอลัมน์ ในโค้ด Fortran ฉันใช้ FORMAT (2(1X,D22.16)) ดังนั้นไฟล์ที่ส่งออกจะมีลักษณะเช่นนี้...
3587 มุมมอง
schedule
24.11.2023
รูปแบบที่ดีสำหรับการประกาศตัวแปรใน Fortran คืออะไร?
ฉันมีตัวแปรมากมายที่ฉันต้องประกาศ และโค้ดต้นฉบับมีลักษณะดังนี้:
DIMENSION energy_t(20000),nrt(20000),npsh(1000),xx(1000),yy(1000),
:step(1000),stepz(1000),r1(1000),rr(1000),ic(1000),diffrr(1000)
และฉันก็เขียนมันใหม่ดังนี้:...
334 มุมมอง
schedule
04.01.2024
ปรับโค้ดตัวอย่าง dgemm เพื่อใช้ sgemm (scalapack)
ฉันต้องทำโปรแกรมต่อไปนี้ (จาก http://www.netlib.org/scalapack/examples/pblas.tgz ) ทำงานร่วมกับ SGEMM ฉันต้องเปลี่ยนแปลงอะไรบ้างเพื่อให้มันใช้งานได้? ความรู้ของฉันเกี่ยวกับ Fortran ค่อนข้างจำกัด...
483 มุมมอง
schedule
30.11.2023
เครื่องกำเนิดสมการเชิงเส้น - เอาต์พุตเป็น 0 ทั้งหมดหรือไม่
ฉันพยายามสร้างเครื่องกำเนิดตัวเลขสุ่มเทียม LCG ขั้นพื้นฐานใน Fortran 77 เพื่อพิมพ์ตัวเลขสุ่ม 1,000 ตัวลงในไฟล์ แต่ด้วยเหตุผลใดก็ตามเอาต์พุตจึงเป็นเพียง 1,000 0s โค้ดทั้งหมดค่อนข้างสั้น ดังนั้นฉันจึงหวีมันหลายครั้งและพยายามเปลี่ยนแปลงบางสิ่ง...
139 มุมมอง
schedule
05.01.2024
จะสร้างและใช้ส่วนขยายประเภทอาร์เรย์ใน Fortran ได้อย่างไร [ทำซ้ำ]
ฉันมีคลาสฐานนามธรรม ฉันสร้างส่วนขยายของคลาสพื้นฐานสองอัน ฉันต้องการจัดเก็บอ็อบเจ็กต์ทั้งหมดที่สร้างจากคลาสขยายอันใดอันหนึ่งไว้ในอาร์เรย์เดียว ฉันเชื่อว่าฉันสามารถทำเช่นนี้ได้โดยการสร้างอาร์เรย์ของพอยน์เตอร์ที่ชี้ไปที่คลาสพื้นฐาน...
57 มุมมอง
schedule
04.01.2024
การพิมพ์อาร์เรย์ที่จัดสรรได้ใน Fortran ด้วย gdb: opcode นิพจน์แคระที่ไม่สามารถจัดการได้ 0x97 [ซ้ำกัน]
ฉันได้แก้ไขโค้ดต่อไปนี้ใน Cygwin และ Eclipse โดยใช้ gdb เป็นตัวดีบักเกอร์:
program codetest
implicit none
integer, parameter :: dp = kind(1.0d0)
integer, parameter :: N = 10
real(dp), dimension(:), allocatable :: vector...
387 มุมมอง
schedule
31.10.2023
คุณสมบัติเชิงวัตถุด้วยคอมไพเลอร์ Solaris/Oracle Fortran
ฉันพยายามรวบรวมโค้ด O-O fortran ด้วยคอมไพเลอร์ fortran ของชุด Oracle Solaris Studio 12.4 (ล่าสุดเท่าที่ฉันรู้) แต่คอมไพเลอร์ขัดข้อง นี่คือปัญหาของฉันในเวอร์ชันที่เรียบง่าย
ฉันกำหนดประเภทง่าย ๆ สองประเภทโดยแต่ละประเภทมีขั้นตอนผูกมัดประเภทเดียว...
119 มุมมอง
schedule
16.11.2023
ปัญหา makefile ของ gfortran
ฉันมีโปรแกรมเขียนบน Fortran ซึ่งใช้ไลบรารี่ สำหรับการคอมไพล์ ฉันใช้เชลล์สคริปต์อย่างง่าย และทุกอย่างเรียบร้อยดี:
gfortran -c mod_input.f90 mod_cpp.f90 mod_data.f90 main.f90 b_coils.f
gfortran -o หลัก mod_input.o mod_cpp.o mod_data.o...
303 มุมมอง
schedule
26.11.2023
จะแก้ไขข้อผิดพลาดทางไวยากรณ์ใน Data Statement ได้อย่างไร
สวัสดี ฉันมาใหม่และต้องการแก้ปัญหานี้:
do k=1,31
Data H(1,k)/0/
End do
do l=1,21
Data H(l,1)/0.5*(l-1)/
End do
do m=31,41
Data H(17,m)/0/
End do
do n=17,21
Data H(n,41)/0.5*(n-17)/
End do
ฉันได้รับข้อผิดพลาดสำหรับ l และ n...
589 มุมมอง
schedule
08.12.2023
ข้อผิดพลาดในการแบ่งส่วนโดยใช้ f2py
นี่คือการติดตามคำถามที่ฉันถามที่นี่... ฉันจะจัดสรรอาร์เรย์อินพุตด้วย f2py ได้อย่างไร
ฉันได้ปรับปรุงคำถามเพื่อตอบสนองต่อความคิดเห็น คำถามเดิมไม่ชัดเจนและยาวเกินไป ดูความคิดเห็นที่ 1 ด้านล่าง
ฉันกำลังพยายามใช้ f2py เพื่อรวมรูทีนย่อย/ฟังก์ชัน...
533 มุมมอง
schedule
10.01.2024
คอมไพล์ด้วย cygwin/gfortran หรือ ifort ด้วย acml
ฉันต้องการรวบรวมโปรแกรมของฉันที่เชื่อมโยงกับ acml ใน cygwin (win64) โดยมี gfortran เชื่อมโยงกับ acml:
$ gfortran empty.f90 -L/c:/cygwin64/acml5.3.1/win64/lib/libacml_dll.lib
อัตราผลตอบแทน
/tmp/ccKechZN.o:empty.f90:(.text+0xff): undefined...
598 มุมมอง
schedule
18.12.2023
การสลายตัวของ 2D FFT (และ inverse fft) ด้วย fftw ใน fortran77
[แก้ไข 1] เพิ่มตัวเลขเพื่อแสดงข้อมูลต้นฉบับและข้อมูลที่ได้รับ
[แก้ไข 2] ฉันพบข้อผิดพลาด ฉันใช้ fftw_measure แทน fftw_estimate ในการเรียก dfftw_plan_many_dft
[แก้ไข 3] แก้ไขการพิมพ์ผิดในรหัส (แทนที่ u ด้วย u2d ใน dfftw_execute_dft_r2c )...
551 มุมมอง
schedule
02.01.2024
ข้อผิดพลาด: ต้องใช้เครื่องหมายจุลภาคหลัง P descriptor ในรูปแบบสตริง
เมื่อฉันพยายามคอมไพล์ calrank.for (CALRANK_v7.0.0_L140912.zip) ด้วย gfortran ฉันได้รับข้อผิดพลาดนี้:
> gfortran -m64 -c -o calrank.o calrank.for
calrank.for:1922:32:
write(io,'(2f10.3,1p1000(2x,e16.7,1x,2a10,1x))')...
295 มุมมอง
schedule
27.12.2023
Fortran 90 สร้างคำตอบที่ไม่ถูกต้อง สร้าง 0 เมื่อคำตอบควรเป็น 1 ฉันกำลังมองข้ามสิ่งใดไปหรือเปล่า? [ทำซ้ำ]
นี่คือรหัส Fortran 90 ที่ฉันเขียน เป็นฟังก์ชันที่ฉันใช้ในโค้ดที่ใหญ่กว่า แต่ฉันได้รับผลลัพธ์ที่ไม่ถูกต้อง
ฉันรู้ว่าในทางคณิตศาสตร์คำตอบควรเป็น 1 (หรือใกล้เคียงกับ 1 เนื่องจากการคำนวณ) แต่ผลลัพธ์คำตอบของฉันคือ 0
ฉันทำอะไรผิดหรือเปล่า?...
72 มุมมอง
schedule
25.11.2023
องค์ประกอบอาเรย์แบบคงที่ใน Fortran สมัยใหม่ [ซ้ำกัน]
เพื่อความสะดวก ฉันต้องการที่จะอ้างอิงองค์ประกอบของตารางการค้นหาแบบคงที่ที่ใช้งานเป็นอาร์เรย์ 1d เช่นนี้
Item = [1,7,2,5,3](I)
โดยที่ I เป็นดัชนีจำนวนเต็มของอาร์เรย์แบบคงที่ในวงเล็บเหลี่ยม...
35 มุมมอง
schedule
05.11.2023
งานล้มเหลวขณะใช้ srun หรือ mpirun ใน slurm
ฉันกำลังพยายามส่งงานกับสเลม อย่างไรก็ตาม งานจะล้มเหลวหากฉันใช้ srun หรือ mpirun อย่างไรก็ตาม มันทำงานได้ดีด้วย mpiexec แม้ว่าจะทำงานด้วยกระบวนการเดียวเท่านั้น แม้ว่าจะมีหลายโหนดและหลายคอร์ก็ตามที่ได้รับการจัดสรร
คำสั่งจริงที่ใช้คือ:
srun...
770 มุมมอง
schedule
26.12.2023
