Changeset 5634


Ignore:
Timestamp:
08/31/10 11:46:12 (15 years ago)
Author:
Mathieu Morlighem
Message:

Removed almost all gaussgrids

File:
1 edited

Legend:

Unmodified
Added
Removed

  • issm/trunk/src/c/objects/Elements/Tria.cpp

    r5633 r5634  
    14701470        int dim;
    14711471        const int numvertices=3;
    1472         double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    14731472        double  vx_values[numvertices];
    14741473        double  maxabsvx;
    … …  
    14781477
    14791478        /*retrive velocity values at nodes */
    1480         inputs->GetParameterValues(&vx_values[0],&gaussgrids[0][0],numvertices,VxEnum);
     1479        GetParameterListOnVertices(&vx_values[0],VxEnum);
    14811480
    14821481        /*process units if requested: */
    … …  
    14991498        int dim;
    15001499        const int numvertices=3;
    1501         double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    15021500        double  vy_values[numvertices];
    15031501        double  maxabsvy;
    … …  
    15071505
    15081506        /*retrive velocity values at nodes */
    1509         inputs->GetParameterValues(&vy_values[0],&gaussgrids[0][0],numvertices,VyEnum);
     1507        GetParameterListOnVertices(&vy_values[0],VyEnum);
    15101508
    15111509        /*process units if requested: */
    … …  
    15281526        int dim;
    15291527        const int numvertices=3;
    1530         double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    15311528        double  vz_values[numvertices];
    15321529        double  maxabsvz;
    … …  
    15361533
    15371534        /*retrive velocity values at nodes */
    1538         inputs->GetParameterValues(&vz_values[0],&gaussgrids[0][0],numvertices,VzEnum);
     1535        GetParameterListOnVertices(&vz_values[0],VzEnum);
    15391536
    15401537        /*process units if requested: */
    … …  
    15571554        int dim;
    15581555        const int numvertices=3;
    1559         double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    15601556        double  vel_values[numvertices];
    15611557        double  maxvel;
    … …  
    15651561
    15661562        /*retrive velocity values at nodes */
    1567         inputs->GetParameterValues(&vel_values[0],&gaussgrids[0][0],numvertices,VelEnum);
     1563        GetParameterListOnVertices(&vel_values[0],VelEnum);
    15681564
    15691565        /*process units if requested: */
    … …  
    15871583        int dim;
    15881584        const int numvertices=3;
    1589         double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    15901585        double  vx_values[numvertices];
    15911586        double  maxvx;
    … …  
    15951590
    15961591        /*retrive velocity values at nodes */
    1597         inputs->GetParameterValues(&vx_values[0],&gaussgrids[0][0],numvertices,VxEnum);
     1592        GetParameterListOnVertices(&vx_values[0],VxEnum);
    15981593
    15991594        /*process units if requested: */
    … …  
    16171612        int dim;
    16181613        const int numvertices=3;
    1619         double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    16201614        double  vy_values[numvertices];
    16211615        double  maxvy;
    … …  
    16251619
    16261620        /*retrive velocity values at nodes */
    1627         inputs->GetParameterValues(&vy_values[0],&gaussgrids[0][0],numvertices,VyEnum);
     1621        GetParameterListOnVertices(&vy_values[0],VyEnum);
    16281622
    16291623        /*process units if requested: */
    … …  
    16471641        int dim;
    16481642        const int numvertices=3;
    1649         double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    16501643        double  vz_values[numvertices];
    16511644        double  maxvz;
    … …  
    16551648
    16561649        /*retrive velocity values at nodes */
    1657         inputs->GetParameterValues(&vz_values[0],&gaussgrids[0][0],numvertices,VzEnum);
     1650        GetParameterListOnVertices(&vz_values[0],VzEnum);
    16581651
    16591652        /*process units if requested: */
    … …  
    16771670        int dim;
    16781671        const int numvertices=3;
    1679         double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    16801672        double  vel_values[numvertices];
    16811673        double  minvel;
    … …  
    16851677
    16861678        /*retrive velocity values at nodes */
    1687         inputs->GetParameterValues(&vel_values[0],&gaussgrids[0][0],numvertices,VelEnum);
     1679        GetParameterListOnVertices(&vel_values[0],VelEnum);
    16881680
    16891681        /*process units if requested: */
    … …  
    17071699        int dim;
    17081700        const int numvertices=3;
    1709         double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    17101701        double  vx_values[numvertices];
    17111702        double  minvx;
    … …  
    17151706
    17161707        /*retrive velocity values at nodes */
    1717         inputs->GetParameterValues(&vx_values[0],&gaussgrids[0][0],numvertices,VxEnum);
     1708        GetParameterListOnVertices(&vx_values[0],VxEnum);
    17181709
    17191710        /*process units if requested: */
    … …  
    17371728        int dim;
    17381729        const int numvertices=3;
    1739         double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    17401730        double  vy_values[numvertices];
    17411731        double  minvy;
    … …  
    17451735
    17461736        /*retrive velocity values at nodes */
    1747         inputs->GetParameterValues(&vy_values[0],&gaussgrids[0][0],numvertices,VyEnum);
     1737        GetParameterListOnVertices(&vy_values[0],VyEnum);
    17481738
    17491739        /*process units if requested: */
    … …  
    17671757        int dim;
    17681758        const int numvertices=3;
    1769         double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    17701759        double  vz_values[numvertices];
    17711760        double  minvz;
    … …  
    17751764
    17761765        /*retrive velocity values at nodes */
    1777         inputs->GetParameterValues(&vz_values[0],&gaussgrids[0][0],numvertices,VzEnum);
     1766        GetParameterListOnVertices(&vz_values[0],VzEnum);
    17781767
    17791768        /*process units if requested: */
    … …  
    18131802        double  gauss_weight;
    18141803        double  gauss_l1l2l3[3];
    1815         double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    18161804
    18171805        /* parameters: */
    … …  
    19041892        double  gauss_weight;
    19051893        double  gauss_l1l2l3[3];
    1906         double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    19071894
    19081895        /* parameters: */
    … …  
    19281915
    19291916        /* Recover input data: */
    1930         inputs->GetParameterValues(&obs_vx_list[0],&gaussgrids[0][0],3,VxObsEnum);
    1931         inputs->GetParameterValues(&obs_vy_list[0],&gaussgrids[0][0],3,VyObsEnum);
    1932         inputs->GetParameterValues(&vx_list[0],&gaussgrids[0][0],3,VxEnum);
    1933         inputs->GetParameterValues(&vy_list[0],&gaussgrids[0][0],3,VyEnum);
    1934         inputs->GetParameterValues(&weights_list[0],&gaussgrids[0][0],3,WeightsEnum);
     1917        GetParameterListOnVertices(&obs_vx_list[0],VxObsEnum);
     1918        GetParameterListOnVertices(&obs_vy_list[0],VyObsEnum);
     1919        GetParameterListOnVertices(&vx_list[0],VxEnum);
     1920        GetParameterListOnVertices(&vy_list[0],VyEnum);
     1921        GetParameterListOnVertices(&weights_list[0],WeightsEnum);
    19351922
    19361923        /*retrieve some parameters: */
    … …  
    20272014        double  gauss_weight;
    20282015        double  gauss_l1l2l3[3];
    2029         double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    20302016
    20312017        /* parameters: */
    … …  
    20552041
    20562042        /* Recover input data: */
    2057         inputs->GetParameterValues(&obs_vx_list[0],&gaussgrids[0][0],3,VxObsEnum);
    2058         inputs->GetParameterValues(&obs_vy_list[0],&gaussgrids[0][0],3,VyObsEnum);
    2059         inputs->GetParameterValues(&vx_list[0],&gaussgrids[0][0],3,VxEnum);
    2060         inputs->GetParameterValues(&vy_list[0],&gaussgrids[0][0],3,VyEnum);
    2061         inputs->GetParameterValues(&weights_list[0],&gaussgrids[0][0],3,WeightsEnum);
     2043        GetParameterListOnVertices(&obs_vx_list[0],VxObsEnum);
     2044        GetParameterListOnVertices(&obs_vy_list[0],VyObsEnum);
     2045        GetParameterListOnVertices(&vx_list[0],VxEnum);
     2046        GetParameterListOnVertices(&vy_list[0],VyEnum);
     2047        GetParameterListOnVertices(&weights_list[0],WeightsEnum);
    20622048
    20632049        /*retrieve some parameters: */
    … …  
    21592145        double  gauss_weight;
    21602146        double  gauss_l1l2l3[3];
    2161         double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    21622147
    21632148        /* parameters: */
    … …  
    21872172
    21882173        /* Recover input data: */
    2189         inputs->GetParameterValues(&obs_vx_list[0],&gaussgrids[0][0],3,VxObsEnum);
    2190         inputs->GetParameterValues(&obs_vy_list[0],&gaussgrids[0][0],3,VyObsEnum);
    2191         inputs->GetParameterValues(&vx_list[0],&gaussgrids[0][0],3,VxEnum);
    2192         inputs->GetParameterValues(&vy_list[0],&gaussgrids[0][0],3,VyEnum);
    2193         inputs->GetParameterValues(&weights_list[0],&gaussgrids[0][0],3,WeightsEnum);
     2174        GetParameterListOnVertices(&obs_vx_list[0],VxObsEnum);
     2175        GetParameterListOnVertices(&obs_vy_list[0],VyObsEnum);
     2176        GetParameterListOnVertices(&vx_list[0],VxEnum);
     2177        GetParameterListOnVertices(&vy_list[0],VyEnum);
     2178        GetParameterListOnVertices(&weights_list[0],WeightsEnum);
    21942179
    21952180        /*retrieve some parameters: */
    … …  
    22892274        double  gauss_weight;
    22902275        double  gauss_l1l2l3[3];
    2291         double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    22922276
    22932277        /* parameters: */
    … …  
    23192303
    23202304        /* Recover input data: */
    2321         inputs->GetParameterValues(&obs_vx_list[0],&gaussgrids[0][0],3,VxObsEnum);
    2322         inputs->GetParameterValues(&obs_vy_list[0],&gaussgrids[0][0],3,VyObsEnum);
    2323         inputs->GetParameterValues(&vx_list[0],&gaussgrids[0][0],3,VxEnum);
    2324         inputs->GetParameterValues(&vy_list[0],&gaussgrids[0][0],3,VyEnum);
    2325         inputs->GetParameterValues(&weights_list[0],&gaussgrids[0][0],3,WeightsEnum);
     2305        GetParameterListOnVertices(&obs_vx_list[0],VxObsEnum);
     2306        GetParameterListOnVertices(&obs_vy_list[0],VyObsEnum);
     2307        GetParameterListOnVertices(&vx_list[0],VxEnum);
     2308        GetParameterListOnVertices(&vy_list[0],VyEnum);
     2309        GetParameterListOnVertices(&weights_list[0],WeightsEnum);
    23262310
    23272311        /*retrieve some parameters: */
    … …  
    24212405        double  gauss_weight;
    24222406        double  gauss_l1l2l3[3];
    2423         double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    24242407
    24252408        /* parameters: */
    … …  
    24522435        /* Recover input data: */
    24532436        inputs->GetParameterValue(&S,SurfaceAreaEnum);
    2454         inputs->GetParameterValues(&obs_vx_list[0],&gaussgrids[0][0],3,VxObsEnum);
    2455         inputs->GetParameterValues(&obs_vy_list[0],&gaussgrids[0][0],3,VyObsEnum);
    2456         inputs->GetParameterValues(&vx_list[0],&gaussgrids[0][0],3,VxEnum);
    2457         inputs->GetParameterValues(&vy_list[0],&gaussgrids[0][0],3,VyEnum);
    2458         inputs->GetParameterValues(&weights_list[0],&gaussgrids[0][0],3,WeightsEnum);
     2437        GetParameterListOnVertices(&obs_vx_list[0],VxObsEnum);
     2438        GetParameterListOnVertices(&obs_vy_list[0],VyObsEnum);
     2439        GetParameterListOnVertices(&vx_list[0],VxEnum);
     2440        GetParameterListOnVertices(&vy_list[0],VyEnum);
     2441        GetParameterListOnVertices(&weights_list[0],WeightsEnum);
    24592442
    24602443        /*retrieve some parameters: */
    … …  
    31373120        double       xyz_list[numvertices][3];
    31383121        int*         doflist=NULL;
    3139         double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    31403122
    31413123        /* gaussian points: */
    … …  
    32013183
    32023184        /*Modify z so that it reflects the surface*/
    3203         surface_input->GetParameterValues(&surface_list[0],&gaussgrids[0][0],3);
     3185        GetParameterListOnVertices(&surface_list[0],SurfaceEnum);
    32043186        for(i=0;i<numvertices;i++) xyz_list[i][2]=surface_list[i];
    32053187
    … …  
    51395121        double       xyz_list[numvertices][3];
    51405122        int*         doflist=NULL;
    5141         double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    51425123
    51435124        /* grid data: */
    … …  
    51905171
    51915172        /* Recover input data: */
    5192         inputs->GetParameterValues(&obs_vx_list[0],&gaussgrids[0][0],3,VxObsEnum);
    5193         inputs->GetParameterValues(&obs_vy_list[0],&gaussgrids[0][0],3,VyObsEnum);
    5194 
    5195         inputs->GetParameterValues(&vx_list[0],&gaussgrids[0][0],3,VxEnum);
    5196         inputs->GetParameterValues(&vy_list[0],&gaussgrids[0][0],3,VyEnum);
    5197 
    5198         inputs->GetParameterValues(&weights_list[0],&gaussgrids[0][0],3,WeightsEnum);
     5173        GetParameterListOnVertices(&obs_vx_list[0],VxObsEnum);
     5174        GetParameterListOnVertices(&obs_vy_list[0],VyObsEnum);
     5175        GetParameterListOnVertices(&vx_list[0],VxEnum);
     5176        GetParameterListOnVertices(&vy_list[0],VyEnum);
     5177        GetParameterListOnVertices(&weights_list[0],WeightsEnum);
    51995178
    52005179        inputs->GetParameterValue(&response,CmResponseEnum);
    … …  
    53785357        double       xyz_list[numvertices][3];
    53795358        int*         doflist=NULL;
    5380         double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    53815359
    53825360        /* grid data: */
    … …  
    54295407
    54305408        /* Recover input data: */
    5431         inputs->GetParameterValues(&obs_vx_list[0],&gaussgrids[0][0],3,VxObsEnum);
    5432         inputs->GetParameterValues(&obs_vy_list[0],&gaussgrids[0][0],3,VyObsEnum);
    5433 
    5434         inputs->GetParameterValues(&vx_list[0],&gaussgrids[0][0],3,VxEnum);
    5435         inputs->GetParameterValues(&vy_list[0],&gaussgrids[0][0],3,VyEnum);
    5436 
    5437         inputs->GetParameterValues(&weights_list[0],&gaussgrids[0][0],3,WeightsEnum);
     5409        GetParameterListOnVertices(&obs_vx_list[0],VxObsEnum);
     5410        GetParameterListOnVertices(&obs_vy_list[0],VyObsEnum);
     5411        GetParameterListOnVertices(&vx_list[0],VxEnum);
     5412        GetParameterListOnVertices(&vy_list[0],VyEnum);
     5413        GetParameterListOnVertices(&weights_list[0],WeightsEnum);
    54385414
    54395415        inputs->GetParameterValue(&response,CmResponseEnum);
    … …  
    60716047        double alpha2,vx,vy;
    60726048        double geothermalflux_value;
    6073         double    alpha2_list[numvertices];                                       //TO BE DELETED
    6074         double    gauss[numvertices][numvertices] = {{1,0,0},{0,1,0},{0,0,1}}; //TO BE DELETED
     6049        double alpha2_list[numvertices];                                    //TO BE DELETED
     6050        double gauss[numvertices][numvertices] = {{1,0,0},{0,1,0},{0,0,1}}; //TO BE DELETED
    60756051        double vx_list[numvertices]; //TO BE DELETED
    60766052        double vy_list[numvertices]; //TO BE DELETED
    … …  
    66536629        double  gauss_weight;
    66546630        double  gauss_l1l2l3[3];
    6655         double  gaussgrids[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    66566631
    66576632        /* parameters: */
    … …  
    69316906        double       thickness[numvertices];
    69326907        double       rho_ice,g;
    6933         double       gauss[numvertices][numvertices]={{1,0,0},{0,1,0},{0,0,1}};
    69346908        Input*       vz_input=NULL;
    69356909        double*      vz_ptr=NULL;
    … …  
    69706944        rho_ice=matpar->GetRhoIce();
    69716945        g=matpar->GetG();
    6972         inputs->GetParameterValues(&thickness[0],&gauss[0][0],3,ThicknessEnum);
     6946        GetParameterListOnVertices(&thickness[0],ThicknessEnum);
    69736947       
    69746948        for(i=0;i<numvertices;i++){
Note: See TracChangeset for help on using the changeset viewer.