| Age | Commit message (Collapse) | Author |
|---|
|
Saves 2 sincos() and 1 atan2() calls.
With the following bench
```
#include "proj.h"
#include <stdlib.h>
#include <stdio.h>
#include <string.h>
int main(int argc, char* argv[])
{
if( argc != 2 ) {
fprintf(stderr, "Usage: bench_inv_cart fwd/inv\n");
exit(1);
}
PJ* p = proj_create(0, "+proj=cart");
const int dir = strcmp(argv[1], "inv") == 0 ? PJ_INV : PJ_FWD;
PJ_COORD coord;
if( dir == PJ_FWD )
{
coord.xyz.x = 3.14159/2;
coord.xyz.y = 3.14159/2;
coord.xyz.z = 100;
}
else
{
coord.xyz.x = 3e6;
coord.xyz.y = 3e6;
coord.xyz.z = 3e6;
}
for(int i = 0; i < 10* 1024 * 1024; i++ )
{
proj_trans(p, dir, coord);
}
proj_destroy(p);
return 0;
}
```
On Intel(R) Core(TM) i7-6700HQ CPU @ 2.60GHz
Time before: 2.37s
Time after: 1.57s
|
|
This should avoid issues with numerical stability as uncovered in
https://github.com/OSGeo/PROJ/issues/1906.
Practically speaking this change isn't going to affect real life
scenarios since the position of the center of the Earth is rarely
expressed in geodetic coordinates.
|
|
needs <limits.h>. Update scripts/reference_exported_symbols.txt and
src/proj_symbol_rename.h.
|
|
In that case, for a non-spherical ellipsoid, a phi = 180deg was returned,
which caused a division by zero in the foward path of moll.cpp
Fixup the latitude to be 0 when that happens.
Fixes https://bugs.chromium.org/p/oss-fuzz/issues/detail?id=14348
Credit to OSS Fuzz
|
|
Introduced by "Merge projects.h into proj_internal.h" 8ab6f683.
|
|
This communicates better what the unit actually is. It also opens the
door for potentially adding PJ_IO_UNITS_DEGREES in the future in order
to coordinates in degrees correctly as well.
|
|
|
|
structures
|
|
|
