+ trinode_list_iterator tn_current, tn_last;
+ tn_current = node_list.begin();
+ tn_last = node_list.end();
+ int counter = 0;
+ for ( ; tn_current != tn_last; ++tn_current ) {
+ in.pointlist[counter++] = tn_current->x();
+ in.pointlist[counter++] = tn_current->y();
+ }
+
+ in.pointattributelist = (REAL *) NULL;
+ in.pointmarkerlist = (int *) NULL;
+
+ // segment list
+ in.numberofsegments = 0;
+
+ tripoly_list_iterator tp_current, tp_last;
+ for ( int i = 0; i < FG_MAX_AREA_TYPES; ++i ) {
+ cout << "area type = " << i << endl;
+ tp_current = polylist[i].begin();
+ tp_last = polylist[i].end();
+ for ( ; tp_current != tp_last; ++tp_current ) {
+ poly = *tp_current;
+ in.numberofsegments += poly.size() + 1;
+ }
+ }
+
+ in.numberofsegments = 0;
+
+ in.numberofholes = 0;
+ in.numberofregions = 1;
+ in.regionlist = (REAL *) malloc(in.numberofregions * 4 * sizeof(REAL));
+ in.regionlist[0] = 0.5;
+ in.regionlist[1] = 5.0;
+ in.regionlist[2] = 7.0; /* Regional attribute (for whole mesh). */
+ in.regionlist[3] = 0.1; /* Area constraint that will not be used. */
+
+ /*
+ tripoly_list_iterator current, last;