Bug in PJ_cea.c

[proj4-bot]

Reported by tharan on 8 Feb 2011 21:56 UTC
Checked out proj4 development version 1962 on 2011/02/04. Running make using cc --version cc (SUSE Linux) 4.3.2 revision 141291 on uname -a Linux snow 2.6.27.54-0.2-default #1 SMP 2010-10-19 18:40:07 +0200 x86_64 x86_64 x86_64 GNU/Linux.

PJ_cea.c - In function ENTRY1(cea, apa) in the first "if" statement, removed "else t = 0.;" and replaced it with ";". Setting "t" to 0 when the if statement is false causes an incorrect value for "t = sin(t)" and the subsequent "P->k /= sqrt(1. - P->es * t * t);"

Refer to http://lists.maptools.org/pipermail/proj/2004-January/001051.html for a description of the bug that this new version of PJ_cea.c apparently fixes.

Migrated-From: https://trac.osgeo.org/proj/ticket/94

[proj4-bot]

Comment by warmerdam on 8 Feb 2011 23:36 UTC
Patch applied rather blindly in trunk (565f820f (SVN rev1967)). I am too lazy and dumb to establish if this is really justified so I'll just trust you on it.

[proj4-bot]

Attachment added by tharan on 8 Feb 2011 21:56 UTC
https://trac.osgeo.org/proj/attachment/ticket/94/PJ_cea.c