regress.c
上传用户:blenddy
上传日期:2007-01-07
资源大小:6495k
文件大小:14k
源码类别:

数据库系统

开发平台:

Unix_Linux

  1. /*
  2.  * $Header: /usr/local/cvsroot/pgsql/src/test/regress/regress.c,v 1.32 1999/05/25 16:15:26 momjian Exp $
  3.  */
  4. #include <float.h> /* faked on sunos */
  5. #include <stdio.h>
  6. #include <string.h> /* for MemSet() */
  7. #include <postgres.h>
  8. #include "utils/geo_decls.h" /* includes <math.h> */
  9. #include "executor/executor.h" /* For GetAttributeByName */
  10. #define P_MAXDIG 12
  11. #define LDELIM '('
  12. #define RDELIM ')'
  13. #define DELIM ','
  14. typedef void *TUPLE;
  15. extern double *regress_dist_ptpath(Point *pt, PATH *path);
  16. extern double *regress_path_dist(PATH *p1, PATH *p2);
  17. extern PATH *poly2path(POLYGON *poly);
  18. extern Point *interpt_pp(PATH *p1, PATH *p2);
  19. extern void regress_lseg_construct(LSEG *lseg, Point *pt1, Point *pt2);
  20. extern char overpaid(TUPLE tuple);
  21. extern int boxarea(BOX *box);
  22. extern char *reverse_name(char *string);
  23. /*
  24. ** Distance from a point to a path
  25. */
  26. double *
  27. regress_dist_ptpath(pt, path)
  28. Point    *pt;
  29. PATH    *path;
  30. {
  31. double    *result;
  32. double    *tmp;
  33. int i;
  34. LSEG lseg;
  35. switch (path->npts)
  36. {
  37. case 0:
  38. result = palloc(sizeof(double));
  39. *result = Abs((double) DBL_MAX); /* +infinity */
  40. break;
  41. case 1:
  42. result = point_distance(pt, &path->p[0]);
  43. break;
  44. default:
  45. /*
  46.  * the distance from a point to a path is the smallest
  47.  * distance from the point to any of its constituent segments.
  48.  */
  49. Assert(path->npts > 1);
  50. result = palloc(sizeof(double));
  51. for (i = 0; i < path->npts - 1; ++i)
  52. {
  53. regress_lseg_construct(&lseg, &path->p[i], &path->p[i + 1]);
  54. tmp = dist_ps(pt, &lseg);
  55. if (i == 0 || *tmp < *result)
  56. *result = *tmp;
  57. pfree(tmp);
  58. }
  59. break;
  60. }
  61. return result;
  62. }
  63. /* this essentially does a cartesian product of the lsegs in the
  64.    two paths, and finds the min distance between any two lsegs */
  65. double *
  66. regress_path_dist(p1, p2)
  67. PATH    *p1;
  68. PATH    *p2;
  69. {
  70. double    *min,
  71.    *tmp;
  72. int i,
  73. j;
  74. LSEG seg1,
  75. seg2;
  76. regress_lseg_construct(&seg1, &p1->p[0], &p1->p[1]);
  77. regress_lseg_construct(&seg2, &p2->p[0], &p2->p[1]);
  78. min = lseg_distance(&seg1, &seg2);
  79. for (i = 0; i < p1->npts - 1; i++)
  80. for (j = 0; j < p2->npts - 1; j++)
  81. {
  82. regress_lseg_construct(&seg1, &p1->p[i], &p1->p[i + 1]);
  83. regress_lseg_construct(&seg2, &p2->p[j], &p2->p[j + 1]);
  84. if (*min < *(tmp = lseg_distance(&seg1, &seg2)))
  85. *min = *tmp;
  86. pfree(tmp);
  87. }
  88. return min;
  89. }
  90. PATH *
  91. poly2path(poly)
  92. POLYGON    *poly;
  93. {
  94. int i;
  95. char    *output = (char *) palloc(2 * (P_MAXDIG + 1) * poly->npts + 64);
  96. char buf[2 * (P_MAXDIG) + 20];
  97. sprintf(output, "(1, %*d", P_MAXDIG, poly->npts);
  98. for (i = 0; i < poly->npts; i++)
  99. {
  100. sprintf(buf, ",%*g,%*g", P_MAXDIG, poly->p[i].x, P_MAXDIG, poly->p[i].y);
  101. strcat(output, buf);
  102. }
  103. sprintf(buf, "%c", RDELIM);
  104. strcat(output, buf);
  105. return path_in(output);
  106. }
  107. /* return the point where two paths intersect. Assumes that they do. */
  108. Point *
  109. interpt_pp(p1, p2)
  110. PATH    *p1;
  111. PATH    *p2;
  112. {
  113. Point    *retval;
  114. int i,
  115. j;
  116. LSEG seg1,
  117. seg2;
  118. #ifdef NOT_USED
  119. LINE    *ln;
  120. #endif
  121. bool found; /* We've found the intersection */
  122. found = false; /* Haven't found it yet */
  123. for (i = 0; i < p1->npts - 1 && !found; i++)
  124. for (j = 0; j < p2->npts - 1 && !found; j++)
  125. {
  126. regress_lseg_construct(&seg1, &p1->p[i], &p1->p[i + 1]);
  127. regress_lseg_construct(&seg2, &p2->p[j], &p2->p[j + 1]);
  128. if (lseg_intersect(&seg1, &seg2))
  129. found = true;
  130. }
  131. #ifdef NOT_USED
  132. ln = line_construct_pp(&seg2.p[0], &seg2.p[1]);
  133. retval = interpt_sl(&seg1, ln);
  134. #endif
  135. retval = lseg_interpt(&seg1, &seg2);
  136. return retval;
  137. }
  138. /* like lseg_construct, but assume space already allocated */
  139. void
  140. regress_lseg_construct(lseg, pt1, pt2)
  141. LSEG    *lseg;
  142. Point    *pt1;
  143. Point    *pt2;
  144. {
  145. lseg->p[0].x = pt1->x;
  146. lseg->p[0].y = pt1->y;
  147. lseg->p[1].x = pt2->x;
  148. lseg->p[1].y = pt2->y;
  149. lseg->m = point_sl(pt1, pt2);
  150. }
  151. char
  152. overpaid(tuple)
  153. TUPLE tuple;
  154. {
  155. bool isnull;
  156. long salary;
  157. salary = (long) GetAttributeByName(tuple, "salary", &isnull);
  158. return salary > 699;
  159. }
  160. /* New type "widget"
  161.  * This used to be "circle", but I added circle to builtins,
  162.  * so needed to make sure the names do not collide. - tgl 97/04/21
  163.  */
  164. typedef struct
  165. {
  166. Point center;
  167. double radius;
  168. } WIDGET;
  169. WIDGET    *widget_in(char *str);
  170. char    *widget_out(WIDGET * widget);
  171. int pt_in_widget(Point *point, WIDGET * widget);
  172. #define NARGS 3
  173. WIDGET *
  174. widget_in(str)
  175. char    *str;
  176. {
  177. char    *p,
  178.    *coord[NARGS],
  179. buf2[1000];
  180. int i;
  181. WIDGET    *result;
  182. if (str == NULL)
  183. return NULL;
  184. for (i = 0, p = str; *p && i < NARGS && *p != RDELIM; p++)
  185. if (*p == ',' || (*p == LDELIM && !i))
  186. coord[i++] = p + 1;
  187. if (i < NARGS - 1)
  188. return NULL;
  189. result = (WIDGET *) palloc(sizeof(WIDGET));
  190. result->center.x = atof(coord[0]);
  191. result->center.y = atof(coord[1]);
  192. result->radius = atof(coord[2]);
  193. sprintf(buf2, "widget_in: read (%f, %f, %f)n", result->center.x,
  194. result->center.y, result->radius);
  195. return result;
  196. }
  197. char *
  198. widget_out(widget)
  199. WIDGET    *widget;
  200. {
  201. char    *result;
  202. if (widget == NULL)
  203. return NULL;
  204. result = (char *) palloc(60);
  205. sprintf(result, "(%g,%g,%g)",
  206. widget->center.x, widget->center.y, widget->radius);
  207. return result;
  208. }
  209. int
  210. pt_in_widget(point, widget)
  211. Point    *point;
  212. WIDGET    *widget;
  213. {
  214. extern double point_dt();
  215. return point_dt(point, &widget->center) < widget->radius;
  216. }
  217. #define ABS(X) ((X) > 0 ? (X) : -(X))
  218. int
  219. boxarea(box)
  220. BOX    *box;
  221. {
  222. int width,
  223. height;
  224. width = ABS(box->high.x - box->low.x);
  225. height = ABS(box->high.y - box->low.y);
  226. return width * height;
  227. }
  228. char *
  229. reverse_name(string)
  230. char    *string;
  231. {
  232. int i;
  233. int len;
  234. char    *new_string;
  235. if (!(new_string = palloc(NAMEDATALEN)))
  236. {
  237. fprintf(stderr, "reverse_name: palloc failedn");
  238. return NULL;
  239. }
  240. MemSet(new_string, 0, NAMEDATALEN);
  241. for (i = 0; i < NAMEDATALEN && string[i]; ++i)
  242. ;
  243. if (i == NAMEDATALEN || !string[i])
  244. --i;
  245. len = i;
  246. for (; i >= 0; --i)
  247. new_string[len - i] = string[i];
  248. return new_string;
  249. }
  250. #include "executor/spi.h" /* this is what you need to work with SPI */
  251. #include "commands/trigger.h" /* -"- and triggers */
  252. static TransactionId fd17b_xid = InvalidTransactionId;
  253. static TransactionId fd17a_xid = InvalidTransactionId;
  254. static int fd17b_level = 0;
  255. static int fd17a_level = 0;
  256. static bool fd17b_recursion = true;
  257. static bool fd17a_recursion = true;
  258. HeapTuple funny_dup17(void);
  259. HeapTuple /* have to return HeapTuple to Executor */
  260. funny_dup17()
  261. {
  262. TransactionId *xid;
  263. int    *level;
  264. bool    *recursion;
  265. Relation rel;
  266. TupleDesc tupdesc;
  267. HeapTuple tuple;
  268. char sql[MAX_QUERY_SIZE];
  269. char    *when;
  270. int inserted;
  271. int selected = 0;
  272. int ret;
  273. tuple = CurrentTriggerData->tg_trigtuple;
  274. rel = CurrentTriggerData->tg_relation;
  275. tupdesc = rel->rd_att;
  276. if (TRIGGER_FIRED_BEFORE(CurrentTriggerData->tg_event))
  277. {
  278. xid = &fd17b_xid;
  279. level = &fd17b_level;
  280. recursion = &fd17b_recursion;
  281. when = "BEFORE";
  282. }
  283. else
  284. {
  285. xid = &fd17a_xid;
  286. level = &fd17a_level;
  287. recursion = &fd17a_recursion;
  288. when = "AFTER ";
  289. }
  290. CurrentTriggerData = NULL;
  291. if (!TransactionIdIsCurrentTransactionId(*xid))
  292. {
  293. *xid = GetCurrentTransactionId();
  294. *level = 0;
  295. *recursion = true;
  296. }
  297. if (*level == 17)
  298. {
  299. *recursion = false;
  300. return tuple;
  301. }
  302. if (!(*recursion))
  303. return tuple;
  304. (*level)++;
  305. SPI_connect();
  306. sprintf(sql, "insert into %s select * from %s where %s = '%s'::%s",
  307. SPI_getrelname(rel), SPI_getrelname(rel),
  308. SPI_fname(tupdesc, 1),
  309. SPI_getvalue(tuple, tupdesc, 1),
  310. SPI_gettype(tupdesc, 1));
  311. if ((ret = SPI_exec(sql, 0)) < 0)
  312. elog(ERROR, "funny_dup17 (fired %s) on level %3d: SPI_exec (insert ...) returned %d",
  313.  when, *level, ret);
  314. inserted = SPI_processed;
  315. sprintf(sql, "select count (*) from %s where %s = '%s'::%s",
  316. SPI_getrelname(rel),
  317. SPI_fname(tupdesc, 1),
  318. SPI_getvalue(tuple, tupdesc, 1),
  319. SPI_gettype(tupdesc, 1));
  320. if ((ret = SPI_exec(sql, 0)) < 0)
  321. elog(ERROR, "funny_dup17 (fired %s) on level %3d: SPI_exec (select ...) returned %d",
  322.  when, *level, ret);
  323. if (SPI_processed > 0)
  324. {
  325. selected = int4in(
  326.   SPI_getvalue(
  327.    SPI_tuptable->vals[0],
  328.    SPI_tuptable->tupdesc,
  329.    1
  330.    )
  331. );
  332. }
  333. elog(NOTICE, "funny_dup17 (fired %s) on level %3d: %d/%d tuples inserted/selected",
  334.  when, *level, inserted, selected);
  335. SPI_finish();
  336. (*level)--;
  337. if (*level == 0)
  338. *xid = InvalidTransactionId;
  339. return tuple;
  340. }
  341. HeapTuple ttdummy(void);
  342. int32 set_ttdummy(int32 on);
  343. extern int4 nextval(struct varlena * seqin);
  344. #define TTDUMMY_INFINITY 999999
  345. static void *splan = NULL;
  346. static bool ttoff = false;
  347. HeapTuple
  348. ttdummy()
  349. {
  350. Trigger    *trigger; /* to get trigger name */
  351. char   **args; /* arguments */
  352. int attnum[2]; /* fnumbers of start/stop columns */
  353. Datum oldon,
  354. oldoff;
  355. Datum newon,
  356. newoff;
  357. Datum    *cvals; /* column values */
  358. char    *cnulls; /* column nulls */
  359. char    *relname; /* triggered relation name */
  360. Relation rel; /* triggered relation */
  361. HeapTuple trigtuple;
  362. HeapTuple newtuple = NULL;
  363. HeapTuple rettuple;
  364. TupleDesc tupdesc; /* tuple description */
  365. int natts; /* # of attributes */
  366. bool isnull; /* to know is some column NULL or not */
  367. int ret;
  368. int i;
  369. if (!CurrentTriggerData)
  370. elog(ERROR, "ttdummy: triggers are not initialized");
  371. if (TRIGGER_FIRED_FOR_STATEMENT(CurrentTriggerData->tg_event))
  372. elog(ERROR, "ttdummy: can't process STATEMENT events");
  373. if (TRIGGER_FIRED_AFTER(CurrentTriggerData->tg_event))
  374. elog(ERROR, "ttdummy: must be fired before event");
  375. if (TRIGGER_FIRED_BY_INSERT(CurrentTriggerData->tg_event))
  376. elog(ERROR, "ttdummy: can't process INSERT event");
  377. if (TRIGGER_FIRED_BY_UPDATE(CurrentTriggerData->tg_event))
  378. newtuple = CurrentTriggerData->tg_newtuple;
  379. trigtuple = CurrentTriggerData->tg_trigtuple;
  380. rel = CurrentTriggerData->tg_relation;
  381. relname = SPI_getrelname(rel);
  382. /* check if TT is OFF for this relation */
  383. if (ttoff) /* OFF - nothing to do */
  384. {
  385. pfree(relname);
  386. return (newtuple != NULL) ? newtuple : trigtuple;
  387. }
  388. trigger = CurrentTriggerData->tg_trigger;
  389. if (trigger->tgnargs != 2)
  390. elog(ERROR, "ttdummy (%s): invalid (!= 2) number of arguments %d",
  391.  relname, trigger->tgnargs);
  392. args = trigger->tgargs;
  393. tupdesc = rel->rd_att;
  394. natts = tupdesc->natts;
  395. CurrentTriggerData = NULL;
  396. for (i = 0; i < 2; i++)
  397. {
  398. attnum[i] = SPI_fnumber(tupdesc, args[i]);
  399. if (attnum[i] < 0)
  400. elog(ERROR, "ttdummy (%s): there is no attribute %s", relname, args[i]);
  401. if (SPI_gettypeid(tupdesc, attnum[i]) != INT4OID)
  402. elog(ERROR, "ttdummy (%s): attributes %s and %s must be of abstime type",
  403.  relname, args[0], args[1]);
  404. }
  405. oldon = SPI_getbinval(trigtuple, tupdesc, attnum[0], &isnull);
  406. if (isnull)
  407. elog(ERROR, "ttdummy (%s): %s must be NOT NULL", relname, args[0]);
  408. oldoff = SPI_getbinval(trigtuple, tupdesc, attnum[1], &isnull);
  409. if (isnull)
  410. elog(ERROR, "ttdummy (%s): %s must be NOT NULL", relname, args[1]);
  411. if (newtuple != NULL) /* UPDATE */
  412. {
  413. newon = SPI_getbinval(newtuple, tupdesc, attnum[0], &isnull);
  414. if (isnull)
  415. elog(ERROR, "ttdummy (%s): %s must be NOT NULL", relname, args[0]);
  416. newoff = SPI_getbinval(newtuple, tupdesc, attnum[1], &isnull);
  417. if (isnull)
  418. elog(ERROR, "ttdummy (%s): %s must be NOT NULL", relname, args[1]);
  419. if (oldon != newon || oldoff != newoff)
  420. elog(ERROR, "ttdummy (%s): you can't change %s and/or %s columns (use set_ttdummy)",
  421.  relname, args[0], args[1]);
  422. if (newoff != TTDUMMY_INFINITY)
  423. {
  424. pfree(relname); /* allocated in upper executor context */
  425. return NULL;
  426. }
  427. }
  428. else if (oldoff != TTDUMMY_INFINITY) /* DELETE */
  429. {
  430. pfree(relname);
  431. return NULL;
  432. }
  433. {
  434. struct varlena *seqname = textin("ttdummy_seq");
  435. newoff = nextval(seqname);
  436. pfree(seqname);
  437. }
  438. /* Connect to SPI manager */
  439. if ((ret = SPI_connect()) < 0)
  440. elog(ERROR, "ttdummy (%s): SPI_connect returned %d", relname, ret);
  441. /* Fetch tuple values and nulls */
  442. cvals = (Datum *) palloc(natts * sizeof(Datum));
  443. cnulls = (char *) palloc(natts * sizeof(char));
  444. for (i = 0; i < natts; i++)
  445. {
  446. cvals[i] = SPI_getbinval((newtuple != NULL) ? newtuple : trigtuple,
  447.  tupdesc, i + 1, &isnull);
  448. cnulls[i] = (isnull) ? 'n' : ' ';
  449. }
  450. /* change date column(s) */
  451. if (newtuple) /* UPDATE */
  452. {
  453. cvals[attnum[0] - 1] = newoff; /* start_date eq current date */
  454. cnulls[attnum[0] - 1] = ' ';
  455. cvals[attnum[1] - 1] = TTDUMMY_INFINITY; /* stop_date eq INFINITY */
  456. cnulls[attnum[1] - 1] = ' ';
  457. }
  458. else
  459. /* DELETE */
  460. {
  461. cvals[attnum[1] - 1] = newoff; /* stop_date eq current date */
  462. cnulls[attnum[1] - 1] = ' ';
  463. }
  464. /* if there is no plan ... */
  465. if (splan == NULL)
  466. {
  467. void    *pplan;
  468. Oid    *ctypes;
  469. char sql[MAX_QUERY_SIZE];
  470. /* allocate ctypes for preparation */
  471. ctypes = (Oid *) palloc(natts * sizeof(Oid));
  472. /*
  473.  * Construct query: INSERT INTO _relation_ VALUES ($1, ...)
  474.  */
  475. sprintf(sql, "INSERT INTO %s VALUES (", relname);
  476. for (i = 1; i <= natts; i++)
  477. {
  478. sprintf(sql + strlen(sql), "$%d%s",
  479. i, (i < natts) ? ", " : ")");
  480. ctypes[i - 1] = SPI_gettypeid(tupdesc, i);
  481. }
  482. /* Prepare plan for query */
  483. pplan = SPI_prepare(sql, natts, ctypes);
  484. if (pplan == NULL)
  485. elog(ERROR, "ttdummy (%s): SPI_prepare returned %d", relname, SPI_result);
  486. pplan = SPI_saveplan(pplan);
  487. if (pplan == NULL)
  488. elog(ERROR, "ttdummy (%s): SPI_saveplan returned %d", relname, SPI_result);
  489. splan = pplan;
  490. }
  491. ret = SPI_execp(splan, cvals, cnulls, 0);
  492. if (ret < 0)
  493. elog(ERROR, "ttdummy (%s): SPI_execp returned %d", relname, ret);
  494. /* Tuple to return to upper Executor ... */
  495. if (newtuple) /* UPDATE */
  496. {
  497. HeapTuple tmptuple;
  498. tmptuple = SPI_copytuple(trigtuple);
  499. rettuple = SPI_modifytuple(rel, tmptuple, 1, &(attnum[1]), &newoff, NULL);
  500. SPI_pfree(tmptuple);
  501. }
  502. else
  503. /* DELETE */
  504. rettuple = trigtuple;
  505. SPI_finish(); /* don't forget say Bye to SPI mgr */
  506. pfree(relname);
  507. return rettuple;
  508. }
  509. int32
  510. set_ttdummy(int32 on)
  511. {
  512. if (ttoff) /* OFF currently */
  513. {
  514. if (on == 0)
  515. return 0;
  516. /* turn ON */
  517. ttoff = false;
  518. return 0;
  519. }
  520. /* ON currently */
  521. if (on != 0)
  522. return 1;
  523. /* turn OFF */
  524. ttoff = true;
  525. return 1;
  526. }