1 28994509 2004-04-21 devnull #include <u.h>
 2 28994509 2004-04-21 devnull #include <libc.h>
 3 28994509 2004-04-21 devnull #include "map.h"
 4 28994509 2004-04-21 devnull 
 5 28994509 2004-04-21 devnull static int
 6 28994509 2004-04-21 devnull Xmollweide(struct place *place, double *x, double *y)
 7 28994509 2004-04-21 devnull {
 8 28994509 2004-04-21 devnull 	double z;
 9 28994509 2004-04-21 devnull 	double w;
10 28994509 2004-04-21 devnull 	z = place->nlat.l;
11 28994509 2004-04-21 devnull 	if(fabs(z)<89.9*RAD)
12 28994509 2004-04-21 devnull 		do {	/*newton for 2z+sin2z=pi*sin(lat)*/
13 28994509 2004-04-21 devnull 			w = (2*z+sin(2*z)-PI*place->nlat.s)/(2+2*cos(2*z));
14 28994509 2004-04-21 devnull 			z -= w;
15 28994509 2004-04-21 devnull 		} while(fabs(w)>=.00001);
16 28994509 2004-04-21 devnull 	*y = sin(z);
17 28994509 2004-04-21 devnull 	*x = - (2/PI)*cos(z)*place->wlon.l;
18 28994509 2004-04-21 devnull 	return(1);
19 28994509 2004-04-21 devnull }
20 28994509 2004-04-21 devnull 
21 28994509 2004-04-21 devnull proj
22 28994509 2004-04-21 devnull mollweide(void)
23 28994509 2004-04-21 devnull {
24 28994509 2004-04-21 devnull 	return(Xmollweide);
25 28994509 2004-04-21 devnull }