Bug #469

Eckert IV runs forever in Linux 32 bits

Added by Eduardo over 3 years ago. Updated over 3 years ago.

Status:ClosedStart date:2013-12-24
Priority:NormalDue date:
Assignee:-% Done:

0%

Category:-
Target version:Candidate for next minor release
Affected version:5.x-svn Platform:Linux

Description

For example, psbasemap is trapped forever in the loop (GMT_eckert4):

do {
sincos (phi, &s, &c);
delta = -(phi + s * c + 2.0 * s - (2.0 + M_PI_2) * s_lat) / (2.0 * c * (1.0 + c));
phi += delta;
}
while (fabs(delta) > GMT_PROJ_CONV_LIMIT && n_iter < GMT_PROJ_MAX_ITERATIONS);

Perphaps because GMT_PROJ_CONV_LIMIT is small for 32 bits or because n_iter never gets incremented.

Associated revisions

Revision 12682
Added by Paul over 3 years ago

Address issue #469

History

#1 Updated by Paul over 3 years ago

  • Status changed from New to Resolved

Fixed in GMT 4 (r10172) and GMT5 (12682). Please share the psbasemap command that exposed this problem.

#2 Updated by Eduardo over 3 years ago

psbasemap -R-180/180/-90/90 -JKf0/18 -Ba0g0f0 -Bwesn+t"Eckert IV" > eckiv.ps

#3 Updated by Paul over 3 years ago

and it is now working correctly for you?

#4 Updated by Eduardo over 3 years ago

Yes, it works for me, and I get the same map in 32 and 64 bits. Thanks!

#5 Updated by Paul over 3 years ago

  • Status changed from Resolved to Closed

Closing this issue

Also available in: Atom PDF