Bug #469
Eckert IV runs forever in Linux 32 bits
| Status: | Closed | Start date: | 2013-12-24 | |
|---|---|---|---|---|
| Priority: | Normal | Due 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
Address issue #469
History
#1
Updated by Paul almost 5 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 almost 5 years ago
psbasemap -R-180/180/-90/90 -JKf0/18 -Ba0g0f0 -Bwesn+t"Eckert IV" > eckiv.ps
#3
Updated by Paul almost 5 years ago
and it is now working correctly for you?
#4
Updated by Eduardo almost 5 years ago
Yes, it works for me, and I get the same map in 32 and 64 bits. Thanks!