Blame

Date:
Message:
Why not?
Actions:
History | Blob | Raw File
  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 <stdio.h>

  4 28994509 2004-04-21 devnull #include "map.h"

  5 28994509 2004-04-21 devnull #include "iplot.h"

  6 28994509 2004-04-21 devnull 

  7 28994509 2004-04-21 devnull #define NSYMBOL 20

  8 28994509 2004-04-21 devnull 

  9 28994509 2004-04-21 devnull enum flag { POINT,ENDSEG,ENDSYM };

 10 28994509 2004-04-21 devnull struct symb {

 11 28994509 2004-04-21 devnull 	double x, y;

 12 28994509 2004-04-21 devnull 	char name[10+1];

 13 28994509 2004-04-21 devnull 	enum flag flag;

 14 28994509 2004-04-21 devnull } *symbol[NSYMBOL];

 15 28994509 2004-04-21 devnull 

 16 28994509 2004-04-21 devnull static int nsymbol;

 17 28994509 2004-04-21 devnull static double halfrange = 1;

 18 28994509 2004-04-21 devnull extern int halfwidth;

 19 28994509 2004-04-21 devnull extern int vflag;

 20 28994509 2004-04-21 devnull 

 21 28994509 2004-04-21 devnull static int	getrange(FILE *);

 22 28994509 2004-04-21 devnull static int	getsymbol(FILE *, int);

 23 28994509 2004-04-21 devnull static void	setrot(struct place *, double, int);

 24 28994509 2004-04-21 devnull static void	dorot(struct symb *, double *, double *);

 25 28994509 2004-04-21 devnull 

 26 28994509 2004-04-21 devnull 

 27 28994509 2004-04-21 devnull void

 28 28994509 2004-04-21 devnull getsyms(char *file)

 29 28994509 2004-04-21 devnull {

 30 28994509 2004-04-21 devnull 	FILE *sf = fopen(file,"r");

 31 28994509 2004-04-21 devnull 	if(sf==0)

 32 28994509 2004-04-21 devnull 		filerror("cannot open", file);

 33 28994509 2004-04-21 devnull 	while(nsymbol<NSYMBOL-1 && getsymbol(sf,nsymbol))

 34 28994509 2004-04-21 devnull 		nsymbol++;

 35 28994509 2004-04-21 devnull 	fclose(sf);

 36 28994509 2004-04-21 devnull }

 37 28994509 2004-04-21 devnull 

 38 28994509 2004-04-21 devnull static int

 39 28994509 2004-04-21 devnull getsymbol(FILE *sf, int n)

 40 28994509 2004-04-21 devnull {

 41 28994509 2004-04-21 devnull 	double x,y;

 42 28994509 2004-04-21 devnull 	char s[2];

 43 28994509 2004-04-21 devnull 	int i;

 44 28994509 2004-04-21 devnull 	struct symb *sp;

 45 28994509 2004-04-21 devnull 	for(;;) {

 46 28994509 2004-04-21 devnull 		if(fscanf(sf,"%1s",s)==EOF)

 47 28994509 2004-04-21 devnull 			return 0;

 48 28994509 2004-04-21 devnull 		switch(s[0]) {

 49 28994509 2004-04-21 devnull 		case ':':

 50 28994509 2004-04-21 devnull 			break;

 51 28994509 2004-04-21 devnull 		case 'o':

 52 28994509 2004-04-21 devnull 		case 'c':	/* cl */

 53 28994509 2004-04-21 devnull 			fscanf(sf,"%*[^\n]");

 54 28994509 2004-04-21 devnull 			continue;

 55 28994509 2004-04-21 devnull 		case 'r':

 56 28994509 2004-04-21 devnull 			if(getrange(sf))

 57 28994509 2004-04-21 devnull 				continue;

 58 28994509 2004-04-21 devnull 		default:

 59 28994509 2004-04-21 devnull 			error("-y file syntax error");

 60 28994509 2004-04-21 devnull 		}

 61 28994509 2004-04-21 devnull 		break;

 62 28994509 2004-04-21 devnull 	}

 63 28994509 2004-04-21 devnull 	sp = (struct symb*)malloc(sizeof(struct symb));

 64 28994509 2004-04-21 devnull 	symbol[n] = sp;

 65 28994509 2004-04-21 devnull 	if(fscanf(sf,"%10s",sp->name)!=1)

 66 28994509 2004-04-21 devnull 		return 0;

 67 28994509 2004-04-21 devnull 	i = 0;

 68 28994509 2004-04-21 devnull 	while(fscanf(sf,"%1s",s)!=EOF) {

 69 28994509 2004-04-21 devnull 		switch(s[0]) {

 70 28994509 2004-04-21 devnull 		case 'r':

 71 28994509 2004-04-21 devnull 			if(!getrange(sf))

 72 28994509 2004-04-21 devnull 				break;

 73 28994509 2004-04-21 devnull 			continue;

 74 28994509 2004-04-21 devnull 		case 'm':

 75 28994509 2004-04-21 devnull 			if(i>0)

 76 28994509 2004-04-21 devnull 				symbol[n][i-1].flag = ENDSEG;

 77 28994509 2004-04-21 devnull 			continue;

 78 28994509 2004-04-21 devnull 		case ':':

 79 28994509 2004-04-21 devnull 			ungetc(s[0],sf);

 80 28994509 2004-04-21 devnull 			break;

 81 28994509 2004-04-21 devnull 		default:

 82 28994509 2004-04-21 devnull 			ungetc(s[0],sf);

 83 28994509 2004-04-21 devnull 		case 'v':

 84 28994509 2004-04-21 devnull 			if(fscanf(sf,"%lf %lf",&x,&y)!=2)

 85 28994509 2004-04-21 devnull 				break;

 86 28994509 2004-04-21 devnull 			sp[i].x = x*halfwidth/halfrange;

 87 28994509 2004-04-21 devnull 			sp[i].y = y*halfwidth/halfrange;

 88 28994509 2004-04-21 devnull 			sp[i].flag = POINT;

 89 28994509 2004-04-21 devnull 			i++;

 90 28994509 2004-04-21 devnull 			sp = symbol[n] = (struct symb*)realloc(symbol[n],

 91 28994509 2004-04-21 devnull 					(i+1)*sizeof(struct symb));

 92 28994509 2004-04-21 devnull 			continue;

 93 28994509 2004-04-21 devnull 		}

 94 28994509 2004-04-21 devnull 		break;

 95 28994509 2004-04-21 devnull 	}

 96 28994509 2004-04-21 devnull 	if(i>0)

 97 28994509 2004-04-21 devnull 		symbol[n][i-1].flag = ENDSYM;

 98 28994509 2004-04-21 devnull 	else

 99 28994509 2004-04-21 devnull 		symbol[n] = 0;

100 28994509 2004-04-21 devnull 	return 1;

101 28994509 2004-04-21 devnull }

102 28994509 2004-04-21 devnull 

103 28994509 2004-04-21 devnull static int

104 28994509 2004-04-21 devnull getrange(FILE *sf)

105 28994509 2004-04-21 devnull {

106 28994509 2004-04-21 devnull 	double x,y,xmin,ymin;

107 28994509 2004-04-21 devnull 	if(fscanf(sf,"%*s %lf %lf %lf %lf",

108 28994509 2004-04-21 devnull 		&xmin,&ymin,&x,&y)!=4)

109 28994509 2004-04-21 devnull 		return 0;

110 28994509 2004-04-21 devnull 	x -= xmin;

111 28994509 2004-04-21 devnull 	y -= ymin;

112 28994509 2004-04-21 devnull 	halfrange = (x>y? x: y)/2;

113 28994509 2004-04-21 devnull 	if(halfrange<=0)

114 28994509 2004-04-21 devnull 		error("bad ra command in -y file");

115 28994509 2004-04-21 devnull 	return 1;

116 28994509 2004-04-21 devnull }

117 28994509 2004-04-21 devnull 

118 28994509 2004-04-21 devnull /* r=0 upright;=1 normal;=-1 reverse*/

119 28994509 2004-04-21 devnull int

120 28994509 2004-04-21 devnull putsym(struct place *p, char *name, double s, int r)

121 28994509 2004-04-21 devnull {

122 28994509 2004-04-21 devnull 	int x,y,n;

123 28994509 2004-04-21 devnull 	struct symb *sp;

124 28994509 2004-04-21 devnull 	double dx,dy;

125 28994509 2004-04-21 devnull 	int conn = 0;

126 28994509 2004-04-21 devnull 	for(n=0; symbol[n]; n++)

127 28994509 2004-04-21 devnull 		if(strcmp(name,symbol[n]->name)==0)

128 28994509 2004-04-21 devnull 			break;

129 28994509 2004-04-21 devnull 	sp = symbol[n];

130 28994509 2004-04-21 devnull 	if(sp==0)

131 28994509 2004-04-21 devnull 		return 0;

132 28994509 2004-04-21 devnull 	if(doproj(p,&x,&y)*vflag <= 0)

133 28994509 2004-04-21 devnull 		return 1;

134 28994509 2004-04-21 devnull 	setrot(p,s,r);

135 28994509 2004-04-21 devnull 	for(;;) {

136 28994509 2004-04-21 devnull 		dorot(sp,&dx,&dy);

137 28994509 2004-04-21 devnull 		conn = cpoint(x+(int)dx,y+(int)dy,conn);

138 28994509 2004-04-21 devnull 		switch(sp->flag) {

139 28994509 2004-04-21 devnull 		case ENDSEG:

140 28994509 2004-04-21 devnull 			conn = 0;

141 28994509 2004-04-21 devnull 		case POINT:

142 28994509 2004-04-21 devnull 			sp++;

143 28994509 2004-04-21 devnull 			continue;

144 28994509 2004-04-21 devnull 		case ENDSYM:

145 28994509 2004-04-21 devnull 			break;

146 28994509 2004-04-21 devnull 		}

147 28994509 2004-04-21 devnull 		break;

148 28994509 2004-04-21 devnull 	}

149 28994509 2004-04-21 devnull 	return 1;

150 28994509 2004-04-21 devnull }

151 28994509 2004-04-21 devnull 

152 28994509 2004-04-21 devnull static double rot[2][2];

153 28994509 2004-04-21 devnull 

154 28994509 2004-04-21 devnull static void

155 28994509 2004-04-21 devnull setrot(struct place *p, double s, int r)

156 28994509 2004-04-21 devnull {

157 28994509 2004-04-21 devnull 	double x0,y0,x1,y1;

158 28994509 2004-04-21 devnull 	struct place up;

159 28994509 2004-04-21 devnull 	up = *p;

160 28994509 2004-04-21 devnull 	up.nlat.l += .5*RAD;

161 28994509 2004-04-21 devnull 	sincos(&up.nlat);

162 28994509 2004-04-21 devnull 	if(r&&(*projection)(p,&x0,&y0)) {

163 28994509 2004-04-21 devnull 		if((*projection)(&up,&x1,&y1)<=0) {

164 28994509 2004-04-21 devnull 			up.nlat.l -= RAD;

165 28994509 2004-04-21 devnull 			sincos(&up.nlat);

166 28994509 2004-04-21 devnull 			if((*projection)(&up,&x1,&y1)<=0)

167 28994509 2004-04-21 devnull 				goto unit;

168 28994509 2004-04-21 devnull 			x1 = x0 - x1;

169 28994509 2004-04-21 devnull 			y1 = y0 - y1;

170 28994509 2004-04-21 devnull 		} else {

171 28994509 2004-04-21 devnull 			x1 -= x0;

172 28994509 2004-04-21 devnull 			y1 -= y0;

173 28994509 2004-04-21 devnull 		}

174 28994509 2004-04-21 devnull 		x1 = r*x1;

175 28994509 2004-04-21 devnull 		s /= hypot(x1,y1);

176 28994509 2004-04-21 devnull 		rot[0][0] = y1*s;

177 28994509 2004-04-21 devnull 		rot[0][1] = x1*s;

178 28994509 2004-04-21 devnull 		rot[1][0] = -x1*s;

179 28994509 2004-04-21 devnull 		rot[1][1] = y1*s;

180 28994509 2004-04-21 devnull 	} else {

181 28994509 2004-04-21 devnull unit:

182 28994509 2004-04-21 devnull 		rot[0][0] = rot[1][1] = s;

183 28994509 2004-04-21 devnull 		rot[0][1] = rot[1][0] = 0;

184 28994509 2004-04-21 devnull 	}

185 28994509 2004-04-21 devnull }

186 28994509 2004-04-21 devnull 

187 28994509 2004-04-21 devnull static void

188 28994509 2004-04-21 devnull dorot(struct symb *sp, double *px, double *py)

189 28994509 2004-04-21 devnull {

190 28994509 2004-04-21 devnull 	*px = rot[0][0]*sp->x + rot[0][1]*sp->y;

191 28994509 2004-04-21 devnull 	*py = rot[1][0]*sp->x + rot[1][1]*sp->y;

192 28994509 2004-04-21 devnull }