| 1 |
/* |
| 2 |
* SLV: Ascend Numeric Solver |
| 3 |
* by Karl Michael Westerberg |
| 4 |
* Created: 2/6/90 |
| 5 |
* Version: $Revision: 1.31 $ |
| 6 |
* Version control file: $RCSfile: var.c,v $ |
| 7 |
* Date last modified: $Date: 1998/02/19 13:31:36 $ |
| 8 |
* Last modified by: $Author: mthomas $ |
| 9 |
* |
| 10 |
* This file is part of the SLV solver. |
| 11 |
* |
| 12 |
* Copyright (C) 1990 Karl Michael Westerberg |
| 13 |
* Copyright (C) 1993 Joseph Zaher |
| 14 |
* Copyright (C) 1994 Joseph Zaher, Benjamin Andrew Allan |
| 15 |
* |
| 16 |
* The SLV solver is free software; you can redistribute |
| 17 |
* it and/or modify it under the terms of the GNU General Public License as |
| 18 |
* published by the Free Software Foundation; either version 2 of the |
| 19 |
* License, or (at your option) any later version. |
| 20 |
* |
| 21 |
* The SLV solver is distributed in hope that it will be |
| 22 |
* useful, but WITHOUT ANY WARRANTY; without even the implied warranty of |
| 23 |
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU |
| 24 |
* General Public License for more details. |
| 25 |
* |
| 26 |
* You should have received a copy of the GNU General Public License |
| 27 |
* along with the program; if not, write to the Free Software Foundation, |
| 28 |
* Inc., 675 Mass Ave, Cambridge, MA 02139 USA. Check the file named |
| 29 |
* COPYING. COPYING is found in ../compiler. |
| 30 |
* |
| 31 |
*/ |
| 32 |
|
| 33 |
#include <utilities/ascConfig.h> |
| 34 |
#include <utilities/ascMalloc.h> |
| 35 |
#include <general/dstring.h> |
| 36 |
#include <general/list.h> |
| 37 |
#include <compiler/compiler.h> |
| 38 |
#include <compiler/symtab.h> |
| 39 |
#include <compiler/instance_enum.h> |
| 40 |
#include <compiler/fractions.h> |
| 41 |
#include <compiler/module.h> |
| 42 |
#include <compiler/library.h> |
| 43 |
#include <compiler/dimen.h> |
| 44 |
#include <compiler/child.h> |
| 45 |
#include <compiler/type_desc.h> |
| 46 |
#include <compiler/atomvalue.h> |
| 47 |
#include <compiler/parentchild.h> |
| 48 |
#include <compiler/instquery.h> |
| 49 |
#include <compiler/instance_io.h> |
| 50 |
#include "mtx.h" |
| 51 |
#include "slv_types.h" |
| 52 |
#include "rel.h" |
| 53 |
#include "var.h" |
| 54 |
#include "discrete.h" |
| 55 |
#include "conditional.h" |
| 56 |
#include "logrel.h" |
| 57 |
#include "bnd.h" |
| 58 |
#include "slv_server.h" |
| 59 |
#include "slv_common.h" |
| 60 |
#include "linsol.h" |
| 61 |
#include "linsolqr.h" |
| 62 |
#include "slv_client.h" |
| 63 |
|
| 64 |
/* useful cast */ |
| 65 |
#define IPTR(i) ((struct Instance *)(i)) |
| 66 |
|
| 67 |
/* useful symbol table things to know */ |
| 68 |
#define FIXED_V g_strings[0] |
| 69 |
#define LOWER_V g_strings[1] |
| 70 |
#define UPPER_V g_strings[2] |
| 71 |
#define RELAXED_V g_strings[3] |
| 72 |
#define NOMINAL_V g_strings[4] |
| 73 |
#define INTERFACE_V g_strings[5] |
| 74 |
#define ODEATOL_V g_strings[6] |
| 75 |
|
| 76 |
/* |
| 77 |
* array of those symbol table entries we need. |
| 78 |
*/ |
| 79 |
static symchar * g_strings[7]; |
| 80 |
|
| 81 |
SlvBackendToken var_instanceF(const struct var_variable *var) |
| 82 |
{ if (var==NULL || var->ratom==NULL) { |
| 83 |
FPRINTF(stderr,"var_instance called on bad var\n"); |
| 84 |
return NULL; |
| 85 |
} |
| 86 |
return var->ratom; |
| 87 |
} |
| 88 |
|
| 89 |
void var_set_instanceF(struct var_variable *var, SlvBackendToken i) |
| 90 |
{ |
| 91 |
if (var==NULL) { |
| 92 |
FPRINTF(stderr,"var_set_instance called on NULL var\n"); |
| 93 |
return; |
| 94 |
} |
| 95 |
var->ratom = i; |
| 96 |
} |
| 97 |
|
| 98 |
|
| 99 |
char *var_make_name(const slv_system_t sys,const struct var_variable *var) |
| 100 |
{ |
| 101 |
return WriteInstanceNameString(IPTR(var->ratom),IPTR(slv_instance(sys))); |
| 102 |
} |
| 103 |
|
| 104 |
char *var_make_xname(const struct var_variable *var) |
| 105 |
{ |
| 106 |
static char name[81]; |
| 107 |
char *res; |
| 108 |
sprintf(name,"x%d",var_sindex(var)); |
| 109 |
res=ASC_NEW_ARRAY(char,strlen(name)+1); |
| 110 |
sprintf(res,"%s",name); |
| 111 |
return res; |
| 112 |
} |
| 113 |
|
| 114 |
void var_write_name(const slv_system_t sys, |
| 115 |
const struct var_variable *var,FILE *fp) |
| 116 |
{ |
| 117 |
if (var == NULL || fp==NULL) return; |
| 118 |
if (sys!=NULL) { |
| 119 |
WriteInstanceName(fp,var_instance(var),slv_instance(sys)); |
| 120 |
} else { |
| 121 |
WriteInstanceName(fp,var_instance(var),NULL); |
| 122 |
} |
| 123 |
} |
| 124 |
|
| 125 |
void var_destroy(struct var_variable *var) |
| 126 |
{ |
| 127 |
if (var==NULL) return; |
| 128 |
var->ratom = NULL; |
| 129 |
ascfree((POINTER)var->incidence); |
| 130 |
} |
| 131 |
|
| 132 |
|
| 133 |
int32 var_mindexF(const struct var_variable *var) |
| 134 |
{ |
| 135 |
if (var==NULL || var->ratom==NULL) { |
| 136 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"bad var"); |
| 137 |
return -1; |
| 138 |
} |
| 139 |
return var->mindex; |
| 140 |
} |
| 141 |
|
| 142 |
void var_set_mindexF(struct var_variable *var, int32 mindex) |
| 143 |
{ |
| 144 |
if (var==NULL || var->ratom==NULL) { |
| 145 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"bad var"); |
| 146 |
return; |
| 147 |
} |
| 148 |
var->mindex = mindex; |
| 149 |
} |
| 150 |
|
| 151 |
int32 var_sindexF(const struct var_variable *var) |
| 152 |
{ |
| 153 |
if (var==NULL || var->ratom==NULL) { |
| 154 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"bad var"); |
| 155 |
return -1; |
| 156 |
} |
| 157 |
return var->sindex; |
| 158 |
} |
| 159 |
|
| 160 |
void var_set_sindexF(struct var_variable *var, int32 sindex) |
| 161 |
{ |
| 162 |
if (var==NULL || var->ratom==NULL) { |
| 163 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"bad var"); |
| 164 |
return; |
| 165 |
} |
| 166 |
var->sindex = sindex; |
| 167 |
} |
| 168 |
|
| 169 |
real64 var_value(const struct var_variable *var) |
| 170 |
{ |
| 171 |
if (var==NULL || var->ratom==NULL) { |
| 172 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"bad var"); |
| 173 |
return 0.0; |
| 174 |
} |
| 175 |
return( RealAtomValue(var->ratom) ); |
| 176 |
} |
| 177 |
|
| 178 |
void var_set_value(struct var_variable *var, real64 value) |
| 179 |
{ |
| 180 |
if (var==NULL || var->ratom==NULL) { |
| 181 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"bad var"); |
| 182 |
return; |
| 183 |
} |
| 184 |
SetRealAtomValue(var->ratom,value,(unsigned)0); |
| 185 |
} |
| 186 |
|
| 187 |
real64 var_nominal(struct var_variable *var) |
| 188 |
{ |
| 189 |
struct Instance *c; |
| 190 |
if (var==NULL || var->ratom==NULL) { |
| 191 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"bad var"); |
| 192 |
return 1.0; |
| 193 |
} |
| 194 |
c = ChildByChar(var->ratom,NOMINAL_V); |
| 195 |
if( c == NULL ) { |
| 196 |
FPRINTF(ASCERR,"no 'nominal' field in variable"); |
| 197 |
/* WriteInstance(stderr,IPTR(var->ratom)); */ |
| 198 |
return 1.0; |
| 199 |
} |
| 200 |
return( RealAtomValue(c) ); |
| 201 |
} |
| 202 |
|
| 203 |
void var_set_nominal(struct var_variable *var, real64 nominal) |
| 204 |
{ |
| 205 |
struct Instance *c; |
| 206 |
if (var==NULL || var->ratom==NULL) { |
| 207 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"bad var"); |
| 208 |
return; |
| 209 |
} |
| 210 |
c = ChildByChar(IPTR(var->ratom),NOMINAL_V); |
| 211 |
if( c == NULL ) { |
| 212 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"no 'nominal' field in var"); |
| 213 |
/* WriteInstance(stderr,IPTR(var->ratom)); */ |
| 214 |
return; |
| 215 |
} |
| 216 |
SetRealAtomValue(c,nominal,(unsigned)0); |
| 217 |
} |
| 218 |
|
| 219 |
|
| 220 |
real64 var_lower_bound(struct var_variable *var) |
| 221 |
{ |
| 222 |
struct Instance *c; |
| 223 |
if (var==NULL || var->ratom==NULL) { |
| 224 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"bad var"); |
| 225 |
return 0.0; |
| 226 |
} |
| 227 |
c = ChildByChar(IPTR(var->ratom),LOWER_V); |
| 228 |
if( c == NULL ) { |
| 229 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"no 'lower_bound' field"); |
| 230 |
WriteInstance(stderr,IPTR(var->ratom)); |
| 231 |
return 0.0; |
| 232 |
} |
| 233 |
return( RealAtomValue(c) ); |
| 234 |
} |
| 235 |
|
| 236 |
void var_set_lower_bound(struct var_variable *var, real64 lower_bound) |
| 237 |
{ |
| 238 |
struct Instance *c; |
| 239 |
if (var==NULL || var->ratom==NULL) { |
| 240 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"bad var"); |
| 241 |
return; |
| 242 |
} |
| 243 |
c = ChildByChar(IPTR(var->ratom),LOWER_V); |
| 244 |
if( c == NULL ) { |
| 245 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"no 'lower_bound' field"); |
| 246 |
/* WriteInstance(stderr,IPTR(var->ratom)); */ |
| 247 |
return; |
| 248 |
} |
| 249 |
SetRealAtomValue(c,lower_bound,(unsigned)0); |
| 250 |
} |
| 251 |
|
| 252 |
|
| 253 |
real64 var_upper_bound(struct var_variable *var) |
| 254 |
{ |
| 255 |
struct Instance *c; |
| 256 |
if (var==NULL || var->ratom==NULL) { |
| 257 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"bad var"); |
| 258 |
return 0.0; |
| 259 |
} |
| 260 |
c = ChildByChar(IPTR(var->ratom),UPPER_V); |
| 261 |
if( c == NULL ) { |
| 262 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"no 'upper_bound' field"); |
| 263 |
/* WriteInstance(stderr,IPTR(var->ratom)); */ |
| 264 |
return 0.0; |
| 265 |
} |
| 266 |
return( RealAtomValue(c) ); |
| 267 |
} |
| 268 |
|
| 269 |
void var_set_upper_bound(struct var_variable *var, real64 upper_bound) |
| 270 |
{ |
| 271 |
struct Instance *c; |
| 272 |
if (var==NULL || var->ratom==NULL) { |
| 273 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"bad var"); |
| 274 |
return; |
| 275 |
} |
| 276 |
c = ChildByChar(IPTR(var->ratom),UPPER_V); |
| 277 |
if( c == NULL ) { |
| 278 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"no 'upper_bound' field"); |
| 279 |
/* WriteInstance(stderr,IPTR(var->ratom)); */ |
| 280 |
return; |
| 281 |
} |
| 282 |
SetRealAtomValue(c,upper_bound,(unsigned)0); |
| 283 |
} |
| 284 |
|
| 285 |
double var_odeatol(struct var_variable *var){ |
| 286 |
struct Instance *c; |
| 287 |
if(var==NULL||var->ratom==NULL){ |
| 288 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"bad var"); |
| 289 |
return -1; |
| 290 |
} |
| 291 |
c = ChildByChar(IPTR(var->ratom),ODEATOL_V); |
| 292 |
if(c==NULL){ |
| 293 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"no '%s' field",SCP(ODEATOL_V)); |
| 294 |
return -1; |
| 295 |
} |
| 296 |
return RealAtomValue(c); |
| 297 |
} |
| 298 |
|
| 299 |
|
| 300 |
uint32 var_flagsF(const struct var_variable *var) |
| 301 |
{ |
| 302 |
return var->flags; |
| 303 |
} |
| 304 |
|
| 305 |
void var_set_flagsF(struct var_variable *var, uint32 flags) |
| 306 |
{ |
| 307 |
var->flags = flags; |
| 308 |
} |
| 309 |
|
| 310 |
uint32 var_fixed(struct var_variable *var) |
| 311 |
{ |
| 312 |
struct Instance *c; |
| 313 |
if (var==NULL || var->ratom==NULL) { |
| 314 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"bad var"); |
| 315 |
return FALSE; |
| 316 |
} |
| 317 |
c = ChildByChar(IPTR(var->ratom),FIXED_V); |
| 318 |
if( c == NULL ) { |
| 319 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"no 'fixed' field"); |
| 320 |
/* WriteInstance(stderr,IPTR(var->ratom)); */ |
| 321 |
return FALSE; |
| 322 |
} |
| 323 |
var_set_flagbit(var,VAR_FIXED,GetBooleanAtomValue(c)); |
| 324 |
return( GetBooleanAtomValue(c) ); |
| 325 |
} |
| 326 |
|
| 327 |
void var_set_fixed(struct var_variable *var, uint32 fixed) |
| 328 |
{ |
| 329 |
struct Instance *c; |
| 330 |
if (var==NULL || var->ratom==NULL) { |
| 331 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"bad var"); |
| 332 |
return; |
| 333 |
} |
| 334 |
c = ChildByChar(IPTR(var->ratom),FIXED_V); |
| 335 |
if( c == NULL ) { |
| 336 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"no 'fixed' field"); |
| 337 |
/* WriteInstance(stderr,IPTR(var->ratom)); */ |
| 338 |
return; |
| 339 |
} |
| 340 |
SetBooleanAtomValue(c,fixed,(unsigned)0); |
| 341 |
var_set_flagbit(var,VAR_FIXED,fixed); |
| 342 |
} |
| 343 |
|
| 344 |
uint32 var_relaxed(struct var_variable *var) |
| 345 |
{ |
| 346 |
struct Instance *c; |
| 347 |
if (var==NULL || var->ratom==NULL) { |
| 348 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"bad var"); |
| 349 |
return FALSE; |
| 350 |
} |
| 351 |
c = ChildByChar((var->ratom),RELAXED_V); |
| 352 |
if( c == NULL ) { |
| 353 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"no 'relaxed' field"); |
| 354 |
/* WriteInstance(stderr,(var->ratom)); */ |
| 355 |
return FALSE; |
| 356 |
} |
| 357 |
var_set_flagbit(var,VAR_RELAXED,GetBooleanAtomValue(c)); |
| 358 |
return( GetBooleanAtomValue(c) ); |
| 359 |
} |
| 360 |
|
| 361 |
void var_set_relaxed(struct var_variable *var, uint32 fixed) |
| 362 |
{ |
| 363 |
struct Instance *c; |
| 364 |
if (var==NULL || var->ratom==NULL) { |
| 365 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"bad var"); |
| 366 |
return; |
| 367 |
} |
| 368 |
c = ChildByChar(IPTR(var->ratom),RELAXED_V); |
| 369 |
if( c == NULL ) { |
| 370 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"no 'relaxed' field"); |
| 371 |
/* WriteInstance(stderr,IPTR(var->ratom)); */ |
| 372 |
return; |
| 373 |
} |
| 374 |
SetBooleanAtomValue(c,fixed,(unsigned)0); |
| 375 |
var_set_flagbit(var,VAR_RELAXED,fixed); |
| 376 |
} |
| 377 |
|
| 378 |
uint32 var_interface(struct var_variable *var) |
| 379 |
{ |
| 380 |
struct Instance *c; |
| 381 |
if (var==NULL || var->ratom==NULL) { |
| 382 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"bad var"); |
| 383 |
return FALSE; |
| 384 |
} |
| 385 |
c = ChildByChar(IPTR(var->ratom),INTERFACE_V); |
| 386 |
if( c == NULL ) { |
| 387 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"no 'interface' field"); |
| 388 |
/* WriteInstance(stderr,IPTR(var->ratom)); */ |
| 389 |
return 0; |
| 390 |
} |
| 391 |
var_set_flagbit(var,VAR_INTERFACE,GetBooleanAtomValue(c)); |
| 392 |
return( GetIntegerAtomValue(c) ); |
| 393 |
} |
| 394 |
|
| 395 |
extern uint32 var_flagbit(const struct var_variable *var,const uint32 one) |
| 396 |
{ |
| 397 |
if (var==NULL || var->ratom == NULL) { |
| 398 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"bad var"); |
| 399 |
return 0; |
| 400 |
} |
| 401 |
return (var->flags & one); |
| 402 |
} |
| 403 |
|
| 404 |
void var_set_flagbit(struct var_variable *var, uint32 field,uint32 one) |
| 405 |
{ |
| 406 |
if (var==NULL || var->ratom == NULL) { |
| 407 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"bad var"); |
| 408 |
return; |
| 409 |
} |
| 410 |
if (one) { |
| 411 |
var->flags |= field; |
| 412 |
} else { |
| 413 |
var->flags &= ~field; |
| 414 |
} |
| 415 |
} |
| 416 |
|
| 417 |
int32 var_apply_filter(const struct var_variable *var, |
| 418 |
const var_filter_t *filter |
| 419 |
){ |
| 420 |
if (var==NULL || filter==NULL || var->ratom == NULL) { |
| 421 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"'%s' miscalled with NULL",__FUNCTION__); |
| 422 |
return FALSE; |
| 423 |
} |
| 424 |
/* AND to mask off irrelevant bits in flags and match value, then compare */ |
| 425 |
return (filter->matchbits & var->flags) == (filter->matchbits & filter->matchvalue); |
| 426 |
} |
| 427 |
|
| 428 |
/** |
| 429 |
Implementation function for var_n_incidences(). Do not call this |
| 430 |
function directly - use var_n_incidences() instead. |
| 431 |
*/ |
| 432 |
int32 var_n_incidencesF(struct var_variable *var) |
| 433 |
{ |
| 434 |
if (var!=NULL) return var->n_incidences; |
| 435 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"NULL var"); |
| 436 |
return 0; |
| 437 |
} |
| 438 |
|
| 439 |
void var_set_incidencesF(struct var_variable *var,int32 n, |
| 440 |
struct rel_relation **i) |
| 441 |
{ |
| 442 |
if(var!=NULL && n >=0) { |
| 443 |
if (n && i==NULL) { |
| 444 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"NULL i"); |
| 445 |
} |
| 446 |
var->n_incidences = n; |
| 447 |
var->incidence = i; |
| 448 |
return; |
| 449 |
} |
| 450 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"NULL var, or n < 0"); |
| 451 |
} |
| 452 |
const struct rel_relation **var_incidence_list( struct var_variable *var) |
| 453 |
{ |
| 454 |
if (var==NULL) return NULL; |
| 455 |
return( (const struct rel_relation **)var->incidence ); |
| 456 |
} |
| 457 |
|
| 458 |
struct rel_relation **var_incidence_list_to_modify( struct var_variable *var) |
| 459 |
{ |
| 460 |
if (var==NULL) return NULL; |
| 461 |
return( (struct rel_relation **)var->incidence ); |
| 462 |
} |
| 463 |
|
| 464 |
|
| 465 |
|
| 466 |
/* |
| 467 |
* global for use with the push function. Sets the ip to the |
| 468 |
* value in g_var_tag; |
| 469 |
* should be using vp instead of a global counter. |
| 470 |
*/ |
| 471 |
static void *g_var_tag = NULL; |
| 472 |
static void * SetVarTags(struct Instance *i,VOIDPTR vp) |
| 473 |
{ |
| 474 |
(void)vp; |
| 475 |
if (i!=NULL && InstanceKind(i)==REAL_ATOM_INST) { |
| 476 |
return g_var_tag; |
| 477 |
} else { |
| 478 |
return NULL; |
| 479 |
} |
| 480 |
} |
| 481 |
|
| 482 |
struct var_variable **var_BackendTokens_to_vars(slv_system_t sys, |
| 483 |
SlvBackendToken *atoms, int32 len) |
| 484 |
{ |
| 485 |
int32 i,vartot,vlen,count=0; |
| 486 |
uint32 apos,ulen; |
| 487 |
struct var_variable **result; |
| 488 |
struct var_variable **vlist; |
| 489 |
struct gl_list_t *oldips; |
| 490 |
if (sys==NULL || atoms == NULL || len < 1) { |
| 491 |
return NULL; |
| 492 |
} |
| 493 |
ulen = (uint32)len; |
| 494 |
result = (struct var_variable **)malloc(len*sizeof(struct var_variable *)); |
| 495 |
if (result == NULL) return result; |
| 496 |
/* init results to null */ |
| 497 |
for (i=0; i<len; i++) result[i] = NULL; |
| 498 |
/* fill ips w/len in all the vars in tree. */ |
| 499 |
g_var_tag = (void *)len; |
| 500 |
vartot = slv_get_num_master_vars(sys) + |
| 501 |
slv_get_num_master_pars(sys) + |
| 502 |
slv_get_num_master_unattached(sys); |
| 503 |
oldips = PushInterfacePtrs(slv_instance(sys),SetVarTags,vartot,0,NULL); |
| 504 |
/* fill ips of wanted atoms with i */ |
| 505 |
for (i=0; i<len; i++) { |
| 506 |
if (GetInterfacePtr(atoms[i])==g_var_tag && |
| 507 |
InstanceKind(atoms[i]) == REAL_ATOM_INST) { |
| 508 |
/* guard a little */ |
| 509 |
SetInterfacePtr((struct Instance *)atoms[i],(void *)i); |
| 510 |
} else { |
| 511 |
/* the odds of g_var_tag being a legal pointer are vanishingly |
| 512 |
small, so if we find an ATOM without g_var_tag we assume it |
| 513 |
is outside the tree and shouldn't have been in the list. */ |
| 514 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"bad token"); |
| 515 |
} |
| 516 |
} |
| 517 |
/* run through the master lists and put the vars with their atoms */ |
| 518 |
vlist = slv_get_master_var_list(sys); |
| 519 |
vlen = slv_get_num_master_vars(sys); |
| 520 |
for (i = 0; i <vlen; i++) { |
| 521 |
apos = (uint32)GetInterfacePtr(var_instance(vlist[i])); |
| 522 |
if ( apos < ulen ) { |
| 523 |
result[apos] = vlist[i]; |
| 524 |
count++; |
| 525 |
} |
| 526 |
} |
| 527 |
vlist = slv_get_master_par_list(sys); |
| 528 |
vlen = slv_get_num_master_pars(sys); |
| 529 |
for (i = 0; i <vlen; i++) { |
| 530 |
apos = (uint32)GetInterfacePtr(var_instance(vlist[i])); |
| 531 |
if ( apos < ulen ) { |
| 532 |
result[apos] = vlist[i]; |
| 533 |
count++; |
| 534 |
} |
| 535 |
} |
| 536 |
vlist = slv_get_master_unattached_list(sys); |
| 537 |
vlen = slv_get_num_master_unattached(sys); |
| 538 |
for (i = 0; i <vlen; i++) { |
| 539 |
apos = (uint32)GetInterfacePtr(var_instance(vlist[i])); |
| 540 |
if ( apos < ulen ) { |
| 541 |
result[apos] = vlist[i]; |
| 542 |
count++; |
| 543 |
} |
| 544 |
} |
| 545 |
if (count < len) { |
| 546 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"found less than expected vars (len = %d, found = %d)",len,count); |
| 547 |
} else { |
| 548 |
ERROR_REPORTER_HERE(ASC_PROG_ERR,"found more than expected vars (len = %d, found = %d)",len,count); |
| 549 |
} |
| 550 |
PopInterfacePtrs(oldips,NULL,NULL); |
| 551 |
return result; |
| 552 |
} |
| 553 |
|
| 554 |
static struct TypeDescription *g_solver_var_type; |
| 555 |
static struct TypeDescription *g_solver_int_type; |
| 556 |
static struct TypeDescription *g_solver_binary_type; |
| 557 |
static struct TypeDescription *g_solver_semi_type; |
| 558 |
|
| 559 |
boolean set_solver_types(void) { |
| 560 |
boolean nerr = 0; |
| 561 |
if( (g_solver_var_type = FindType(AddSymbol(SOLVER_VAR_STR))) == NULL ) { |
| 562 |
ERROR_REPORTER_HERE(ASC_PROG_FATAL,"'solver_var' not defined"); |
| 563 |
nerr++; |
| 564 |
} |
| 565 |
if( (g_solver_int_type = FindType(AddSymbol(SOLVER_INT_STR))) == NULL ) { |
| 566 |
ERROR_REPORTER_HERE(ASC_PROG_ERROR,"'solver_int' not defined: MPS will not work"); |
| 567 |
nerr++; |
| 568 |
} |
| 569 |
g_solver_binary_type = FindType(AddSymbol(SOLVER_BINARY_STR)); |
| 570 |
if( g_solver_binary_type == NULL) { |
| 571 |
ERROR_REPORTER_HERE(ASC_PROG_FATAL,"'solver_binary' not defined: MPS will not work"); |
| 572 |
nerr++; |
| 573 |
} |
| 574 |
if( (g_solver_semi_type = FindType(AddSymbol(SOLVER_SEMI_STR))) == NULL ) { |
| 575 |
ERROR_REPORTER_HERE(ASC_PROG_FATAL,"'solver_semi' not defined: MPS will not work"); |
| 576 |
nerr++; |
| 577 |
} |
| 578 |
|
| 579 |
LOWER_V = AddSymbolL("lower_bound",11); |
| 580 |
UPPER_V = AddSymbolL("upper_bound",11); |
| 581 |
RELAXED_V = AddSymbolL("relaxed",7); |
| 582 |
NOMINAL_V = AddSymbolL("nominal",7); |
| 583 |
FIXED_V = AddSymbolL("fixed",5); |
| 584 |
INTERFACE_V = AddSymbolL("interface",9); |
| 585 |
ODEATOL_V = AddSymbol("ode_atol"); |
| 586 |
return nerr; |
| 587 |
} |
| 588 |
|
| 589 |
boolean solver_var( SlvBackendToken inst) |
| 590 |
{ |
| 591 |
struct TypeDescription *type; |
| 592 |
|
| 593 |
if (!g_solver_var_type) return FALSE; |
| 594 |
type = InstanceTypeDesc(IPTR(inst)); |
| 595 |
return( type == MoreRefined(type,g_solver_var_type) ); |
| 596 |
} |
| 597 |
|
| 598 |
boolean solver_int( SlvBackendToken inst) |
| 599 |
{ |
| 600 |
struct TypeDescription *type; |
| 601 |
|
| 602 |
if (!g_solver_int_type) return FALSE; |
| 603 |
type = InstanceTypeDesc(IPTR(inst)); |
| 604 |
return( type == MoreRefined(type,g_solver_int_type) ); |
| 605 |
} |
| 606 |
|
| 607 |
boolean solver_binary( SlvBackendToken inst) |
| 608 |
{ |
| 609 |
struct TypeDescription *type; |
| 610 |
|
| 611 |
if (!g_solver_binary_type) return FALSE; |
| 612 |
type = InstanceTypeDesc(IPTR(inst)); |
| 613 |
return( type == MoreRefined(type,g_solver_binary_type) ); |
| 614 |
} |
| 615 |
|
| 616 |
boolean solver_semi( SlvBackendToken inst) |
| 617 |
{ |
| 618 |
struct TypeDescription *type; |
| 619 |
|
| 620 |
if (!g_solver_semi_type) return FALSE; |
| 621 |
type = InstanceTypeDesc(IPTR(inst)); |
| 622 |
return( type == MoreRefined(type,g_solver_semi_type) ); |
| 623 |
} |
| 624 |
|
Parent Directory
| 
Revision Log