/*  * This file is part of the UCB release of Plan 9. It is subject to the license * terms in the LICENSE file found in the top-level directory of this * distribution and at http://akaros.cs.berkeley.edu/files/Plan9License. No * part of the UCB release of Plan 9, including this file, may be copied, * modified, propagated, or distributed except according to the terms contained * in the LICENSE file. */ #include <u.h>#include <libc.h>#include "map.h" static intXlagrange(struct place *place, double *x, double *y){	double z1,z2;	double w1,w2,t1,t2;	struct place p;	copyplace(place,&p);	if(place->nlat.l<0) {		p.nlat.l = -p.nlat.l;		p.nlat.s = -p.nlat.s;	}	Xstereographic(&p,&z1,&z2);	csqrt(-z2/2,z1/2,&w1,&w2);	cdiv(w1-1,w2,w1+1,w2,&t1,&t2);	*y = -t1;	*x = t2;	if(place->nlat.l<0)		*y = -*y;	return(1);} projlagrange(void){	return(Xlagrange);}