Unproper functionality of usuall things in my program on Intel Visual Fortran 10.1.019 64 bit

Les Neilson
Total Points:
5,816
Status Points:
5,316
Brown Belt
March 25, 2008 4:08 AM PDT
Rate
 
#1

Comparing floating-point numbers for equality is your problem.

You should really use a tolerance and check the abs of the difference between the two numbers.

eg :

real(8) :: diff
real(8) :: tolerance      ! set to some appropriate acceptable value

diff = omega - (omega_center-4.0d0*gamma)

if ( abs(diff) <= tolerance ) then
       ! the two numbers are considered to be equal

Les

Intel Software Network Forums Statistics

8458 users have contributed to 31571 threads and 100533 posts to date.
In the past 24 hours, we have 17 new thread(s) 131 new posts(s), and 152 new user(s).

In the past 3 days, the most popular thread for everyone has been gemm(A,A,A) like possible? The most posts were made to gemm(A,A,A) like possible? The post with the most views is Quoting - rase if (k.eq.0

Please welcome our newest member soundmyth