Problem with C functions ?

Hi there,

I compiled a C-funtion to compute the azimuth of a LSEG but the function
never
returns the correct value, but elog( ... ) report the correct value.
So there seems to be a problem returning doubles or floats ??

compile the folling function with
gcc -fPIC -O2 -shared  -o lseg_azimuth.o -c lseg_azimuth.c
-I/opt/local/postgres/include -L/opt/local/postgres/lib
ld -G -Bdynamic -o libgeom.so   lseg_azimuth.o

#include <math.h>

#include "postgres.h"
#include "utils/geo_decls.h"

double lseg_azimuth( LSEG* lseg );
float add_one( float arg );

double lseg_azimuth( LSEG* lseg )
{
  double x1, y1;
  double x2, y2;
  double ori;

  x1 = lseg->p[0].x;
  y1 = lseg->p[0].y;

  x2 = lseg->p[1].x;
  y2 = lseg->p[1].y;

  ori = atan2( (x2-x1), (y2-y1));
  ori = (ori / M_PI) * 200.0;

  if ( ori < 0.0 )
    ori += 400.0;

  elog( NOTICE, "ori = %8.3lf", ori );

  return( ori );
}


then use:

create function lseg_azimuth(lseg) returns float8
  as ' ... path to the shared lib you just created ...' language 'c';


create table xxx (geometry lseg);

insert into xxx (geometry) values
('[(4148.675,5414.255),(4151.405,5414.255)]');

select lseg_azimuth(geometry) from xxx;

will return:

NOTICE:  ori =  100.000
lseg_azimuth
------------
    4148.675
(1 row)

But that's obviously not what it should be ...

Maybe you have any suggestions ?

many thanks
Manuel
--
Manuel Weindorf (weindorf@ipf.bau-verm.uni-karlsruhe.de)
Institut fuer Photogrammetrie und Fernerkundung, Universitaet Karlsruhe
Postfach 6980 D-76128 Karlsruhe Tel. +49721 6086010 Fax +49721 694568
http://www-ipf.bau-verm.uni-karlsruhe.de