Cải thiện độ chính xác trên các biến được xác định bởi số nguyên thương

Question

Nói rằng tôi có chương trình sau đây:

program derp 
    implicit none 
    integer, parameter :: ikind = selected_real_kind(18) 
    real (kind = ikind) :: a = 2.0/3.0 
    print*, a 
end program derp 

Chương trình derp đầu ra 0.6666666865348815917, mà rõ ràng là không 18 chữ số chính xác. Tuy nhiên, nếu tôi xác định a=2.0 và b=3.0 sử dụng cùng một phương pháp và thì xác định c=a/b Tôi nhận được kết quả là 0.666666666666666666685, điều này là tốt. Làm thế nào để tôi chỉ xác định một biến như là một thương số nguyên và có nó lưu trữ tất cả các chữ số chính xác tôi muốn từ selected_real_kind?

Answer 1

Hãy thử: real (kind = ikind) :: a = 2.0_ikind/3.0_ikind

Lý do là trong khi LHS là độ chính xác cao, RHS trong ví dụ mã của bạn, 2.0/3.0, không phải là. Fortran thực hiện phép tính đó với độ chính xác đơn và sau đó gán kết quả cho LHS. Phía RHS không được tính chính xác cao hơn vì LHS có độ chính xác cao. digits_kind là cách chỉ định loại hằng số digits.