Changeset 255 for trunk/ADOL-C


Ignore:
Timestamp:
Jun 30, 2011 9:20:56 AM (9 years ago)
Author:
awalther
Message:

change nlf to csod due to overwrites, different handling required

Location:
trunk/ADOL-C/src
Files:
2 deleted
3 edited

Legend:

Unmodified
Added
Removed
  • trunk/ADOL-C/src/Makefile.am

    r247 r255  
    6060                       int_forward_s.c int_forward_t.c \
    6161                       indopro_forward_s.c indopro_forward_t.c \
    62                        nlf_forward_s.c nlf_forward_t.c \
     62                       csod_forward_s.c csod_forward_t.c \
    6363                       nonl_ind_forward_s.c nonl_ind_forward_t.c \
    6464                       int_reverse_s.c int_reverse_t.c
  • trunk/ADOL-C/src/interfaces.h

    r247 r255  
    316316/*                                             NONLINEAR FRONTIER UTILITIES */
    317317/*--------------------------------------------------------------------------*/
    318 /*                                                                NLF, SAFE */
    319 /* nlf_forward_safe(tag, m, n, p, x[n], *crs[m])                            */
    320 
    321 ADOLC_DLL_EXPORT int nlf_forward_safe
     318/*                                                               CSOD, SAFE */
     319/* csod_forward_safe(tag, m, n, p, x[n], *crs[m])                            */
     320
     321ADOLC_DLL_EXPORT int csod_forward_safe
    322322(short, int, int, const double*, unsigned int**);
    323323
     
    325325/*--------------------------------------------------------------------------*/
    326326/*                                                           INDOPRO, TIGHT */
    327 /* nlf_forward_tight(tag, m, n,  x[n], *crs[m])                             */
    328 
    329 ADOLC_DLL_EXPORT int nlf_forward_tight
     327/* csod_forward_tight(tag, m, n,  x[n], *crs[m])                             */
     328
     329ADOLC_DLL_EXPORT int csod_forward_tight
    330330(short, int, int, const double*, unsigned int**);
    331331
  • trunk/ADOL-C/src/uni5_for.c

    r247 r255  
    141141#endif
    142142
    143 #if defined(_NLF_)
     143#if defined(_CSOD_)
    144144#if defined(_TIGHT_)
    145 #define GENERATED_FILENAME "nlf_forward_t"
     145#define GENERATED_FILENAME "csod_forward_t"
    146146#endif
    147147#if defined(_NTIGHT_)
    148 #define GENERATED_FILENAME "nlf_forward_s"
     148#define GENERATED_FILENAME "csod_forward_s"
    149149#endif
    150150typedef struct IndexElement {
     
    594594#endif
    595595#else
    596 #if defined(_NLF_)
     596#if defined(_CSOD_)
    597597#if defined(_TIGHT_)
    598598/****************************************************************************/
    599599/* First Order Vector version of the forward mode for bit patterns, tight   */
    600600/****************************************************************************/
    601 int nlf_forward_tight(
     601int csod_forward_tight(
    602602    short             tnum,        /* tape id                              */
    603603    int               depcheck,    /* consistency chk on # of dependents   */
     
    611611/* First Order Vector version of the forward mode, bit pattern, safe        */
    612612/****************************************************************************/
    613 int nlf_forward_safe(
     613int csod_forward_safe(
    614614    short             tnum,        /* tape id                              */
    615615    int               depcheck,    /* consistency chk on # of dependents   */
     
    779779    int l=0;
    780780    int max_ind_dom;
    781 #if defined (_NLF_)
     781#if defined (_CSOD_)
    782782    IndexElement** csod;
    783783    IndexElement*  temp;
     
    786786    int index1;
    787787    int num;
     788    int num_eq_op = 0;
     789    int offset;
    788790    locint entry;
    789791    int ii;
     
    836838    /* index domains */
    837839    locint** ind_dom;
    838 #if defined(_INDO_)
     840#if defined(_CSOD_)
    839841    locint** nlf;
     842    locint** fod;
    840843#endif
    841844#else
     
    10201023#else                                                                /* INDOPRO */
    10211024#if defined(_INDO_)
    1022     ind_dom = (locint **)  malloc(sizeof(locint*) * ADOLC_CURRENT_TAPE_INFOS.stats[NUM_MAX_LIVES]);
    1023 #if defined(_NLF_)
    1024     nlf = (locint **)  malloc(sizeof(locint*) * ADOLC_CURRENT_TAPE_INFOS.stats[NUM_MAX_LIVES]);
    1025     csod = (struct IndexElement**) malloc(sizeof(struct IndexElement*) * indcheck);
    1026     for(i=0;i<indcheck;i++)
    1027       {
    1028         csod[i] = (struct IndexElement*) malloc(sizeof(struct IndexElement));
    1029         csod[i]->entry = 0;
    1030         csod[i]->next = NULL;
    1031       } 
    1032 
    1033 #endif
    1034     for(i=0;i<ADOLC_CURRENT_TAPE_INFOS.stats[NUM_MAX_LIVES];i++)
     1025#if !defined(_CSOD_)
     1026        ind_dom = (locint **)  malloc(sizeof(locint*) * ADOLC_CURRENT_TAPE_INFOS.stats[NUM_MAX_LIVES]);
     1027        max_ind_dom = ADOLC_CURRENT_TAPE_INFOS.stats[NUM_MAX_LIVES];
     1028#else
     1029        ind_dom = (locint **)  malloc(sizeof(locint*) * (ADOLC_CURRENT_TAPE_INFOS.stats[NUM_MAX_LIVES]+ADOLC_CURRENT_TAPE_INFOS.stats[NUM_EQ_PROD]));
     1030        nlf = (locint **)  malloc(sizeof(locint*) * (ADOLC_CURRENT_TAPE_INFOS.stats[NUM_MAX_LIVES]+ADOLC_CURRENT_TAPE_INFOS.stats[NUM_EQ_PROD]));
     1031        fod = (locint **)  malloc(sizeof(locint*) * (ADOLC_CURRENT_TAPE_INFOS.stats[NUM_MAX_LIVES]+ADOLC_CURRENT_TAPE_INFOS.stats[NUM_EQ_PROD]));   
     1032        csod = (struct IndexElement**) malloc(sizeof(struct IndexElement*) * indcheck);
     1033        for(i=0;i<indcheck;i++)
     1034          {
     1035            csod[i] = (struct IndexElement*) malloc(sizeof(struct IndexElement));
     1036            csod[i]->entry = 0;
     1037            csod[i]->next = NULL;
     1038          }
     1039        max_ind_dom = ADOLC_CURRENT_TAPE_INFOS.stats[NUM_MAX_LIVES]+ADOLC_CURRENT_TAPE_INFOS.stats[NUM_EQ_PROD];
     1040        offset = ADOLC_CURRENT_TAPE_INFOS.stats[NUM_MAX_LIVES];
     1041#endif
     1042    for(i=0;i<max_ind_dom;i++)
    10351043    {
    1036         ind_dom[i] = (locint *)  malloc(sizeof(locint) * (NUMNNZ+2));
    1037         ind_dom[i][0] = 0;
    1038         ind_dom[i][1] = NUMNNZ;
    1039 #if defined(_NLF_)
    1040         nlf[i] = (locint *)  malloc(sizeof(locint) * (NUMNNZ+2));
    1041         nlf[i][0] = 0;
    1042         nlf[i][1] = NUMNNZ;
     1044      // structure of ind_dom[i]:
     1045      // ind_dom[i][0] = number of entries
     1046      // ind_dom[i][1] = length of array-2  = place for entries
     1047      ind_dom[i] = (locint *)  malloc(sizeof(locint) * (NUMNNZ+2));
     1048      ind_dom[i][0] = 0;
     1049      ind_dom[i][1] = NUMNNZ;
     1050#if defined(_CSOD_)
     1051      // structure of nlf: same as ind_dom
     1052      nlf[i] = (locint *)  malloc(sizeof(locint) * (NUMNNZ+2));
     1053      nlf[i][0] = 0;
     1054      nlf[i][1] = NUMNNZ;
     1055      // structure of fod[i]
     1056      // fod[i][0] = number of ind_doms contained in fod[i]
     1057      // fod[i][1] = length of array-3  = place for entries
     1058      // fod[i][2] = number of entries
     1059      fod[i] = (locint *)  malloc(sizeof(locint) * (103));
     1060      fod[i][0] = 0;
     1061      fod[i][1] = 100;
     1062      fod[i][2] = 0;
    10431063#endif
    10441064    }
    1045 
    1046     max_ind_dom = ADOLC_CURRENT_TAPE_INFOS.stats[NUM_MAX_LIVES];
    10471065#if defined(_NONLIND_)
    10481066
     
    12801298#if defined(_INDO_)
    12811299                copy_index_domain(res, arg, ind_dom);
    1282 #if defined(_NLF_)
    1283                 copy_index_domain(res, arg, nlf);
     1300#if defined(_CSOD_)
     1301                for(i=0;i<fod[arg][2]+3;i++)
     1302                  {
     1303                    fod[res][i] = fod[arg][i];
     1304                  }
     1305                printf("fod \n");
     1306                for(i=0;i<fod[res][2]+3;i++)
     1307                  printf(" %d ",fod[res][i]);
     1308                printf("\n");
     1309
    12841310#endif
    12851311#else
     
    13081334#if defined(_INDO_)
    13091335                ind_dom[res][0]=0;
    1310 #if defined(_NLF_)
     1336#if defined(_CSOD_)
    13111337                nlf[res][0]=0;
    13121338#endif
     
    13341360#if defined(_INDO_)
    13351361                ind_dom[res][0]=0;
    1336 #if defined(_NLF_)
     1362#if defined(_CSOD_)
    13371363                nlf[res][0]=0;
    13381364#endif
     
    13601386#if defined(_INDO_)
    13611387                ind_dom[res][0]=0;
    1362 #if defined(_NLF_)
     1388#if defined(_CSOD_)
    13631389                nlf[res][0]=0;
    13641390#endif
     
    14311457                }
    14321458#endif
    1433 #if defined(_NLF_)
     1459#if defined(_CSOD_)
    14341460                               
    1435                 for(l=2;l<=nlf[res][0]+1;l++)
     1461                int startfod = 3;
     1462                for(l=0;l<fod[res][0];l++)
    14361463                  {
    1437                     int index = nlf[res][l];
    1438                     for(ii=2;ii<=ind_dom[index][0]+1;ii++)
     1464                    for(ii=1;ii<=fod[res][startfod];ii++)
    14391465                      {
    1440                         index1 = ind_dom[index][ii];
    1441                         for(jj=2;jj<=ind_dom[index][0]+1;jj++)
     1466                        index1 = fod[res][startfod+ii];
     1467                        for(jj=1;jj<=fod[res][startfod];jj++)
    14421468                          {
    14431469                            temp1 = csod[index1];
     
    14511477                                    temp = (struct IndexElement*) malloc(sizeof(struct IndexElement));
    14521478                                    temp1->next = temp;
    1453                                     temp->entry = ind_dom[index][jj];
     1479                                    temp->entry = fod[res][startfod+jj];
    14541480                                    temp->next = NULL;
    14551481                                    csod[index1]->entry++;
     
    14581484                                else
    14591485                                  {
    1460                                     if (temp->entry < ind_dom[index][jj])
     1486                                    if (temp->entry < fod[res][startfod+jj])
    14611487                                      {
    14621488                                        temp1 = temp;
     
    14651491                                    else
    14661492                                      {
    1467                                         if (temp->entry == ind_dom[index][jj])
     1493                                        if (temp->entry == fod[res][startfod+jj])
    14681494                                          i = 1;
    14691495                                        else
    14701496                                          {
    14711497                                            temp = (struct IndexElement*) malloc(sizeof(struct IndexElement));
    1472                                             temp->entry = ind_dom[index][jj];
     1498                                            temp->entry = fod[res][startfod+jj];
    14731499                                            csod[index1]->entry++;
    14741500                                            temp->next = temp1->next;
     
    14821508                          }
    14831509                      }
     1510                    startfod += (fod[res][startfod]+1);
    14841511                  }
    14851512#endif
     
    15341561#if defined(_INDO_)
    15351562                merge_2_index_domains(res, arg, ind_dom);
    1536 #if defined(_NLF_)
     1563#if defined(_CSOD_)
    15371564                nlf[res][0] = 1;
    15381565                nlf[res][2] = res;
     
    15811608#if defined(_INDO_)
    15821609                merge_2_index_domains(res, arg, ind_dom);
    1583 #if defined(_NLF_)
     1610#if defined(_CSOD_)
    15841611                merge_2_index_domains(res, arg, nlf);
    15851612#endif
     
    16361663#if defined(_INDO_)
    16371664                merge_2_index_domains(res, arg, ind_dom);
    1638 #if defined(_NLF_)
     1665#if defined(_CSOD_)
    16391666                nlf[res][0] = 1;
    16401667                nlf[res][2] = res;
     
    17121739
    17131740#if !defined(_NTIGHT_)
    1714                 dp_T0[res] = dp_T0[arg1] +
    1715                                                dp_T0[arg2];
     1741                dp_T0[res] = dp_T0[arg1] + dp_T0[arg2];
    17161742#endif /* !_NTIGHT_ */
    17171743
    17181744#if defined(_INDO_)
    17191745                combine_2_index_domains(res, arg1, arg2, ind_dom);
    1720 #if defined(_NLF_)
    1721                 combine_2_index_domains(res, arg1, arg2, nlf);
     1746#if defined(_CSOD_)
     1747                for(i=0;i<fod[arg1][2]+3;i++)
     1748                  {
     1749                    fod[res][i] = fod[arg1][i];
     1750                  }
     1751                index = fod[res][2]+3;
     1752                for(i=3;i<fod[arg2][2]+3;i++)
     1753                  {
     1754                    fod[res][index++] = fod[arg2][i];
     1755                  }
     1756                fod[res][2]+=fod[arg2][2];
     1757                fod[res][0]+=fod[arg2][0];
     1758
     1759                printf("fod \n");
     1760                for(i=0;i<fod[res][2]+3;i++)
     1761                  printf(" %d ",fod[res][i]);
     1762                printf("\n");
    17221763#endif
    17231764#else
     
    17531794#if defined(_INDO_)
    17541795                copy_index_domain(res, arg, ind_dom);
    1755 #if defined(_NLF_)
    1756                 copy_index_domain(res, arg, nlf);
     1796#if defined(_CSOD_)
     1797                for(i=0;i<fod[arg][2]+3;i++)
     1798                  {
     1799                    fod[res][i] = fod[arg][i];
     1800                  }
     1801                printf("fod \n");
     1802                for(i=0;i<fod[res][2]+3;i++)
     1803                  printf(" %d ",fod[res][i]);
     1804                printf("\n");
    17571805#endif
    17581806#else
     
    17841832#if defined(_INDO_)
    17851833                combine_2_index_domains(res, arg1, arg2, ind_dom);
    1786 #if defined(_NLF_)
     1834#if defined(_CSOD_)
    17871835                combine_2_index_domains(res, arg1, arg2, nlf);
    17881836#endif
     
    18191867#if defined(_INDO_)
    18201868                copy_index_domain(res, arg, ind_dom);
    1821 #if defined(_NLF_)
    1822 #endif
     1869#if defined(_CSOD_)
    18231870                copy_index_domain(res, arg, nlf);
     1871#endif
    18241872#else
    18251873#if !defined(_ZOS_) /* BREAK_ZOS */
     
    18481896#if defined(_INDO_)
    18491897                combine_2_index_domains(res, arg1, arg2, ind_dom);
    1850 #if defined(_NLF_)
    1851                 nlf[res][0] = 1;
    1852                 nlf[res][2] = res;
     1898#if defined(_CSOD_)
     1899                fod[res][0] = 0;
     1900                index=3;
     1901                if (ind_dom[res][0] > 0)
     1902                  {
     1903                    fod[res][0]++;
     1904                    fod[res][index++] = ind_dom[res][0];
     1905                    for(l=2;l<ind_dom[res][0]+2;l++)
     1906                      fod[res][index++]=ind_dom[res][l];
     1907                    fod[res][2]+=ind_dom[res][0]+1;
     1908                   
     1909                  }
     1910                printf("fod \n");
     1911                for(i=0;i<fod[res][2]+3;i++)
     1912                  printf(" %d ",fod[res][i]);
     1913                printf("\n");
     1914
    18531915#endif
    18541916#if defined(_NONLIND_)
     
    19041966#if defined(_INDO_)
    19051967                merge_3_index_domains(res, arg1, arg2, ind_dom);
    1906 #if defined(_NLF_)
    1907                 /* copy_index_domain(res, arg1, nlf); */
     1968#if defined(_CSOD_)
     1969                // take
     1970                combine_2_index_domains(offset, arg1, arg2, ind_dom);
     1971                index = fod[res][2]+3;
     1972
     1973                if (ind_dom[offset][0] > 0)
     1974                  {
     1975                    fod[res][0]++;
     1976                    fod[res][index++] = ind_dom[offset][0];
     1977                    for(l=2;l<ind_dom[offset][0]+2;l++)
     1978                      fod[res][index++]=ind_dom[offset][l];
     1979                    fod[res][2]+=ind_dom[offset][0]+1;
     1980                   
     1981                  }
     1982                printf("fod \n");
     1983                for(i=0;i<fod[res][2]+3;i++)
     1984                  printf(" %d ",fod[res][i]);
     1985                printf("\n");
    19081986#endif
    19091987#if defined(_NONLIND_)
     
    19582036#if defined(_INDO_)
    19592037                merge_3_index_domains(res, arg1, arg2, ind_dom);
    1960 #if defined(_NLF_)
     2038#if defined(_CSOD_)
    19612039                /* copy_index_domain(res, arg1, nlf); */
    19622040#endif
     
    20182096#if defined(_INDO_)
    20192097                copy_index_domain(res, arg, ind_dom);
    2020 #if defined(_NLF_)
    2021                 copy_index_domain(res, arg, nlf);
     2098#if defined(_CSOD_)
     2099
     2100                for(i=0;i<fod[arg][2]+3;i++)
     2101                  {
     2102                    fod[res][i] = fod[arg][i];
     2103                  }
     2104                printf("fod \n");
     2105                for(i=0;i<fod[res][2]+3;i++)
     2106                  printf(" %d ",fod[res][i]);
     2107                printf("\n");
    20222108#endif
    20232109#else
     
    20572143#if defined(_INDO_)
    20582144                combine_2_index_domains(res, arg1, arg2, ind_dom);
    2059 #if defined(_NLF_)
     2145#if defined(_CSOD_)
    20602146                nlf[res][0] = 1;
    20612147                nlf[res][2] = res;
     
    21222208#if defined(_INDO_)
    21232209                copy_index_domain(res, arg, ind_dom);
    2124 #if defined(_NLF_)
     2210#if defined(_CSOD_)
    21252211                nlf[res][0] = 1;
    21262212                nlf[res][2] = res;
     
    21782264#if defined(_INDO_)
    21792265                copy_index_domain(res, arg, ind_dom);
    2180 #if defined(_NLF_)
     2266#if defined(_CSOD_)
    21812267                copy_index_domain(res, arg, nlf);
    21822268#endif
     
    22052291#if defined(_INDO_)
    22062292                copy_index_domain(res, arg, ind_dom);
    2207 #if defined(_NLF_)
     2293#if defined(_CSOD_)
    22082294                copy_index_domain(res, arg, nlf);
    22092295#endif
     
    22432329#if defined(_INDO_)
    22442330                copy_index_domain(res, arg, ind_dom);
    2245 #if defined(_NLF_)
     2331#if defined(_CSOD_)
    22462332                nlf[res][0] = 1;
    22472333                nlf[res][2] = res;
     
    23032389#if defined(_INDO_)
    23042390                copy_index_domain(res, arg1, ind_dom);
    2305 #if defined(_NLF_)
     2391#if defined(_CSOD_)
    23062392                nlf[res][0] = 1;
    23072393                nlf[res][2] = res;
     
    23732459#if defined(_INDO_)
    23742460                copy_index_domain(res, arg1, ind_dom);
    2375 #if defined(_NLF_)
     2461#if defined(_CSOD_)
    23762462                nlf[res][0] = 1;
    23772463                nlf[res][2] = res;
     
    24402526#if defined(_INDO_)
    24412527                copy_index_domain(res, arg1, ind_dom);
    2442 #if defined(_NLF_)
     2528#if defined(_CSOD_)
    24432529                nlf[res][0] = 1;
    24442530                nlf[res][2] = res;
     
    24992585#if defined(_INDO_)
    25002586                copy_index_domain(res, arg1, ind_dom);
    2501 #if defined(_NLF_)
     2587#if defined(_CSOD_)
    25022588                nlf[res][0] = 1;
    25032589                nlf[res][2] = res;
     
    25962682#if defined(_INDO_)
    25972683                copy_index_domain(res, arg1, ind_dom);
    2598 #if defined(_NLF_)
     2684#if defined(_CSOD_)
    25992685                nlf[res][0] = 1;
    26002686                nlf[res][2] = res;
     
    26952781#if defined(_INDO_)
    26962782                copy_index_domain(res, arg1, ind_dom);
    2697 #if defined(_NLF_)
     2783#if defined(_CSOD_)
    26982784                nlf[res][0] = 1;
    26992785                nlf[res][2] = res;
     
    27542840#if defined(_INDO_)
    27552841                copy_index_domain(res, arg1, ind_dom);
    2756 #if defined(_NLF_)
     2842#if defined(_CSOD_)
    27572843                nlf[res][0] = 1;
    27582844                nlf[res][2] = res;
     
    28352921#if defined(_INDO_)
    28362922                copy_index_domain(res, arg1, ind_dom);
    2837 #if defined(_NLF_)
     2923#if defined(_CSOD_)
    28382924                nlf[res][0] = 1;
    28392925                nlf[res][2] = res;
     
    29373023#if defined(_INDO_)
    29383024                copy_index_domain(res, arg1, ind_dom);
    2939 #if defined(_NLF_)
     3025#if defined(_CSOD_)
    29403026                nlf[res][0] = 1;
    29413027                nlf[res][2] = res;
     
    29883074#if defined(_INDO_)
    29893075                copy_index_domain(res, arg, ind_dom);
    2990 #if defined(_NLF_)
     3076#if defined(_CSOD_)
    29913077                nlf[res][0] = 1;
    29923078                nlf[res][2] = res;
     
    30703156#if defined(_INDO_)
    30713157                copy_index_domain(res, arg, ind_dom);
    3072 #if defined(_NLF_)
     3158#if defined(_CSOD_)
    30733159                nlf[res][0] = 1;
    30743160                nlf[res][2] = res;
     
    32093295#if defined(_INDO_)
    32103296                copy_index_domain(res, arg, ind_dom);
    3211 #if defined(_NLF_)
     3297#if defined(_CSOD_)
    32123298                nlf[res][0] = 1;
    32133299                nlf[res][2] = res;
     
    33753461                    combine_2_index_domains(res, arg1, arg2, ind_dom);
    33763462#endif
    3377 #if defined(_NLF_)
     3463#if defined(_CSOD_)
    33783464#ifdef _TIGHT_
    33793465                    if (dp_T0[arg1] < dp_T0[arg2])
     
    35073593#if defined(_INDO_)
    35083594                copy_index_domain(res, arg, ind_dom);
    3509 #if defined(_NLF_)
     3595#if defined(_CSOD_)
    35103596                copy_index_domain(res, arg, nlf);
    35113597#endif
     
    35843670#if defined(_INDO_)
    35853671                copy_index_domain(res, arg, ind_dom);
    3586 #if defined(_NLF_)
     3672#if defined(_CSOD_)
    35873673                copy_index_domain(res, arg, nlf);
    35883674#endif
     
    36143700#if defined(_INDO_)
    36153701                copy_index_domain(res, arg, ind_dom);
    3616 #if defined(_NLF_)
     3702#if defined(_CSOD_)
    36173703                copy_index_domain(res, arg, nlf);
    36183704#endif
     
    36503736
    36513737                    combine_2_index_domains(res, arg1, arg2, ind_dom);
    3652 #if defined(_NLF_)
     3738#if defined(_CSOD_)
    36533739                    combine_2_index_domains(res, arg1, arg2, nlf);
    36543740#endif
    36553741#else
    36563742                    copy_index_domain(res, arg1, ind_dom);
    3657 #if defined(_NLF_)
     3743#if defined(_CSOD_)
    36583744                    copy_index_domain(res, arg1, nlf);
    36593745#endif
     
    36683754
    36693755                        combine_2_index_domains(res, arg1, arg2, ind_dom);
    3670 #if defined(_NLF_)
     3756#if defined(_CSOD_)
    36713757                    combine_2_index_domains(res, arg1, arg2, nlf);
    36723758#endif
     
    36743760#else
    36753761                        copy_index_domain(res, arg2, ind_dom);
    3676 #if defined(_NLF_)
     3762#if defined(_CSOD_)
    36773763                    copy_index_domain(res, arg1, nlf);
    36783764#endif
     
    37493835#if defined(_INDO_)
    37503836                    copy_index_domain(res, arg1, ind_dom);
    3751 #if defined(_NLF_)
     3837#if defined(_CSOD_)
    37523838#endif
    37533839#else
     
    39244010        } /* endswitch */
    39254011
     4012     printf("op %d \n",operation);
     4013
     4014#if defined (_CSOD_)
     4015      printf(" ind_dom: \n");
     4016      for(l=0;l<max_ind_dom;l++)
     4017        {
     4018          printf(" %d : ",l);
     4019          for(i=2;i<ind_dom[l][0]+2;i++)
     4020            printf(" %d ",ind_dom[l][i]);
     4021          printf("\n");
     4022        }
     4023      printf(" nlf: \n");
     4024      for(l=0;l<max_ind_dom;l++)
     4025        {
     4026          printf(" %d : ",l);
     4027          for(i=2;i<nlf[l][0]+2;i++)
     4028            printf(" %d ",nlf[l][i]);
     4029          printf("\n");
     4030        }
     4031#endif
     4032
    39264033        /* Read the next operation */
    39274034        operation=get_op_f();
     
    39724079      }
    39734080    free(ind_dom);
    3974 #if defined(_NLF_)
     4081#if defined(_CSOD_)
    39754082
    39764083    for( i=0; i < indcheck; i++) {
Note: See TracChangeset for help on using the changeset viewer.