Changes

#172 (Nov 12, 2022, 7:17:37 AM)

  1. Add print statements for debugging — dongjoon / detail
  2. add more debugOutput statements to show timing in SLCP — dongjoon / detail
  3. Fix skipping geometry Design_Variable in TACS — galbramc / detail
  4. More attempts at fixing the exodus makefile — galbramc / detail
  5. Stanalizer fix — galbramc / detail
  6. Finally got netcdf.h check correct on both macOS and Linux — galbramc / detail
  7. Another sanitizer fix — galbramc / detail
  8. sanitizer fixes — galbramc / detail
  9. Check for exodusWriter when running plato examples — galbramc / detail
  10. Write out surface to volume index map for plato — galbramc / detail
  11. Lint fix for mses — galbramc / detail
  12. Fix output redirect for detecing netcdf.h — galbramc / detail
  13. Update fun3d test to Mesh_Morph. Add first plato mesh generation examples. — galbramc / detail
  14. Allow tetgen to write MultiDomain mesh files. Updates to Exodus mesh writer. — galbramc / detail
  15. More checks when reading mapbc files in aimMesh — galbramc / detail
  16. Fun3D lint fixes — galbramc / detail
  17. fix - unused variables — ryan / detail
  18. fix - null catch in avlaim — ryan / detail
  19. fix fun3d morphing — ryan / detail
  20. Lint fixes — galbramc / detail
  21. Warning fix — galbramc / detail
  22. Fix last commit — galbramc / detail
  23. Fix waring and add missing files — galbramc / detail
  24. Support tetgen generated regions — galbramc / detail
  25. More lint fixes — galbramc / detail
  26. Fix lint warnings — galbramc / detail
  27. Update Makefile for plato. Add element group names for aflr3 and tetgen. — galbramc / detail
  28. Add plato skeleton AIM — galbramc / detail
  29. Initial cut at an exodus writer — galbramc / detail
  30. Fix for incorrect Design_Variable in masstranAIM — galbramc / detail
  31. Allow unconverged SmoothUV status to continue during parameterization — haimes / detail
  32. Fix masstran derivatives with multiple bodies — galbramc / detail
  33. Fix closeBound to write bound and vertexset files for journaling — haimes / detail
  34. fix typo in corsairlite slcp — dongjoon / detail
  35. delete redundant kulfan csm and finalize multifidelity scripts — dongjoon / detail
  36. Fix warnings from last commit — galbramc / detail
  37. small fix on Kulfan MSES wrapper and removing remnants of full corsair — dongjoon / detail
  38. Remove some depricated checks — galbramc / detail
  39. Journal the writing of files — haimes / detail
  40. Move up setting outLevel — galbramc / detail
  41. Use ocsmLoadFromModel for static geometry. Expost stepSize for finite differnce in pyCAPS. — galbramc / detail
  42. Add the ability to control the sensitivity finite-difference step size from CAPS proper — haimes / detail
  43. Remove old html doc before copying — galbramc / detail
  44. Update cart3dAIM dox — galbramc / detail
  45. Update pyCAPS dox dependencies — galbramc / detail
  46. update corsairlite files to python3.8 compatible — dongjoon / detail
  47. Use TIM viewer in pyCAPS view methods — galbramc / detail
  48. Allow ftype=None for FieldIn/FieldOut data sets — galbramc / detail
  49. fix bug in serveESP that did not properly update the build status message after File->Edit->Save; implement analytic sensitivities in udpWaffle; add designX* test cases; make sure userPointer gets reset after perturbation is removed — jfdannen / detail
  50. add ocsmAdjoint and testOcsmAdjoint (in serveESP); allow timPlotter to have two y-axes; add plotter3.py test case — jfdannen / detail
  51. update adjoint test to remove (incorrect) compiler warning — jfdannen / detail
  52. remove data/ide (because it is redundant with data/tutorial6); eliminate automatic build when entering CAPS mode (depend on the lazy evaluation in CAPS); enable sensitivity calculations for scribes and trims; add designW* test cases; fix bug that put wrong _edgeID on scribed Edges; update attrScribe, scribeWing, and scribedSphere test cases; create splitEdges.udc; add splitEdges* test cases — jfdannen / detail
  53. remember name of Pyscript file for next Tool->Pyscript; clean up ESP messages after builds; remove code in esp.timSetCaps that removed the _finalize method; fix memory problem in timViewer for very long scene-graph meta-data buffers; reset magic number to zero to catch cases when a MODL is tried to be used after ocsmFree; modify IMPORT so that it properly marks NODE_BODYs; force tessellation if needed in timViewer; fix bug in ocsmLoadFromEgads that did not properly initialize the .ileft and .irite Branches; remove underscore attributes in .egads files in ocsmLoadFromEgads; force python garbage collection at end of each pyscript execution; add interactive testScript — jfdannen / detail
  54. add PLOT_TESSSENS and PLOT_TEMPCOMPS flags to help in debugging tessellation sensitivities; remove tessellation velocity calculation shortcut to make Edge and Face tessellation velocities consistent; refactor signalError() and signalError2(); allow plotTypes =7 for x-component, =8 for y-component, =9 for z-component of sensitivity — jfdannen / detail
  55. Add simple figure for tri ordering with mixed quad faces — galbramc / detail
  56. Allow unconverged SmoothUV status to continue during parameterization — haimes / detail
  57. More general ubuntu suppression — galbramc / detail
  58. Ubuntu valgrid supression — galbramc / detail
  59. scipy 1.9.2. and 1.9.3 have valgrind issues — galbramc / detail
  60. Add seacas to macys — galbramc / detail
  61. Add SEACAS environment variable for exodus — galbramc / detail
  62. Upgrade cvxopt build for gcc 12 on viggen — galbramc / detail

#171 (Oct 8, 2022, 2:51:45 PM)

  1. update corsairlite slcp wrapper with speed-up updates — dongjoon / detail
  2. Move DESPMTR check into just TACS — galbramc / detail
  3. Don't write DESPMTR to bdf files — galbramc / detail
  4. Fix segFault from last commit, again — haimes / detail
  5. Fix segFault from last commit — haimes / detail
  6. First cut at using MD5 checksums to check API inputs for continuation mode — haimes / detail
  7. Fix masstran memory leak — galbramc / detail
  8. masstran AIM updates — galbramc / detail
  9. Another OCSM FD sensitivity update — haimes / detail
  10. The calling sequence (SetDtime) was wrong and should be after setting the parameter when calculating sensitivities — haimes / detail
  11. add some debugging output capabilities to corsairlite and properly pass problemObj in Kulfan.py — dongjoon / detail
  12. Compile fix for Windoze — galbramc / detail
  13. Fix continuation with a function that returned an error — galbramc / detail
  14. Fix memory leak — galbramc / detail
  15. Fix masstran units and warning — galbramc / detail
  16. Update aim_copyValue — galbramc / detail
  17. Add geometric sensitvities to masstran AIM — galbramc / detail
  18. Suppress 'could not create compact unwind on M1' for now. — galbramc / detail
  19. Allow for FD sensitivities within the AIM sensitivity functions — haimes / detail
  20. add torch tests for kulfan corsairlite examples — dongjoon / detail
  21. update problem.__init__ to work if running pyscript but not from capsMode — jfdannen / detail
  22. Simplify multifidelity examples. Fig unit bug in caps.py. Support pint dimensionless to caps.py. Support udunits import in corsairlite. — galbramc / detail
  23. Fix typo in comment — haimes / detail
  24. Fix multifidelity corsairlite examples — dongjoon / detail
  25. Another cart3d twist example fix — galbramc / detail
  26. Cart3D example fixes — galbramc / detail
  27. Support Cart3D preSpec inputs and bodies. — galbramc / detail
  28. Forgot the pickle — galbramc / detail
  29. One more makefile fix — galbramc / detail
  30. corsair lite makefile example fix — galbramc / detail
  31. Dissable failing corsair examples again — galbramc / detail
  32. Fix corsairlite execution script — galbramc / detail
  33. Fix uninitialized journaling — galbramc / detail
  34. Run all phasing examples for corsairlite — galbramc / detail
  35. corsair fix for Python 3.10. Minor pyCAPS warning fix. — galbramc / detail
  36. Fix scan-build warning — galbramc / detail
  37. Fix for aim_relPath. Fix for sheet bodies with aim_storeMeshRef. — galbramc / detail
  38. fix morph example — ryan / detail
  39. Fix last commit — galbramc / detail
  40. Fix pyCAPS getValue for string and pointers — galbramc / detail
  41. Fix destroy_aimStorage to resolve segfault — galbramc / detail
  42. lint fixes — galbramc / detail
  43. I think we have mighty morphing powers — galbramc / detail
  44. test/fix morph — ryan / detail
  45. Skip wake with morphing for now. — galbramc / detail
  46. Fixes for mighty morphing powers, still needs some more testing. — galbramc / detail
  47. fix missing file — ryan / detail
  48. fix missing file — ryan / detail
  49. // This software has been cleared for public release on 05 Nov 2020, case number 88ABW-2020-3462.

    #include <string.h>
    #include <stdio.h>
    #include <errno.h>
    #include <ctype.h>

    #include "egads.h"        // Bring in egads utilss
    #include "capsTypes.h"    // Bring in CAPS types
    #include "aimUtil.h"      // Bring in AIM utils
    #include "aimMesh.h"      // Bring in AIM meshing

    #include "miscUtils.h"    // Bring in misc. utility functions
    #include "meshUtils.h"    // Bring in meshing utility functions
    #include "cfdTypes.h"     // Bring in cfd specific types
    #include "cfdUtils.h"
    #include "tecplotUtils.h" // Bring in tecplot utility functions
    #include "fun3dUtils.h"   // Bring in fun3d utility header
    #include "fun3dInputs.h"
    #include "ugridWriter.h"


    #ifdef WIN32
    #define strcasecmp  stricmp
    #define strncasecmp _strnicmp
    #endif


    // Extract the FEPOINT Tecoplot data from a FUN3D Aero-Loads file (connectivity is ignored) - dataMatrix = [numVariable][numDataPoint]
    int fun3d_readAeroLoad(void *aimInfo, char *filename, int *numVariable, char **variableName[],
                           int *numDataPoint, double ***dataMatrix)
    {

        int status = CAPS_SUCCESS; // Function return
        int i, j; // Indexing

        size_t linecap = 0;

        char *line = NULL; // Temporary line holder
        char *tempStr = NULL, *tempStr2 = NULL; // Temporary strings
        int stringLen = 0; // Length of string holder

        FILE *fp = NULL; // File pointer

        // Open file
        fp = aim_fopen(aimInfo, filename, "r");
        if (fp == NULL) {
            AIM_ERROR(aimInfo, "Unable to open file: %s\n", filename);
            return CAPS_IOERR;
        }

        printf("Reading FUN3D AeroLoad File - %s!!!!!!\n", filename);

        *numVariable = 0;
        *numDataPoint = 0;
        // Loop through file line by line until we have determined how many variables and data points there are
        while (*numVariable == 0 || *numDataPoint == 0) {

            // Get line from file
            status = getline(&line, &linecap, fp);
            if ((status < 0) || (line == NULL)) break;

            // Get variable list if available in file line
            if (strncmp("variables=", line, strlen("variables=")) == 0) {

                // Pull out substring at first occurrence of "
                tempStr = strstr(line, "\"");

                // Create a temperory string of the variables in a the folling format - ["a","ae"]
                stringLen  = strlen(tempStr)-1 + 2;

                tempStr2 = (char *) EG_alloc((stringLen +1)*sizeof(char *));
                if (tempStr2 == NULL) {
                    fclose(fp);
                    if (line != NULL) EG_free(line);
                    return EGADS_MALLOC;
                }

                tempStr2[0] = '[';
                strncpy(tempStr2+1, tempStr, strlen(tempStr)-1);
                tempStr2[stringLen-1] = ']';
                tempStr2[stringLen] = '\0';

                // Sort string into an array of strings
                status =  string_toStringDynamicArray(tempStr2, numVariable, variableName);

                if (tempStr2 != NULL) EG_free(tempStr2);
                tempStr2 = NULL;

                if (status != CAPS_SUCCESS) goto cleanup;

                // Print out list of variables found in load file
                printf("Variables found in file %s:\n",filename);
                for (i = 0; i < *numVariable; i++) printf("Variable %d = %s\n", i, (*variableName)[i]);
            }

            // Get the number of data points in file if available in file line
            if (strncmp("zone t=", line, strlen("zone t=")) == 0) {

                // Pull out substring at first occurrence of i=
                tempStr = strstr(line, "i=");

                // Retrieve the i= value
                sscanf(&tempStr[2], "%d", numDataPoint);

                // Print out the number of data points found in load file
                printf("Number of data points = %d, in file %s\n", *numDataPoint, filename);
            }

        }

        if (*numVariable != 0 && *numDataPoint != 0) {

            // Allocate dataMatrix array
            AIM_FREE(*dataMatrix);

            AIM_ALLOC(*dataMatrix, (*numVariable), double *, aimInfo, status);
            for (i = 0; i < *numVariable; i++) (*dataMatrix)[i] = NULL;

            for (i = 0; i < *numVariable; i++) {
                AIM_ALLOC((*dataMatrix)[i], (*numDataPoint), double, aimInfo, status);
            }

            // Loop through the file and fill up the data matrix
            for (j = 0; j < *numDataPoint; j++) {
                for (i = 0; i < *numVariable; i++) {
                    fscanf(fp, "%lf", &(*dataMatrix)[i][j]);
                }
            }

            // Output the first row of the dataMatrix
            //for (i = 0; i < *numVariable; i++) printf("Variable %d - %.6f\n", i, (*dataMatrix)[i][0]);

        } else {

            printf("No data values extracted from file - %s",filename);
            status = CAPS_BADVALUE;
        }

    cleanup:
        if (fp != NULL) fclose(fp);

        if (status != CAPS_SUCCESS) {
          if (*dataMatrix != NULL) {
            for (j = 0; j < *numVariable; j++) {
              AIM_FREE((*dataMatrix)[j]);
            }
            AIM_FREE((*dataMatrix));
          }
        }

        EG_free(line);
        return status;
    }


    static int
    fun3d_read2DBinaryUgrid(void *aimInfo, FILE *fp, meshStruct *surfaceMesh)
    {
      int    status = CAPS_SUCCESS;

      int    numNode, numLine, numTriangle, numQuadrilateral;
      int    numTetrahedral, numPyramid, numPrism, numHexahedral;
      int    i, elementIndex, numPoint, bcID;
      double *coords = NULL;

      /* we get a binary UGRID file */
      status = fread(&numNode,          sizeof(int), 1, fp);
      if (status != 1) { status = CAPS_IOERR; AIM_STATUS(aimInfo, status); }
      status = fread(&numTriangle,      sizeof(int), 1, fp);
      if (status != 1) { status = CAPS_IOERR; AIM_STATUS(aimInfo, status); }
      status = fread(&numQuadrilateral, sizeof(int), 1, fp);
      if (status != 1) { status = CAPS_IOERR; AIM_STATUS(aimInfo, status); }
      status = fread(&numTetrahedral,   sizeof(int), 1, fp);
      if (status != 1) { status = CAPS_IOERR; AIM_STATUS(aimInfo, status); }
      status = fread(&numPyramid,       sizeof(int), 1, fp);
      if (status != 1) { status = CAPS_IOERR; AIM_STATUS(aimInfo, status); }
      status = fread(&numPrism,         sizeof(int), 1, fp);
      if (status != 1) { status = CAPS_IOERR; AIM_STATUS(aimInfo, status); }
      status = fread(&numHexahedral,    sizeof(int), 1, fp);
      if (status != 1) { status = CAPS_IOERR; AIM_STATUS(aimInfo, status); }

      if ( numTetrahedral   +
           numPyramid       +
           numPrism         +
           numHexahedral  > 0) {
        AIM_ERROR(aimInfo, "Expecting a 2D ugrid file!!!");
        status = CAPS_IOERR;
        goto cleanup;
      }

      /*
      printf("\n Header from UGRID file: %d  %d %d  %d %d %d %d\n", numNode,
             numTriangle, numQuadrilateral, numTetrahedral, numPyramid, numPrism,
             numHexahedral);
       */

      AIM_ALLOC(coords, 3*numNode, double, aimInfo, status);

      /* read all of the vertices */
      status = fread(coords, sizeof(double), 3*numNode, fp);
      if (status != 3*numNode) { status = CAPS_IOERR; AIM_STATUS(aimInfo, status); }

      surfaceMesh->analysisType = UnknownMeshAnalysis;

      // Set that this is a volume mesh
      surfaceMesh->meshType = Surface2DMesh;

      // Numbers
      surfaceMesh->numNode = numNode;
      surfaceMesh->numElement = numTriangle + numQuadrilateral;

      surfaceMesh->meshQuickRef.useStartIndex = (int) true;

      surfaceMesh->meshQuickRef.numTriangle      = numTriangle;
      surfaceMesh->meshQuickRef.numQuadrilateral = numQuadrilateral;

      // Nodes - allocate
      AIM_ALLOC(surfaceMesh->node, surfaceMesh->numNode, meshNodeStruct, aimInfo, status);

      // Initialize
      for (i = 0; i < surfaceMesh->numNode; i++) {
          status = initiate_meshNodeStruct(&surfaceMesh->node[i],
                                            surfaceMesh->analysisType);
          AIM_STATUS(aimInfo, status);
      }

      // Nodes - set
      for (i = 0; i < surfaceMesh->numNode; i++) {

        // Copy node data
        surfaceMesh->node[i].nodeID = i+1;

        surfaceMesh->node[i].xyz[0] = coords[3*i+0];
        surfaceMesh->node[i].xyz[1] = coords[3*i+1];
        surfaceMesh->node[i].xyz[2] = coords[3*i+2];
      }
      AIM_FREE(coords);

      // Elements - allocate
      AIM_ALLOC(surfaceMesh->element, surfaceMesh->numElement, meshElementStruct, aimInfo, status);

      // Initialize
      for (i = 0; i < surfaceMesh->numElement; i++ ) {
          status = initiate_meshElementStruct(&surfaceMesh->element[i],
                                               surfaceMesh->analysisType);
          AIM_STATUS(aimInfo, status);
      }

      // Start of element index
      elementIndex = 0;

      // Elements -Set triangles
      if (numTriangle > 0)
        surfaceMesh->meshQuickRef.startIndexTriangle = elementIndex;

      numPoint = mesh_numMeshConnectivity(Triangle);
      for (i = 0; i < numTriangle; i++) {

          surfaceMesh->element[elementIndex].elementType = Triangle;
          surfaceMesh->element[elementIndex].elementID   = elementIndex+1;

          status = mesh_allocMeshElementConnectivity(&surfaceMesh->element[elementIndex]);
          if (status != CAPS_SUCCESS) goto cleanup;

          // read the element connectivity
          status = fread(surfaceMesh->element[elementIndex].connectivity,
                         sizeof(int), numPoint, fp);
          if (status != numPoint) { status = CAPS_IOERR; AIM_STATUS(aimInfo, status); }

          elementIndex += 1;
      }

      // Elements -Set quadrilateral
      if (numQuadrilateral > 0)
        surfaceMesh->meshQuickRef.startIndexQuadrilateral = elementIndex;

      numPoint = mesh_numMeshConnectivity(Quadrilateral);
      for (i = 0; i < numQuadrilateral; i++) {

          surfaceMesh->element[elementIndex].elementType = Quadrilateral;
          surfaceMesh->element[elementIndex].elementID   = elementIndex+1;

          status = mesh_allocMeshElementConnectivity(&surfaceMesh->element[elementIndex]);
          if (status != CAPS_SUCCESS) goto cleanup;

          // read the element connectivity
          status = fread(surfaceMesh->element[elementIndex].connectivity,
                         sizeof(int), numPoint, fp);
          if (status != numPoint) { status = CAPS_IOERR; AIM_STATUS(aimInfo, status); }

          elementIndex += 1;
      }

      // skip face ID section of the file
      status = fseek(fp, (numTriangle + numQuadrilateral)*sizeof(int), SEEK_CUR);
      if (status != 0) { status = CAPS_IOERR; AIM_STATUS(aimInfo, status); }

      // Get the number of Line elements
      status = fread(&numLine, sizeof(int), 1, fp);
      if (status != 1) { status = CAPS_IOERR; AIM_STATUS(aimInfo, status); }

      // Elements - re-allocate with Line elements
      surfaceMesh->meshQuickRef.numLine = numLine;
      AIM_REALL(surfaceMesh->element, surfaceMesh->numElement+numLine, meshElementStruct, aimInfo, status);

      surfaceMesh->meshQuickRef.startIndexLine = elementIndex;

      // Initialize
      for (i = surfaceMesh->numElement; i < surfaceMesh->numElement+numLine; i++ ) {
          status = initiate_meshElementStruct(&surfaceMesh->element[i],
                                               surfaceMesh->analysisType);
          AIM_STATUS(aimInfo, status);
      }
      surfaceMesh->numElement += numLine;

      numPoint = mesh_numMeshConnectivity(Line);
      for (i = 0; i < numLine; i++) {

        surfaceMesh->element[elementIndex].elementType = Line;
        surfaceMesh->element[elementIndex].elementID   = elementIndex+1;

        status = mesh_allocMeshElementConnectivity(&surfaceMesh->element[elementIndex]);
        AIM_STATUS(aimInfo, status);

        // read the element connectivity
        status = fread(surfaceMesh->element[elementIndex].connectivity,
                       sizeof(int), numPoint, fp);
        if (status != numPoint) { status = CAPS_IOERR; AIM_STATUS(aimInfo, status); }
        status = fread(&bcID, sizeof(int), 1, fp);
        if (status != 1) { status = CAPS_IOERR; AIM_STATUS(aimInfo, status); }

        surfaceMesh->element[elementIndex].markerID = bcID;

        elementIndex += 1;
      }

      status = CAPS_SUCCESS;

      cleanup:
          if (status != CAPS_SUCCESS)
            printf("Premature exit in getUGRID status = %d\n", status);

          EG_free(coords); coords = NULL;

          return status;
    }


    // Create a 3D BC for FUN3D from a 2D mesh
    int fun3d_2DBC(void *aimInfo,
                   cfdBoundaryConditionStruct *bcProps)
    {
        int status; // Function return status

        int i; // Indexing

        int faceBCIndex = -1, extrusionBCIndex = -1;

        // Find the faceBCIndex for the symmetry plane
        for (i = 0; i < bcProps->numSurfaceProp-1; i++) {
          if (bcProps->surfaceProp[i].surfaceType == Symmetry) {
            faceBCIndex = bcProps->surfaceProp[i].bcID;
            break;
          }
        }

        if (faceBCIndex == -1) {
          // Add plane boundary condition
          AIM_REALL(bcProps->surfaceProp, bcProps->numSurfaceProp+1, cfdSurfaceStruct, aimInfo, status);
          bcProps->numSurfaceProp += 1;

          status = initiate_cfdSurfaceStruct(&bcProps->surfaceProp[bcProps->numSurfaceProp-1]);
          AIM_STATUS(aimInfo, status);

          bcProps->surfaceProp[bcProps->numSurfaceProp-1].surfaceType = Symmetry;
          bcProps->surfaceProp[bcProps->numSurfaceProp-1].symmetryPlane = 1;

          // Find largest index value for bcID and set it plus 1 to the new surfaceProp
          for (i = 0; i < bcProps->numSurfaceProp-1; i++) {
              if (bcProps->surfaceProp[i].bcID >= faceBCIndex) {
                faceBCIndex = bcProps->surfaceProp[i].bcID + 1;
              }
          }
          bcProps->surfaceProp[bcProps->numSurfaceProp-1].bcID = faceBCIndex;
        }


        // Add extruded plane boundary condition
        AIM_REALL(bcProps->surfaceProp, bcProps->numSurfaceProp+1, cfdSurfaceStruct, aimInfo, status);
        bcProps->numSurfaceProp += 1;

        status = initiate_cfdSurfaceStruct(&bcProps->surfaceProp[bcProps->numSurfaceProp-1]);
        AIM_STATUS(aimInfo, status);

        bcProps->surfaceProp[bcProps->numSurfaceProp-1].surfaceType = Symmetry;
        bcProps->surfaceProp[bcProps->numSurfaceProp-1].symmetryPlane = 2;

        // Find largest index value for bcID and set it plus 1 to the new surfaceProp
        for (i = 0; i < bcProps->numSurfaceProp-1; i++) {
            if (bcProps->surfaceProp[i].bcID >= extrusionBCIndex) {
              extrusionBCIndex = bcProps->surfaceProp[i].bcID + 1;
            }
        }
        bcProps->surfaceProp[bcProps->numSurfaceProp-1].bcID = extrusionBCIndex;

        status = CAPS_SUCCESS;

    cleanup:
        return status;
    }


    // Create a 3D mesh for FUN3D from a 2D mesh
    int fun3d_2DMesh(void *aimInfo,
                     aimMeshRef *meshRef,
                     const char *projectName,
                     const mapAttrToIndexStruct *groupMap)
    {

        int status; // Function return status

        int i; // Indexing

        int faceBCIndex = -1, extrusionBCIndex = -1;

        double extrusion = -1.0; // Extrusion length

        // Flip coordinates
        int xMeshConstant = (int) true, yMeshConstant = (int) true, zMeshConstant= (int) true; // 2D mesh checks
        double tempCoord;

        meshStruct surfaceMesh;
        meshStruct volumeMesh;

        char filename[PATH_MAX];
    //    int elementIndex;
        FILE *fp = NULL;

        status = initiate_meshStruct(&surfaceMesh);
        AIM_STATUS(aimInfo, status);

        status = initiate_meshStruct(&volumeMesh);
        AIM_STATUS(aimInfo, status);

        sprintf(filename, "%s%s", meshRef->fileName, MESHEXTENSION);

        fp = fopen(filename, "rb");
        if (fp == NULL) {
          AIM_ERROR(aimInfo, "Cannot open file: %s\n", filename);
          status = CAPS_IOERR;
          goto cleanup;
        }

        status = fun3d_read2DBinaryUgrid(aimInfo, fp, &surfaceMesh);
        AIM_STATUS(aimInfo, status);

        // add boundary elements if they are missing
        if (surfaceMesh.meshQuickRef.numLine == 0) {
            status = mesh_addTess2Dbc(aimInfo, &surfaceMesh, groupMap);
            if (status != CAPS_SUCCESS) goto cleanup;
        }

        // Set the symmetry index for all Tri/Quad
        for (i = 0; i < surfaceMesh.numElement; i++) {

            if (surfaceMesh.element[i].elementType != Triangle &&
                surfaceMesh.element[i].elementType != Quadrilateral) {
                continue;
            }

            surfaceMesh.element[i].markerID = faceBCIndex;
        }

    #ifdef I_DONT_THINK_WE_NEED_THIS
        // Determine a suitable boundary index of the extruded plane
        *extrusionBCIndex = faceBCIndex;
        for (i = 0; i < surfaceMesh.meshQuickRef.numLine; i++) {

            if (surfaceMesh.meshQuickRef.startIndexLine >= 0) {
                elementIndex = surfaceMesh.meshQuickRef.startIndexLine + i;
            } else {
                elementIndex = surfaceMesh.meshQuickRef.listIndexLine[i];
            }

            marker = surfaceMesh.element[elementIndex].markerID;

            if (marker > *extrusionBCIndex)  {
                *extrusionBCIndex = marker;
            }
        }
        *extrusionBCIndex += 1;
    #endif

        // Check to make sure the face is on the y = 0 plane
        for (i = 0; i < surfaceMesh.numNode; i++) {

            if (surfaceMesh.node[i].xyz[1] != 0.0) {
                printf("\nSurface mesh is not on y = 0.0 plane, FUN3D could fail during execution for this 2D mesh!!!\n");
                break;
            }
        }

        // Constant x?
        for (i = 0; i < surfaceMesh.numNode; i++) {
            if ((surfaceMesh.node[i].xyz[0] - surfaceMesh.node[0].xyz[0]) > 1E-7) {
                xMeshConstant = (int) false;
                break;
            }
        }

        // Constant y?
        for (i = 0; i < surfaceMesh.numNode; i++) {
            if ((surfaceMesh.node[i].xyz[1] - surfaceMesh.node[0].xyz[1] ) > 1E-7) {
                yMeshConstant = (int) false;
                break;
            }
        }

        // Constant z?
        for (i = 0; i < surfaceMesh.numNode; i++) {
            if ((surfaceMesh.node[i].xyz[2] - surfaceMesh.node[0].xyz[2]) > 1E-7) {
                zMeshConstant = (int) false;
                break;
            }
        }

        if (yMeshConstant != (int) true) {
            printf("FUN3D expects 2D meshes be in the x-z plane... attempting to rotate mesh!\n");

            if (xMeshConstant == (int) true && zMeshConstant == (int) false) {
                printf("Swapping y and x coordinates!\n");
                for (i = 0; i < surfaceMesh.numNode; i++) {
                    tempCoord = surfaceMesh.node[i].xyz[0];
                    surfaceMesh.node[i].xyz[0] = surfaceMesh.node[i].xyz[1];
                    surfaceMesh.node[i].xyz[1] = tempCoord;
                }

            } else if(xMeshConstant == (int) false && zMeshConstant == (int) true) {

                printf("Swapping y and z coordinates!\n");
                for (i = 0; i < surfaceMesh.numNode; i++) {
                    tempCoord = surfaceMesh.node[i].xyz[2];
                    surfaceMesh.node[i].xyz[2] = surfaceMesh.node[i].xyz[1];
                    surfaceMesh.node[i].xyz[1] = tempCoord;
                }

            } else {
                AIM_ERROR(aimInfo, "Unable to rotate mesh!\n");
                status = CAPS_NOTFOUND;
                goto cleanup;
            }
        }

        status = extrude_SurfaceMesh(extrusion, extrusionBCIndex, &surfaceMesh, &volumeMesh);
        AIM_STATUS(aimInfo, status);

        strcpy(filename, projectName);

        // Write AFLR3
    /*@-nullpass@*/
         status = mesh_writeAFLR3(aimInfo, filename,
                                  0, // write binary file
                                  &volumeMesh,
                                  1.0);
    /*@+nullpass@*/
         AIM_STATUS(aimInfo, status);

        status = CAPS_SUCCESS;

    cleanup:
        if (status != CAPS_SUCCESS) {
            printf("Error: Premature exit in fun3d_2DMesh status = %d\n", status);
        }

    /*@-dependenttrans@*/
        if (fp != NULL) fclose(fp);
    /*@+dependenttrans@*/

        // Destroy meshes
        (void) destroy_meshStruct(&surfaceMesh);
        (void) destroy_meshStruct(&volumeMesh);

        return status;
    }


    // Remove unused nodes from dataMatrix and update connectivity matrix
    static int fun3d_removeUnused(void *aimInfo, int numVariable, int *numNode,  int used[],
                                  double ***data, /*@null@*/ int *numConnect,
                                  /*@null@*/ int *dataConnectMatrix)
    {
        int i, j, k; //Indexing
        int status;

        int *usedNode=NULL; // Freeable

        double **dataMatrix;

        dataMatrix = *data;

        // Copy used node array
        usedNode = (int *) EG_alloc(*numNode*sizeof(int));
        if (usedNode == NULL) {
            status = EGADS_MALLOC;
            goto cleanup;
        }
        for (i = 0; i < *numNode; i++) usedNode[i] = used[i];

        // Remove unused nodes
        j = 0;
        for (i = 0; i< *numNode; i++ ) {
            if (usedNode[i] == (int) false || usedNode[i] < 0) continue;

            usedNode[i] = j+1; // Set i-th node to essentially the node ID,  1-bias

            j +=1;
        }

        j = 0;
        for (i = 0; i< *numNode; i++ ) {
            if (usedNode[i] == (int) false || usedNode[i] < 0) continue;

            for (k = 0; k < numVariable; k++) {
                dataMatrix[k][j] = dataMatrix[k][i]; //Re-order dataMatrix - bubbling i'th index to j
            }

            j += 1;
        }

        *numNode = j; // New number of nodes

        // Redo connectivity
        if (dataConnectMatrix != NULL) {
            AIM_NOTNULL(numConnect, aimInfo, status);
            j = 0;
            for (i = 0; i < *numConnect; i++) {

                if (usedNode[dataConnectMatrix[4*i+ 0]-1] == (int) false) continue;
                if (usedNode[dataConnectMatrix[4*i+ 1]-1] == (int) false) continue;
                if (usedNode[dataConnectMatrix[4*i+ 2]-1] == (int) false) continue;
                if (usedNode[dataConnectMatrix[4*i+ 3]-1] == (int) false) continue;

                for (k = 0; k < 4; k++) {
                   dataConnectMatrix[4*j+ k] = usedNode[dataConnectMatrix[4*i+ k]-1];
                }

                j += 1;
            }

            *numConnect = j; // New number of elements
        }

        status = CAPS_SUCCESS;

    cleanup:
        if (status != CAPS_SUCCESS)
            printf("Error: Premature exit in fun3d_removeUnused status = %d\n",
                   status);

        if (usedNode != NULL) EG_free(usedNode);
        return status;
    }


    // Write FUN3D data transfer files
    int fun3d_dataTransfer(void *aimInfo,
                           const char *projectName,
                           const mapAttrToIndexStruct *groupMap,
                           const cfdBoundaryConditionStruct bcProps,
                           aimMeshRef *meshRef,
                           /*@null@*/ cfdModalAeroelasticStruct *eigenVector)
    {

        /*! \page dataTransferFUN3D FUN3D Data Transfer
         *
         * \section dataToFUN3D Data transfer to FUN3D (FieldIn)
         *
         * <ul>
         *  <li> <B>"Displacement"</B> </li> <br>
         *   Retrieves nodal displacements (as from a structural solver)
         *   and updates FUN3D's surface mesh; a new [project_name]_body1.dat file is written out which may
         *   be loaded into FUN3D to update the surface mesh/move the volume mesh using the FUN3D command line option
         *   -\-read_surface_from_file
         * </ul>
         *
         * <ul>
         *  <li> <B>"EigenVector_#"</B> </li> <br>
         *   Retrieves modal eigen-vectors from a structural solver, where "#" should be replaced by the
         *   corresponding mode number for the eigen-vector (eg. EigenVector_3 would correspond to the third mode,
         *   while EigenVector_6 would be the sixth mode) . A [project_name]_body1_mode#.dat file is written
         *   out for each mode.
         * </ul>
         *
         */

        int status; // Function return status
        int i, j, ibound, ibody, iface, iglobal, eigenIndex; // Indexing

        int stringLength = 0;

        char *filename = NULL;

        // Discrete data transfer variables
        capsDiscr *discr;
        char **boundName = NULL;
        int numBoundName = 0;
        enum capsdMethod dataTransferMethod;
        int numDataTransferPoint;
        int dataTransferRank;
        double *dataTransferData;
        char   *units;

        int state, nGlobal, *globalOffset=NULL, nFace;
        ego body, *faces=NULL;

        int alen, ntri, atype, itri, ielem;
        const double *face_xyz, *face_uv, *reals;
        const int *face_ptype, *face_pindex, *face_tris, *face_tric, *nquad=NULL;
        const char *string, *groupName = NULL;

        // Variables used in global node mapping
        int ptype, pindex;
        double xyz[3];

        // Data transfer Out variables
        const char *dataOutName[] = {"x","y","z", "id", "dx", "dy", "dz"};
        const int dataOutFormat[] = {Double, Double, Double, Integer, Double, Double, Double};

        double **dataOutMatrix = NULL;
        int *dataConnectMatrix = NULL;

        int numOutVariable = 7;
        int numOutDataPoint = 0;
        int numOutDataConnect = 0;

        const char fileExtBody[] = "_body1";
        const char fileExt[] = ".dat";
        const char fileExtMode[] = "_mode";

        int foundDisplacement = (int) false, foundEigenVector = (int) false;

        int marker;

        int numUsedNode = 0, numUsedConnectivity = 0, usedElems;
        int *usedNode = NULL;

        status = aim_getBounds(aimInfo, &numBoundName, &boundName);
        AIM_STATUS(aimInfo, status);

        foundDisplacement = foundEigenVector = (int) false;
        for (ibound = 0; ibound < numBoundName; ibound++) {
          AIM_NOTNULL(boundName, aimInfo, status);

          status = aim_getDiscr(aimInfo, boundName[ibound], &discr);
          if (status != CAPS_SUCCESS) continue;

          status = aim_getDataSet(discr,
                                  "Displacement",
                                  &dataTransferMethod,
                                  &numDataTransferPoint,
                                  &dataTransferRank,
                                  &dataTransferData,
                                  &units);

          if (status == CAPS_SUCCESS) { // If we do have data ready is the rank correct

            foundDisplacement = (int) true;

            if (dataTransferRank != 3) {
              AIM_ERROR(aimInfo, "Displacement transfer data found however rank is %d not 3!!!!\n", dataTransferRank);
              status = CAPS_BADRANK;
              goto cleanup;
            }
            break;
          }

          if (eigenVector != NULL) {
            for (eigenIndex = 0; eigenIndex < eigenVector->numEigenValue; eigenIndex++) {

              status = aim_getDataSet(discr,
                                      eigenVector->eigenValue[eigenIndex].name,
                                      &dataTransferMethod,
                                      &numDataTransferPoint,
                                      &dataTransferRank,
                                      &dataTransferData,
                                      &units);

              if (status == CAPS_SUCCESS) { // If we do have data ready is the rank correct

                foundEigenVector = (int) true;

                if (dataTransferRank != 3) {
                  AIM_ERROR(aimInfo, "EigenVector transfer data found however rank is %d not 3!!!!\n", dataTransferRank);
                  status = CAPS_BADRANK;
                  goto cleanup;
                }
                break;
              }
            } // Loop through EigenValues

            if (foundEigenVector == (int) true) break;

          } // If eigen-vectors provided
        } // Loop through transfer names

        if (foundDisplacement != (int) true && foundEigenVector != (int) true) {
          printf("Info: No recognized data transfer names found.\n");
          status = CAPS_NOTFOUND;
          goto cleanup;
        }

        // Ok looks like we have displacements/EigenVectors to get so lets continue
        printf("Writing FUN3D data transfer files\n");

        // Allocate data arrays that are going to be output
        AIM_ALLOC(dataOutMatrix, numOutVariable, double*, aimInfo, status);
        for (i = 0; i < numOutVariable; i++) dataOutMatrix[i] = NULL;

        AIM_ALLOC(globalOffset, meshRef->nmap+1, int, aimInfo, status);
        numOutDataPoint = 0;
        ielem = 0;
        globalOffset[0] = 0;
        for (i = 0; i < meshRef->nmap; i++) {
          if (meshRef->maps[i].tess == NULL) continue;

          status = EG_statusTessBody(meshRef->maps[i].tess, &body, &state, &nGlobal);
          AIM_STATUS(aimInfo, status);

          // re-allocate data arrays
          for (j = 0; j < numOutVariable; j++) {
            AIM_REALL(dataOutMatrix[j], numOutDataPoint + nGlobal, double, aimInfo, status);
          }
          AIM_REALL(usedNode, numOutDataPoint + nGlobal, int, aimInfo, status);
          for (j = globalOffset[0]; j < numOutDataPoint + nGlobal; j++) usedNode[j] = (int) false;


          for (iglobal = 0; iglobal < nGlobal; iglobal++) {
            status = EG_getGlobal(meshRef->maps[i].tess,
                                  iglobal+1, &ptype, &pindex, xyz);
            AIM_STATUS(aimInfo, status);

            // First just set the Coordinates
            dataOutMatrix[0][globalOffset[i]+iglobal] = xyz[0];
            dataOutMatrix[1][globalOffset[i]+iglobal] = xyz[1];
            dataOutMatrix[2][globalOffset[i]+iglobal] = xyz[2];

            // Volume mesh node ID
            dataOutMatrix[3][globalOffset[i]+iglobal] = meshRef->maps[i].map[iglobal];

            // Delta displacements
            dataOutMatrix[4][i] = 0;
            dataOutMatrix[5][i] = 0;
            dataOutMatrix[6][i] = 0;
          }

          // check if the tessellation has a mixture of quad and tess
          status = EG_attributeRet(meshRef->maps[i].tess, ".mixed",
                                   &atype, &alen, &nquad, &reals, &string);
          if (status != EGADS_SUCCESS &&
              status != EGADS_NOTFOUND) AIM_STATUS(aimInfo, status);

          status = EG_getBodyTopos(body, NULL, FACE, &nFace, &faces);
          AIM_STATUS(aimInfo, status);
          for (iface = 0; iface < nFace; iface++) {
            // get the face tessellation
            status = EG_getTessFace(meshRef->maps[i].tess, iface+1, &alen, &face_xyz, &face_uv,
                                    &face_ptype, &face_pindex, &ntri, &face_tris, &face_tric);
            AIM_STATUS(aimInfo, status);
            AIM_NOTNULL(faces, aimInfo, status);

            status = retrieve_CAPSGroupAttr(faces[iface], &groupName);
            if (status == EGADS_SUCCESS) {
              AIM_NOTNULL(groupName, aimInfo, status);
              status = get_mapAttrToIndexIndex(groupMap, groupName, &marker);
              if (status != CAPS_SUCCESS) {
                AIM_ERROR(aimInfo, "No capsGroup \"%s\" not found in attribute map", groupName);
                goto cleanup;
              }
            } else {
              AIM_ERROR(aimInfo, "No capsGroup on face %d", iface+1);
              print_AllAttr(aimInfo, faces[iface]);
              goto cleanup;
            }

            //  To keep with the moving_bodying input we will assume used nodes are all inviscid and viscous surfaces instead
            //                        usedNode[k] = (int) true;
            usedElems = (int) false;
            for (j = 0; j < bcProps.numSurfaceProp; j++) {
              if (marker != bcProps.surfaceProp[j].bcID) continue;

              if (bcProps.surfaceProp[j].surfaceType == Viscous ||
                  bcProps.surfaceProp[j].surfaceType == Inviscid) {
                  usedElems = (int) true;
              }
              break;
            }

            if (nquad == NULL) { // all triangles

              // re-allocate data arrays
              AIM_REALL(dataConnectMatrix, 4*(numOutDataConnect + ntri), int, aimInfo, status);

              for (itri = 0; itri < ntri; itri++, ielem++) {
                for (j = 0; j < 3; j++) {
                  status = EG_localToGlobal(meshRef->maps[i].tess, iface+1, face_tris[3*itri+j], &iglobal);
                  AIM_STATUS(aimInfo, status);
                  dataConnectMatrix[4*ielem+j] = globalOffset[i] + iglobal;
                  usedNode[globalOffset[i]+iglobal-1] = usedElems;
                }
                // repeat the last node for triangles
                dataConnectMatrix[4*ielem+3] = dataConnectMatrix[4*ielem+2];
              }

              numOutDataConnect += ntri;

            } else { // mixture of tri and quad elements

              // re-allocate data arrays
              AIM_REALL(dataConnectMatrix, 4*(numOutDataConnect + ntri-nquad[iface]), int, aimInfo, status);

              // process triangles
              for (itri = 0; itri < ntri-2*nquad[iface]; itri++, ielem++) {
                for (j = 0; j < 3; j++) {
                  status = EG_localToGlobal(meshRef->maps[i].tess, iface+1, face_tris[3*itri+j], &iglobal);
                  AIM_STATUS(aimInfo, status);
                  dataConnectMatrix[4*ielem+j] = globalOffset[i] + iglobal;
                  usedNode[globalOffset[i]+iglobal-1] = usedElems;
                }
                // repeat the last node for triangle
                dataConnectMatrix[4*ielem+3] = dataConnectMatrix[4*ielem+2];
              }
              // process quads
              for (; itri < ntri; itri++, ielem++) {
                for (j = 0; j < 3; j++) {
                  status = EG_localToGlobal(meshRef->maps[i].tess, iface+1, face_tris[3*itri+j], &iglobal);
                  AIM_STATUS(aimInfo, status);
                  dataConnectMatrix[4*ielem+j] = globalOffset[i] + iglobal;
                  usedNode[globalOffset[i]+iglobal-1] = usedElems;
                }

                // add the last node from the 2nd triangle to make the quad
                itri++;
                status = EG_localToGlobal(meshRef->maps[i].tess, iface+1, face_tris[3*itri+2], &iglobal);
                AIM_STATUS(aimInfo, status);
                dataConnectMatrix[4*ielem+3] = globalOffset[i] + iglobal;
                usedNode[globalOffset[i]+iglobal-1] = usedElems;
              }

              numOutDataConnect += ntri-nquad[iface];
            }

          }
          AIM_FREE(faces);

          numOutDataPoint += nGlobal;
          globalOffset[i+1] = globalOffset[i] + nGlobal;
        }

        // Re-loop through transfers - if we are doing displacements
        if (foundDisplacement == (int) true) {

          for (ibound = 0; ibound < numBoundName; ibound++) {
            AIM_NOTNULL(boundName, aimInfo, status);

            status = aim_getDiscr(aimInfo, boundName[ibound], &discr);
            if (status != CAPS_SUCCESS) continue;

            status = aim_getDataSet(discr,
                                    "Displacement",
                                    &dataTransferMethod,
                                    &numDataTransferPoint,
                                    &dataTransferRank,
                                    &dataTransferData,
                                    &units);
            if (status != CAPS_SUCCESS) continue; // If no elements in this object skip to next transfer name

            if (numDataTransferPoint != discr->nPoints &&
                numDataTransferPoint > 1) {
              AIM_ERROR(aimInfo, "Developer error!! %d != %d", numDataTransferPoint, discr->nPoints);
              status = CAPS_MISMATCH;
              goto cleanup;
            }

            for (i = 0; i < discr->nPoints; i++) {

              ibody   = discr->tessGlobal[2*i+0];
              iglobal = discr->tessGlobal[2*i+1];

              status = EG_getGlobal(discr->bodys[ibody-1].tess,
                                    iglobal, &ptype, &pindex, xyz);
              AIM_STATUS(aimInfo, status);

              // Find the disc tessellation in the original list of tessellations
              for (j = 0; j < meshRef->nmap; j++) {
                if (discr->bodys[ibody-1].tess == meshRef->maps[j].tess) {
                  break;
                }
              }
              if (j == meshRef->nmap) {
                AIM_ERROR(aimInfo, "Could not find matching tessellation!");
                status = CAPS_MISMATCH;
                goto cleanup;
              }

              if (numDataTransferPoint == 1) {
                // A single point means this is an initialization phase

                // Apply delta displacements
                dataOutMatrix[0][globalOffset[j]+iglobal-1] += dataTransferData[0];
                dataOutMatrix[1][globalOffset[j]+iglobal-1] += dataTransferData[1];
                dataOutMatrix[2][globalOffset[j]+iglobal-1] += dataTransferData[2];

                // save delta displacements
                dataOutMatrix[4][globalOffset[j]+iglobal-1] = dataTransferData[0];
                dataOutMatrix[5][globalOffset[j]+iglobal-1] = dataTransferData[1];
                dataOutMatrix[6][globalOffset[j]+iglobal-1] = dataTransferData[2];

              } else {
                // Apply delta displacements
                dataOutMatrix[0][globalOffset[j]+iglobal-1] += dataTransferData[3*i+0];
                dataOutMatrix[1][globalOffset[j]+iglobal-1] += dataTransferData[3*i+1];
                dataOutMatrix[2][globalOffset[j]+iglobal-1] += dataTransferData[3*i+2];

                dataOutMatrix[4][globalOffset[j]+iglobal-1] = dataTransferData[3*i+0];
                dataOutMatrix[5][globalOffset[j]+iglobal-1] = dataTransferData[3*i+1];
                dataOutMatrix[6][globalOffset[j]+iglobal-1] = dataTransferData[3*i+2];
              }
            }
          } // End numBoundName loop

          // Remove unused nodes
          numUsedNode = numOutDataPoint;
          numUsedConnectivity = numOutDataConnect;
          AIM_NOTNULL(usedNode, aimInfo, status);
          status = fun3d_removeUnused(aimInfo, numOutVariable, &numUsedNode, usedNode,
                                      &dataOutMatrix, &numUsedConnectivity,
                                      dataConnectMatrix);
          AIM_STATUS(aimInfo, status);

          stringLength = strlen(projectName) + strlen(fileExtBody) + strlen(fileExt) + 1;
          AIM_ALLOC(filename, stringLength+1, char, aimInfo, status);

          strcpy(filename, projectName);
          strcat(filename, fileExtBody);
          strcat(filename, fileExt);
          filename[stringLength] = '\0';

          // Write out displacement in tecplot file
          /*@-nullpass@*/
          status = tecplot_writeFEPOINT(aimInfo, filename,
                                        "FUN3D AeroLoads",
                                        NULL,
                                        numOutVariable,
                                        (char **)dataOutName,
                                        numUsedNode, // numOutDataPoint,
                                        dataOutMatrix,
                                        dataOutFormat,
                                        numUsedConnectivity, //numOutDataConnect, // numConnectivity
                                        dataConnectMatrix, // connectivity matrix
                                        NULL); // Solution time
          /*@+nullpass@*/
          AIM_STATUS(aimInfo, status);
          AIM_FREE(filename);
        } // End if found displacements

        // Re-loop through transfers - if we are doing eigen-vectors
        if ((foundEigenVector == (int) true) && (eigenVector != NULL)) {

          for (eigenIndex = 0; eigenIndex < eigenVector->numEigenValue; eigenIndex++) {

            // Zero out the eigen-vectors each time we are writing out a new one
            for (i = 0; i < numOutDataPoint; i++ ) {

              // Delta eigen-vectors
              dataOutMatrix[4][i] = 0;
              dataOutMatrix[5][i] = 0;
              dataOutMatrix[6][i] = 0;
            }

            for (ibound = 0; ibound < numBoundName; ibound++) {
              AIM_NOTNULL(boundName, aimInfo, status);

              status = aim_getDiscr(aimInfo, boundName[ibound], &discr);
              if (status != CAPS_SUCCESS) continue;

              status = aim_getDataSet(discr,
                                      eigenVector->eigenValue[eigenIndex].name,
                                      &dataTransferMethod,
                                      &numDataTransferPoint,
                                      &dataTransferRank,
                                      &dataTransferData,
                                      &units);
              if (status != CAPS_SUCCESS) continue; // If no elements in this object skip to next transfer name

              if (numDataTransferPoint != discr->nPoints &&
                  numDataTransferPoint > 1) {
                AIM_ERROR(aimInfo, "Developer error!! %d != %d", numDataTransferPoint, discr->nPoints);
                status = CAPS_MISMATCH;
                goto cleanup;
              }

              for (i = 0; i < discr->nPoints; i++) {

                ibody   = discr->tessGlobal[2*i+0];
                iglobal = discr->tessGlobal[2*i+1];

                // Find the disc tessellation in the original list of tessellations
                for (j = 0; j < meshRef->nmap; j++) {
                  if (discr->bodys[ibody-1].tess == meshRef->maps[j].tess) {
                    break;
                  }
                }
                if (j == meshRef->nmap) {
                  AIM_ERROR(aimInfo, "Could not find matching tessellation!");
                  status = CAPS_MISMATCH;
                  goto cleanup;
                }

                if (numDataTransferPoint == 1) {
                  // A single point means this is an initialization phase

                  // save Eigen-vector
                  dataOutMatrix[4][globalOffset[j]+iglobal-1] = dataTransferData[0];
                  dataOutMatrix[5][globalOffset[j]+iglobal-1] = dataTransferData[1];
                  dataOutMatrix[6][globalOffset[j]+iglobal-1] = dataTransferData[2];

                } else {
                  // save Eigen-vector
                  dataOutMatrix[4][globalOffset[j]+iglobal-1] = dataTransferData[3*i+0];
                  dataOutMatrix[5][globalOffset[j]+iglobal-1] = dataTransferData[3*i+1];
                  dataOutMatrix[6][globalOffset[j]+iglobal-1] = dataTransferData[3*i+2];
                }
              }
            } // End dataTransferDiscreteObj loop

            // Remove unused nodes
            numUsedNode = numOutDataPoint;
            AIM_NOTNULL(usedNode, aimInfo, status);

            if (eigenIndex == 0) {
              numUsedConnectivity = numOutDataConnect;
              status = fun3d_removeUnused(aimInfo, numOutVariable, &numUsedNode,
                                          usedNode, &dataOutMatrix,
                                          &numUsedConnectivity, dataConnectMatrix);
              AIM_STATUS(aimInfo, status);
            } else {
              status = fun3d_removeUnused(aimInfo, numOutVariable, &numUsedNode,
                                          usedNode, &dataOutMatrix, NULL, NULL);
              AIM_STATUS(aimInfo, status);
            }

            stringLength = strlen(projectName) +
                           strlen(fileExtBody) +
                           strlen(fileExtMode) +
                           strlen(fileExt) + 5;

            AIM_ALLOC(filename, stringLength+1, char, aimInfo, status);

            sprintf(filename, "%s%s%s%d%s",
                    projectName,
                    fileExtBody,
                    fileExtMode,  // Change modeNumber so it always starts at 1!
                    eigenIndex+1, // eigenVector->eigenValue[eigenIndex].modeNumber,
                    fileExt);

            // Write out eigen-vector in tecplot file
            /*@-nullpass@*/
            status = tecplot_writeFEPOINT(aimInfo,
                                          filename,
                                          "FUN3D Modal",
                                          NULL,
                                          numOutVariable,
                                          (char **)dataOutName,
                                          numUsedNode, //numOutDataPoint,
                                          dataOutMatrix,
                                          dataOutFormat,
                                          numUsedConnectivity, //numOutDataConnect, // numConnectivity
                                          dataConnectMatrix, // connectivity matrix
                                          NULL); // Solution time
            /*@+nullpass@*/
            AIM_STATUS(aimInfo, status);
            AIM_FREE(filename);

          } // End eigenvector names
        } // End if found eigenvectors

        status = CAPS_SUCCESS;

        // Clean-up
    cleanup:

        if (status != CAPS_SUCCESS &&
            status != CAPS_NOTFOUND) printf("Error: Premature exit in fun3d_dataTransfer status = %d\n", status);

        if (dataOutMatrix != NULL) {
          for (i = 0; i < numOutVariable; i++) {
            AIM_FREE(dataOutMatrix[i]);
          }
        }

        AIM_FREE(faces);
        AIM_FREE(dataOutMatrix);
        AIM_FREE(dataConnectMatrix);

        AIM_FREE(filename);
        AIM_FREE(boundName);
        AIM_FREE(usedNode);

        return status;
    }


    // Write FUN3D fun3d.nml file
    int fun3d_writeNML(void *aimInfo, capsValue *aimInputs, cfdBoundaryConditionStruct bcProps)
    {

        int status; // Function return status

        int i; // Indexing

        FILE *fnml = NULL;
        char filename[PATH_MAX];
        char fileExt[] ="fun3d.nml";

        printf("Writing fun3d.nml\n");
        if (aimInputs[Design_Functional-1].nullVal == NotNull ||
            aimInputs[Design_SensFile-1].vals.integer == (int)true) {
    #ifdef WIN32
            snprintf(filename, PATH_MAX, "Flow\\%s", fileExt);
    #else
            snprintf(filename, PATH_MAX, "Flow/%s", fileExt);
    #endif
        } else {
            strcpy(filename, fileExt);
        }

        fnml = aim_fopen(aimInfo, filename, "w");
        if (fnml == NULL) {
            AIM_ERROR(aimInfo, "Unable to open file - %s\n", filename);
            status = CAPS_IOERR;
            goto cleanup;
        }

        // &project
        fprintf(fnml,"&project\n");
        fprintf(fnml," project_rootname = \"%s\"\n",
                aimInputs[Proj_Name-1].vals.string);
        fprintf(fnml,"/\n\n");

        // &raw_grid
        fprintf(fnml,"&raw_grid\n");
        //fprintf(fnml," grid_format = \"%s\"\n",
        //        aimInputs[Mesh_Format-1].vals.string);

    //    if (aimInputs[Mesh_ASCII_Flag-1].vals.integer == (int) true) {
    //        fprintf(fnml," data_format = \"ascii\"\n");
    //    } else fprintf(fnml," data_format = \"stream\"\n");

        fprintf(fnml," grid_format = \"AFLR3\"\n");
        fprintf(fnml," data_format = \"stream\"\n");

        if (aimInputs[Two_Dimensional-1].vals.integer == (int) true) {
            fprintf(fnml," twod_mode = .true.\n");
          //fprintf(fnml," ignore_euler_number = .true.\n");
        }

        fprintf(fnml,"/\n\n");

        // &reference_physical_properties
        fprintf(fnml,"&reference_physical_properties\n");

        if (aimInputs[Mach-1].nullVal != IsNull) {
            fprintf(fnml," mach_number = %f\n", aimInputs[Mach-1].vals.real);
        }

        if (aimInputs[Re-1].nullVal != IsNull) {
            fprintf(fnml," reynolds_number = %f\n", aimInputs[Re-1].vals.real);
        }

        if (aimInputs[Alpha-1].nullVal != IsNull) {
            fprintf(fnml," angle_of_attack = %f\n", aimInputs[Alpha-1].vals.real);
        }

        if (aimInputs[Beta-1].nullVal != IsNull) {
            fprintf(fnml," angle_of_yaw = %f\n", aimInputs[Beta-1].vals.real);
        }

        if (aimInputs[Reference_Temperature-1].nullVal != IsNull) {
                fprintf(fnml," temperature = %f\n", aimInputs[Reference_Temperature-1].vals.real);

                if (aimInputs[Reference_Temperature-1].units != NULL) {
                    fprintf(fnml," temperature_units = \'%s\'\n", aimInputs[Reference_Temperature-1].units);
                }
        }

        fprintf(fnml,"/\n\n");

        // &governing_equations
        fprintf(fnml,"&governing_equations\n");

        if (aimInputs[Viscoux-1].nullVal != IsNull) {
            fprintf(fnml," viscous_terms = \"%s\"\n", aimInputs[Viscoux-1].vals.string);
        }

        if (aimInputs[Equation_Type-1].nullVal != IsNull) {
            fprintf(fnml," eqn_type = \"%s\"\n", aimInputs[Equation_Type-1].vals.string);
        }

        fprintf(fnml,"/\n\n");

        // &nonlinear_solver_parameters
        fprintf(fnml,"&nonlinear_solver_parameters\n");

        if (aimInputs[Time_Accuracy-1].nullVal != IsNull) {
            fprintf(fnml," time_accuracy = \"%s\"\n",
                    aimInputs[Time_Accuracy-1].vals.string);
        }

        if (aimInputs[Time_Step-1].nullVal != IsNull) {
            fprintf(fnml," time_step_nondim = %f\n", aimInputs[Time_Step-1].vals.real);
        }

        if (aimInputs[Num_Subiter-1].nullVal != IsNull) {
            fprintf(fnml," subiterations = %d\n",
                    aimInputs[Num_Subiter-1].vals.integer);
        }

        if (aimInputs[Temporal_Error-1].nullVal != IsNull) {

            fprintf(fnml," temporal_err_control = .true.\n");
            fprintf(fnml," temporal_err_floor = %f\n",
                    aimInputs[Temporal_Error-1].vals.real);

        }

        if (aimInputs[CFL_Schedule-1].nullVal != IsNull) {
            fprintf(fnml," schedule_cfl = %f %f\n",
                    aimInputs[CFL_Schedule-1].vals.reals[0],
                    aimInputs[CFL_Schedule-1].vals.reals[1]);
        }

        if (aimInputs[CFL_Schedule_Iter-1].nullVal != IsNull) {
            fprintf(fnml," schedule_iteration = %d %d\n",
                    aimInputs[CFL_Schedule_Iter-1].vals.integers[0],
                    aimInputs[CFL_Schedule_Iter-1].vals.integers[1]);
        }

        fprintf(fnml,"/\n\n");

        // &code_run_control
        fprintf(fnml,"&code_run_control\n");

        if (aimInputs[Num_Iter-1].nullVal != IsNull) {
            fprintf(fnml," steps = %d\n", aimInputs[Num_Iter-1].vals.integer);
        }

        if (aimInputs[Restart_Read-1].nullVal != IsNull) {
            fprintf(fnml," restart_read = '%s'\n",
                    aimInputs[Restart_Read-1].vals.string);
        }


        fprintf(fnml,"/\n\n");

        //&force_moment_integ_properties
        fprintf(fnml,"&force_moment_integ_properties\n");

        if (aimInputs[Reference_Area-1].nullVal != IsNull) {
            fprintf(fnml," area_reference = %f\n",
                    aimInputs[Reference_Area-1].vals.real);
        }

        if (aimInputs[Moment_Length-1].nullVal != IsNull) {
            fprintf(fnml," x_moment_length = %f\n",
                    aimInputs[Moment_Length-1].vals.reals[0]);

            fprintf(fnml," y_moment_length = %f\n",
                    aimInputs[Moment_Length-1].vals.reals[1]);
        }

        if (aimInputs[Moment_Center-1].nullVal != IsNull) {
            fprintf(fnml," x_moment_center = %f\n",
                    aimInputs[Moment_Center-1].vals.reals[0]);

            fprintf(fnml," y_moment_center = %f\n",
                    aimInputs[Moment_Center-1].vals.reals[1]);

            fprintf(fnml," z_moment_center = %f\n",
                    aimInputs[Moment_Center-1].vals.reals[2]);
        }

        fprintf(fnml,"/\n\n");

        //&boundary_conditions
        fprintf(fnml,"&boundary_conditions\n");

        // Loop through boundary conditions
        for (i = 0; i < bcProps.numSurfaceProp ; i++) {

            // Temperature
            if (bcProps.surfaceProp[i].wallTemperatureFlag == (int) true) {
                fprintf(fnml," wall_temperature(%d) = %f\n",bcProps.surfaceProp[i].bcID,
                        bcProps.surfaceProp[i].wallTemperature);
                fprintf(fnml," wall_temp_flag(%d) = .true.\n",bcProps.surfaceProp[i].bcID);
            }

            // Total pressure and temperature
            if (bcProps.surfaceProp[i].surfaceType == SubsonicInflow) {

                fprintf(fnml, " total_pressure_ratio(%d) = %f\n", bcProps.surfaceProp[i].bcID,
                        bcProps.surfaceProp[i].totalPressure);

                fprintf(fnml, " total_temperature_ratio(%d) = %f\n", bcProps.surfaceProp[i].bcID,
                        bcProps.surfaceProp[i].totalTemperature);
            }

            // Static pressure
            if (bcProps.surfaceProp[i].surfaceType == BackPressure ||
                    bcProps.surfaceProp[i].surfaceType == SubsonicOutflow) {

                fprintf(fnml, " static_pressure_ratio(%d) = %f\n", bcProps.surfaceProp[i].bcID,
                        bcProps.surfaceProp[i].staticPressure);
            }

            // Mach number
            if (bcProps.surfaceProp[i].surfaceType == MachOutflow ||
                    bcProps.surfaceProp[i].surfaceType == MassflowOut) {

                fprintf(fnml, " mach_bc(%d) = %f\n", bcProps.surfaceProp[i].bcID,
                        bcProps.surfaceProp[i].machNumber);
            }

            // Massflow
            if (bcProps.surfaceProp[i].surfaceType == MassflowIn ||
                    bcProps.surfaceProp[i].surfaceType == MassflowOut) {

                fprintf(fnml, " massflow(%d) = %f\n", bcProps.surfaceProp[i].bcID,
                        bcProps.surfaceProp[i].massflow);
            }

            // Fixed inflow and outflow
            /*if (bcProps.surfaceProp[i].surfaceType == FixedInflow ||
                bcProps.surfaceProp[i].surfaceType == FixedOutflow) {

                fprintf(fnml, " qset(%d,1) = %f\n", bcProps.surfaceProp[i].bcID,
                                                    bcProps.surfaceProp[i].staticDensity);

                fprintf(fnml, " qset(%d,2) = %f\n", bcProps.surfaceProp[i].bcID,
                                                    bcProps.surfaceProp[i].uVelocity);

                fprintf(fnml, " qset(%d,3) = %f\n", bcProps.surfaceProp[i].bcID,
                                                    bcProps.surfaceProp[i].vVelocity);

                fprintf(fnml, " qset(%d,4) = %f\n", bcProps.surfaceProp[i].bcID,
                                                    bcProps.surfaceProp[i].wVelocity);

                fprintf(fnml, " qset(%d,5) = %f\n", bcProps.surfaceProp[i].bcID,
                                                    bcProps.surfaceProp[i].staticDensity);
            }*/
        }

        fprintf(fnml,"/\n\n");

        // &noninertial_reference_frame
        fprintf(fnml,"&noninertial_reference_frame\n");


        if (aimInputs[NonInertial_Rotation_Rate-1].nullVal   != IsNull ||
            aimInputs[NonInertial_Rotation_Center-1].nullVal != IsNull) {

            fprintf(fnml," noninertial = .true.\n");
        }

        if (aimInputs[NonInertial_Rotation_Center-1].nullVal != IsNull) {
            fprintf(fnml," rotation_center_x = %f\n",
                    aimInputs[NonInertial_Rotation_Center-1].vals.reals[0]);
            fprintf(fnml," rotation_center_y = %f\n",
                    aimInputs[NonInertial_Rotation_Center-1].vals.reals[1]);
            fprintf(fnml," rotation_center_z = %f\n",
                    aimInputs[NonInertial_Rotation_Center-1].vals.reals[2]);
        }

        if (aimInputs[NonInertial_Rotation_Rate-1].nullVal != IsNull) {
            fprintf(fnml," rotation_rate_x = %f\n",
                    aimInputs[NonInertial_Rotation_Rate-1].vals.reals[0]);
            fprintf(fnml," rotation_rate_y = %f\n",
                    aimInputs[NonInertial_Rotation_Rate-1].vals.reals[1]);
            fprintf(fnml," rotation_rate_z = %f\n",
                    aimInputs[NonInertial_Rotation_Rate-1].vals.reals[2]);

        }

        fprintf(fnml,"/\n\n");

        status = CAPS_SUCCESS;

    cleanup:

        if (status != CAPS_SUCCESS)
            printf("Error: Premature exit in fun3d_writeNML status = %d\n", status);

        if (fnml != NULL) fclose(fnml);

        return status;
    }


    // Write FUN3D movingbody.input file
    int fun3d_writeMovingBody(void *aimInfo, double fun3dVersion, cfdBoundaryConditionStruct bcProps,
                              cfdModalAeroelasticStruct *modalAeroelastic)
    {

        int status; // Function return status

        int i, eigenIndex; // Indexing
        int counter = 0;
        FILE *fp = NULL;
        char *filename = NULL;
        char fileExt[] ="moving_body.input";

        int bodyIndex = 1;
        int stringLength;

        printf("Writing moving_body.input");

        stringLength = strlen(fileExt) + 1;

        AIM_ALLOC(filename,stringLength +1, char, aimInfo, status);

        strcpy(filename, fileExt);
        filename[stringLength] = '\0';

        fp = aim_fopen(aimInfo, filename, "w");
        if (fp == NULL) {
            AIM_ERROR(aimInfo, "Unable to open file - %s\n", filename);
            status = CAPS_IOERR;
            goto cleanup;
        }

        // &body_definitions
        fprintf(fp,"&body_definitions\n");

        fprintf(fp," n_moving_bodies = %d\n", bodyIndex);

        counter = 0;
        for (i = 0; i < bcProps.numSurfaceProp; i++) {

            if (bcProps.surfaceProp[i].surfaceType == Viscous ||
                    bcProps.surfaceProp[i].surfaceType == Inviscid) {

                fprintf(fp," defining_bndry(%d,%d) = %d\n", counter+1,
                        bodyIndex,
                        bcProps.surfaceProp[i].bcID);

                counter += 1;
            }
        }

        fprintf(fp," n_defining_bndry(%d) = %d\n", bodyIndex, counter);

        fprintf(fp," motion_driver(%d) = ", bodyIndex);
        if (modalAeroelastic != NULL) {
            fprintf(fp,"\"aeroelastic\"\n");
        }

        fprintf(fp," mesh_movement(%d) = ", bodyIndex);
        if (modalAeroelastic != NULL) {
            fprintf(fp,"\"deform\"\n");
        }

        fprintf(fp,"/\n\n");

        if (modalAeroelastic != NULL) {

            // &aeroelastic_modal_data
            fprintf(fp,"&aeroelastic_modal_data\n");

            fprintf(fp," nmode(%d) = %d\n", bodyIndex, modalAeroelastic->numEigenValue);

            if (fun3dVersion < 13.1) {
                fprintf(fp," uinf(%d) = %f\n", bodyIndex, modalAeroelastic->freestreamVelocity);
                fprintf(fp," qinf(%d) = %f\n", bodyIndex, modalAeroelastic->freestreamDynamicPressure);
                fprintf(fp," grefl(%d) = %f\n", bodyIndex, modalAeroelastic->lengthScaling);
            } else {
                fprintf(fp," uinf = %f\n", modalAeroelastic->freestreamVelocity);
                fprintf(fp," qinf = %f\n", modalAeroelastic->freestreamDynamicPressure);
                fprintf(fp," grefl = %f\n",modalAeroelastic->lengthScaling);
            }

            fprintf(fp, "\n");


            for (i = 0; i < modalAeroelastic->numEigenValue; i++) {

                eigenIndex = i + 1; // Change mode number so that it always starts at 1
                // modalAeroelastic->eigenValue[i].modeNumber

                fprintf(fp, " ! Mode %d of %d (structural mode %d)\n", eigenIndex,
                        modalAeroelastic->numEigenValue,
                        modalAeroelastic->eigenValue[i].modeNumber);

                fprintf(fp," freq(%d,%d) = %f\n", eigenIndex, bodyIndex, modalAeroelastic->eigenValue[i].frequency);
                fprintf(fp," damp(%d,%d) = %f\n", eigenIndex, bodyIndex, modalAeroelastic->eigenValue[i].damping);

                fprintf(fp," gmass(%d,%d) = %f\n"  , eigenIndex, bodyIndex, modalAeroelastic->eigenValue[i].generalMass);
                fprintf(fp," gdisp0(%d,%d) = %f\n" , eigenIndex, bodyIndex, modalAeroelastic->eigenValue[i].generalDisplacement);
                fprintf(fp," gvel0(%d,%d) = %f\n"  , eigenIndex, bodyIndex, modalAeroelastic->eigenValue[i].generalVelocity);
                fprintf(fp," gforce0(%d,%d) = %f\n", eigenIndex, bodyIndex, modalAeroelastic->eigenValue[i].generalForce);

                fprintf(fp, "\n");
            }
            fprintf(fp,"/\n\n");
        }

        status = CAPS_SUCCESS;

    cleanup:
        if (status != CAPS_SUCCESS)
            printf("Error: Premature exit in fun3d_writeMovingBody status = %d\n",
                   status);

        if (fp != NULL) fclose(fp);

        if (filename != NULL) EG_free(filename);

        return status;
    }


    // Write FUN3D parameterization/sensitivity file
    // Will not calculate shape sensitivities if there are no geometry design variable; will
    // simple check and dump out the body meshes in model.tec files
    int  fun3d_writeParameterization(void *aimInfo,
                                     int numDesignVariable,
                                     cfdDesignVariableStruct designVariable[],
                                     aimMeshRef *meshRef)
    {

        int status; // Function return status

        int i, j, k, m, row, col; // Indexing

        int stringLength = 7;

        // Data transfer Out variables
        char **dataOutName= NULL;

        double ***dataOutMatrix = NULL;
        int *dataOutFormat = NULL;
        int *dataConnectMatrix = NULL;

        int numOutVariable = 4; // x, y, z, id, ... + 3* active GeomIn
        int numOutDataPoint = 0;
        int numOutDataConnect = 0;

        // Variables used in global node mapping
        int ptype, pindex;
        double xyz[3];

        const char *geomInName;
        int numPoint;
        double *dxyz = NULL;

        int index;
        int iface, iglobal;
        int state, nFace;
        ego body, *faces=NULL;

        int alen, ntri, atype, itri, ielem;
        const double *face_xyz, *face_uv, *reals;
        const int *face_ptype, *face_pindex, *face_tris, *face_tric, *nquad=NULL;
        const char *string;

        char message[100];
        char filePre[] = "model.tec.";
        char fileExt[] = ".sd1";
        char *filename = NULL;
        char folder[]  = "Rubberize";
        char zoneTitle[100];

        capsValue *geomInVal;
        int *geomSelect = NULL;

        AIM_ALLOC(geomSelect, numDesignVariable, int, aimInfo, status);

        // Determine number of geometry input variables

        for (i = 0; i < numDesignVariable; i++) {
            geomSelect[i] = (int) false;

            printf("DesignVariable = %s\n", designVariable[i].name);

            index = aim_getIndex(aimInfo, designVariable[i].name, GEOMETRYIN);
            if (index == CAPS_NOTFOUND) continue;
            if (index < CAPS_SUCCESS ) {
              status = index;
              AIM_STATUS(aimInfo, status);
            }

            if(aim_getGeomInType(aimInfo, index) != 0) {
                AIM_ERROR(aimInfo, "GeometryIn value %s is a configuration parameter and not a valid design parameter - can't get sensitivity\n",
                          designVariable[i].name);
                status = CAPS_BADVALUE;
                goto cleanup;
            }


            status = aim_getValue(aimInfo, index, GEOMETRYIN, &geomInVal);
            AIM_STATUS(aimInfo, status);

            numOutVariable += 3*geomInVal->length; // xD1, yD1, zD1, ...

            // Don't compute sensitivities if not needed
            geomSelect[i] = (int) true;
        }

        // No need to write Rubberize files without GeometryIn design variables
        if (numOutVariable == 4) {
            status = CAPS_SUCCESS;
            goto cleanup;
        }

        if (numOutVariable > 99999999) {
            AIM_ERROR(aimInfo, "Array of design variable names will be over-run!");
            status = CAPS_RANGEERR;
            goto cleanup;
        }

        // Allocate our names
        AIM_ALLOC(dataOutName, numOutVariable, char*, aimInfo, status);
        for (i = 0; i < numOutVariable; i++) dataOutName[i] = NULL;

        stringLength = 11;
        j = 1;
        k = 1;
        for (i = 0; i < numOutVariable; i++) {
            AIM_ALLOC(dataOutName[i], stringLength+1, char, aimInfo, status);

            // Set names
            if      (i == 0) sprintf(dataOutName[i], "%s", "x");
            else if (i == 1) sprintf(dataOutName[i], "%s", "y");
            else if (i == 2) sprintf(dataOutName[i], "%s", "z");
            else if (i == 3) sprintf(dataOutName[i], "%s", "id");
            else {

                if      (j == 1) sprintf(dataOutName[i], "%s%d", "xD", k);
                else if (j == 2) sprintf(dataOutName[i], "%s%d", "yD", k);
                else if (j == 3) {
                    sprintf(dataOutName[i], "%s%d", "zD", k);
                    j = 0;
                    k += 1;
                }
                j += 1;
            }
        }

        AIM_ALLOC(dataOutFormat, numOutVariable, int, aimInfo, status);

        // Set data out formatting
        for (i = 0; i < numOutVariable; i++) {
            if (strcasecmp(dataOutName[i], "id") == 0) {
                dataOutFormat[i] = (int) Integer;
            } else {
                dataOutFormat[i] = (int) Double;
            }
        }

        // Allocate data arrays that are going to be output
        AIM_ALLOC(dataOutMatrix, meshRef->nmap, double**, aimInfo, status);
        for (i = 0; i < meshRef->nmap; i++) dataOutMatrix[i] = NULL;

        for (i = 0; i < meshRef->nmap; i++) {
          AIM_ALLOC(dataOutMatrix[i], numOutVariable, double*, aimInfo, status);
          for (j = 0; j < numOutVariable; j++) dataOutMatrix[i][j] = NULL;
        }

        for (i = 0; i < meshRef->nmap; i++) {
          if (meshRef->maps[i].tess == NULL) continue;

          status = EG_statusTessBody(meshRef->maps[i].tess, &body, &state, &numOutDataPoint);
          AIM_STATUS(aimInfo, status);

          // allocate data arrays
          for (j = 0; j < numOutVariable; j++) {
            AIM_ALLOC(dataOutMatrix[i][j], numOutDataPoint, double, aimInfo, status);
          }

          for (iglobal = 0; iglobal < numOutDataPoint; iglobal++) {
            status = EG_getGlobal(meshRef->maps[i].tess,
                                  iglobal+1, &ptype, &pindex, xyz);
            AIM_STATUS(aimInfo, status);

            // First just set the Coordinates
            dataOutMatrix[i][0][iglobal] = xyz[0];
            dataOutMatrix[i][1][iglobal] = xyz[1];
            dataOutMatrix[i][2][iglobal] = xyz[2];

            // Volume mesh node ID
            dataOutMatrix[i][3][iglobal] = meshRef->maps[i].map[iglobal];
          }
        }

        // Loop over the geometry in values and compute sensitivities for all bodies
        m = 4;
        for (j = 0; j < numDesignVariable; j++) {

          if (geomSelect[j] == (int) false) continue;

          geomInName = designVariable[j].name;
          index = aim_getIndex(aimInfo, geomInName, GEOMETRYIN);

          status = aim_getValue(aimInfo, index, GEOMETRYIN, &geomInVal);
          AIM_STATUS(aimInfo, status);

          for (row = 0; row < geomInVal->nrow; row++) {
            for (col = 0; col < geomInVal->ncol; col++) {

              for (i = 0; i < meshRef->nmap; i++) {
                if (meshRef->maps[i].tess == NULL) continue;
                status = aim_tessSensitivity(aimInfo,
                                             geomInName,
                                             row+1, col+1, // row, col
                                             meshRef->maps[i].tess,
                                             &numPoint, &dxyz);
                AIM_STATUS(aimInfo, status, "Sensitivity for: %s\n", geomInName);
                AIM_NOTNULL(dxyz, aimInfo, status);

                for (k = 0; k < numPoint; k++) {
                  dataOutMatrix[i][m+0][k] = dxyz[3*k + 0]; // dx/dGeomIn
                  dataOutMatrix[i][m+1][k] = dxyz[3*k + 1]; // dy/dGeomIn
                  dataOutMatrix[i][m+2][k] = dxyz[3*k + 2]; // dz/dGeomIn
                }
                AIM_FREE(dxyz);
              }
              m += 3;
            }
          }
        }

        // Write sensitivity files for each body tessellation

        for (i = 0; i < meshRef->nmap; i++) {
          if (meshRef->maps[i].tess == NULL) continue;
          status = EG_statusTessBody(meshRef->maps[i].tess, &body, &state, &numOutDataPoint);
          AIM_STATUS(aimInfo, status);

          // check if the tessellation has a mixture of quad and tess
          status = EG_attributeRet(meshRef->maps[i].tess, ".mixed",
                                   &atype, &alen, &nquad, &reals, &string);
          if (status != EGADS_SUCCESS &&
              status != EGADS_NOTFOUND) AIM_STATUS(aimInfo, status);

          status = EG_getBodyTopos(body, NULL, FACE, &nFace, &faces);
          AIM_STATUS(aimInfo, status);

          ielem = 0;
          numOutDataConnect = 0;
          for (iface = 0; iface < nFace; iface++) {
            // get the face tessellation
            status = EG_getTessFace(meshRef->maps[i].tess, iface+1, &alen, &face_xyz, &face_uv,
                                    &face_ptype, &face_pindex, &ntri, &face_tris, &face_tric);
            AIM_STATUS(aimInfo, status);

            if (nquad == NULL) { // all triangles

              // re-allocate data arrays
              AIM_REALL(dataConnectMatrix, 4*(numOutDataConnect + ntri), int, aimInfo, status);

              for (itri = 0; itri < ntri; itri++, ielem++) {
                for (j = 0; j < 3; j++) {
                  status = EG_localToGlobal(meshRef->maps[i].tess, iface+1, face_tris[3*itri+j], &iglobal);
                  AIM_STATUS(aimInfo, status);
                  dataConnectMatrix[4*ielem+j] = iglobal;
                }
                // repeat the last node for triangles
                dataConnectMatrix[4*ielem+3] = dataConnectMatrix[4*ielem+2];
              }

              numOutDataConnect += ntri;

            } else { // mixture of tri and quad elements

              // re-allocate data arrays
              AIM_REALL(dataConnectMatrix, 4*(numOutDataConnect + ntri-nquad[iface]), int, aimInfo, status);

              // process triangles
              for (itri = 0; itri < ntri-2*nquad[iface]; itri++, ielem++) {
                for (j = 0; j < 3; j++) {
                  status = EG_localToGlobal(meshRef->maps[i].tess, iface+1, face_tris[3*itri+j], &iglobal);
                  AIM_STATUS(aimInfo, status);
                  dataConnectMatrix[4*ielem+j] = iglobal;
                }
                // repeat the last node for triangle
                dataConnectMatrix[4*ielem+3] = dataConnectMatrix[4*ielem+2];
              }
              // process quads
              for (; itri < ntri; itri++, ielem++) {
                for (j = 0; j < 3; j++) {
                  status = EG_localToGlobal(meshRef->maps[i].tess, iface+1, face_tris[3*itri+j], &iglobal);
                  AIM_STATUS(aimInfo, status);
                  dataConnectMatrix[4*ielem+j] = iglobal;
                }

                // add the last node from the 2nd triangle to make the quad
                itri++;
                status = EG_localToGlobal(meshRef->maps[i].tess, iface+1, face_tris[3*itri+2], &iglobal);
                AIM_STATUS(aimInfo, status);
                dataConnectMatrix[4*ielem+3] = iglobal;
              }

              numOutDataConnect += ntri-nquad[iface];
            }

          }
          AIM_FREE(faces);

          sprintf(message,"%s %d,", "sensitivity file for body", i+1);

          stringLength = strlen(folder) + 1 + strlen(filePre) + 7 + strlen(fileExt) + 1 ;

          AIM_REALL(filename, stringLength, char, aimInfo, status);

    #ifdef WIN32
          sprintf(filename, "%s\\%s%d%s", folder, filePre, i+1, fileExt);
    #else
          sprintf(filename, "%s/%s%d%s",  folder, filePre, i+1, fileExt);
    #endif

          sprintf(zoneTitle, "%s_%d", "Body", i+1);
          /*@-nullpass@*/
          status = tecplot_writeFEPOINT(aimInfo,
                                        filename,
                                        message,
                                        zoneTitle,
                                        numOutVariable,
                                        dataOutName,
                                        numOutDataPoint,
                                        dataOutMatrix[i],
                                        dataOutFormat,
                                        numOutDataConnect,
                                        dataConnectMatrix,
                                        NULL);
          /*@+nullpass@*/
          AIM_STATUS(aimInfo, status);
        }

        status = CAPS_SUCCESS;

    cleanup:
        (void) string_freeArray(numOutVariable, &dataOutName);
    #ifdef S_SPLINT_S
        EG_free(dataOutName);
    #endif

        if (dataOutMatrix != NULL) {
            for (i = 0; i < meshRef->nmap; i++) {
              if (dataOutMatrix[i] != NULL) {
                  for (j = 0; j < numOutVariable; j++) {
                      AIM_FREE(dataOutMatrix[i][j]);
                  }
              }
              AIM_FREE(dataOutMatrix[i]);
            }
        }

        AIM_FREE(dataOutMatrix);
        AIM_FREE(dataOutFormat);
        AIM_FREE(dataConnectMatrix);
        AIM_FREE(faces);

        AIM_FREE(dxyz);
        AIM_FREE(filename);
        AIM_FREE(geomSelect);

        return status;
    }


    static int _writeFunctinoalComponent(void *aimInfo, FILE *fp, cfdDesignFunctionalCompStruct *comp)
    {
        int status = CAPS_SUCCESS;

        const char *names[] = {"cl", "cd",                   // Lift, drag coefficients
                               "clp", "cdp",                 // Lift, drag coefficients: pressure contributions
                               "clv", "cdv",                 // Lift, drag coefficients: shear contributions
                               "cmx", "cmy", "cmz",          // x/y/z-axis moment coefficients
                               "cmxp", "cmyp", "cmzp",       // x/y/z-axis moment coefficients: pressure contributions
                               "cmxv", "cmyv", "cmzv",       // x/y/z-axis moment coefficients: shear contributions
                               "cx", "cy", "cz",             // x/y/z-axis force coefficients
                               "cxp", "cyp", "czp",          // x/y/z-axis force coefficients: pressure contributions
                               "cxv", "cyv", "czv",          // x/y/z-axis force coefficients: shear contributions
                               "powerx", "powery", "powerz", // x/y/z-axis power coefficients
                               "clcd",                       // Lift-to-drag ratio
                               "fom"    ,                    // Rotorcraft figure of merit
                               "propeff",                    // Rotorcraft propulsive efficiency
                               "rtr"    ,                    // thrust Rotorcraft thrust function
                               "pstag"  ,                    // RMS of stagnation pressure in cutting plane disk
                               "distort",                    // Engine inflow distortion
                               "boom"   ,                    // targ Near-field p/p∞ pressure target
                               "sboom"  ,                    // Coupled sBOOM ground-based noise metrics
                               "ae"     ,                    // Supersonic equivalent area target distribution
                               "press"  ,                    // box RMS of pressure in user-defined box, also pointwise dp/dt, dρ/dt
                               "cpstar" ,                    // Target pressure distributions
                               "sgen"                        // Entropy generation
                              };

        int i, found = (int)false;
        char *function=NULL, *c=NULL;

        for (i = 0; i < sizeof(names)/sizeof(char*); i++)
            if ( strcasecmp(names[i], comp->name) == 0 ) {
                found = (int)true;
                break;
            }

        if (found == (int) false) {
            AIM_ERROR(aimInfo, "Unknown function: '%s'", comp->name);
            AIM_ADDLINE(aimInfo, "Available functions:");
            for (i = 0; i < sizeof(names)/sizeof(char*); i++)
                AIM_ADDLINE(aimInfo, "'%s'", names[i]);
            status = CAPS_BADVALUE;
            goto cleanup;
        }

        // make the name lower case
        AIM_STRDUP(function, comp->name, aimInfo, status);
        for (c = function; *c != '\0'; ++c) *c = tolower(*c);

        // fprintf(fp, "Components of function   1: boundary id (0=all)/name/value/weight/target/power\n");
        fprintf(fp, " %d %s          %f %f %f %f\n", comp->bcID,
                                                     function,
                                                     0.0,
                                                     comp->weight,
                                                     comp->target,
                                                     comp->power);

        status = CAPS_SUCCESS;

    cleanup:
        AIM_FREE(function);
        return status;
    }


    // Write FUN3D  rubber.data file
    // Will not write shape entries unless explicitly told to check if they are need
    int fun3d_writeRubber(void *aimInfo,
                          cfdDesignStruct design,
                          int checkGeomShape,
                          double fun3dVersion,
                          aimMeshRef *meshRef)
    {
        int status; // Function return status

        int i, j, k, m; // Indexing

        int numBody = 0, numShapeVar = 0;

        char file[] = "rubber.data";
        FILE *fp = NULL;

        printf("Writing %s \n", file);
        if (meshRef == NULL) return CAPS_NULLVALUE;

        // Open file
        fp = aim_fopen(aimInfo, file, "w");
        if (fp == NULL) {
            AIM_ERROR(aimInfo, "Unable to open file: %s\n", file);
            status = CAPS_IOERR;
            goto cleanup;
        }

        numBody = meshRef->nmap;

        if (checkGeomShape == (int) true) {
            // Determine number of geometry input variables
            for (i = 0; i < design.numDesignVariable; i++) {

                if (design.designVariable[i].type != DesignVariableGeometry) continue;

                numShapeVar += design.designVariable[i].var->length; // xD1, yD1, zD1, ...
            }
        }

        // Don't write out body information if there are no geometry design variables
        if (numShapeVar == 0) numBody = 0;


        fprintf(fp, "################################################################################\n");
        fprintf(fp, "########################### Design Variable Information ########################\n");
        fprintf(fp, "################################################################################\n");
        fprintf(fp, "Global design variables (Mach number, AOA, Yaw, Noninertial rates)\n");
        fprintf(fp, "Var Active         Value               Lower Bound            Upper Bound\n");

        for (i = 0; i < design.numDesignVariable; i++) {

            if (strcasecmp(design.designVariable[i].name, "Mach") != 0) continue;

            if (design.designVariable[i].var->length < 1) {
                status = CAPS_RANGEERR;
                goto cleanup;
            }

            fprintf(fp, "Mach    1   %.15E  %.15E  %.15E\n", design.designVariable[i].value[0],
                                                             design.designVariable[i].lowerBound[0],
                                                             design.designVariable[i].upperBound[0]);
            break;
        }

        if (i >= design.numDesignVariable) fprintf(fp, "Mach    0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00\n");

        for (i = 0; i < design.numDesignVariable; i++) {

            if (strcasecmp(design.designVariable[i].name, "Alpha") != 0) continue;

            if (design.designVariable[i].var->length < 1) {
                status = CAPS_RANGEERR;
                goto cleanup;
            }

            fprintf(fp, "AOA     1   %.15E  %.15E  %.15E\n", design.designVariable[i].value[0],
                                                             design.designVariable[i].lowerBound[0],
                                                             design.designVariable[i].upperBound[0]);
            break;
        }

        if (i >= design.numDesignVariable) fprintf(fp, "AOA     0   0.000000000000000E+00  0.000000000000000E+00  0.00000000000000E+00\n");


        if (fun3dVersion > 12.4) {
            // FUN3D version 13.1 - version 12.4 doesn't have these available

            for (i = 0; i < design.numDesignVariable; i++) {

                if (strcasecmp(design.designVariable[i].name, "Beta") != 0) continue;

                if (design.designVariable[i].var->length < 1) {
                    status = CAPS_RANGEERR;
                    goto cleanup;
                }

                fprintf(fp, "Yaw     1   %.15E  %.15E  %.15E\n", design.designVariable[i].value[0],
                                                                 design.designVariable[i].lowerBound[0],
                                                                 design.designVariable[i].upperBound[0]);
                break;
            }

            if (i >= design.numDesignVariable) fprintf(fp, "Yaw     0   0.000000000000000E+00  0.000000000000000E+00  0.00000000000000E+00\n");

            for (i = 0; i < design.numDesignVariable; i++) {

                if (strcasecmp(design.designVariable[i].name, "NonInertial_Rotation_Rate") != 0) continue;

                if (design.designVariable[i].var->length < 3) {
                    status = CAPS_RANGEERR;
                    goto cleanup;
                }

                fprintf(fp, "xrate   1   %.15E  %.15E  %.15E\n", design.designVariable[i].value[0],
                                                                 design.designVariable[i].lowerBound[0],
                                                                 design.designVariable[i].upperBound[0]);
                break;
            }

            if (i >= design.numDesignVariable) fprintf(fp, "xrate   0   0.000000000000000E+00  0.000000000000000E+00  0.00000000000000E+00\n");

            for (i = 0; i < design.numDesignVariable; i++) {

                if (strcasecmp(design.designVariable[i].name, "NonInertial_Rotation_Rate") != 0) continue;

                if (design.designVariable[i].var->length < 3) {
                    status = CAPS_RANGEERR;
                    goto cleanup;
                }

                fprintf(fp, "yrate   1   %.15E  %.15E  %.15E\n", design.designVariable[i].value[1],
                                                                 design.designVariable[i].lowerBound[1],
                                                                 design.designVariable[i].upperBound[1]);
                break;
            }

            if (i >= design.numDesignVariable) fprintf(fp, "yrate   0   0.000000000000000E+00  0.000000000000000E+00  0.00000000000000E+00\n");

            for (i = 0; i < design.numDesignVariable; i++) {

                if (strcasecmp(design.designVariable[i].name, "NonInertial_Rotation_Rate") != 0) continue;

                if (design.designVariable[i].var->length < 3) {
                    status = CAPS_RANGEERR;
                    goto cleanup;
                }

                fprintf(fp, "zrate   1   %.15E  %.15E  %.15E\n", design.designVariable[i].value[2],
                                                                 design.designVariable[i].lowerBound[2],
                                                                 design.designVariable[i].upperBound[2]);
                break;
            }

            if (i >= design.numDesignVariable) fprintf(fp, "zrate   0   0.000000000000000E+00  0.000000000000000E+00  0.00000000000000E+00\n");
        }

        fprintf(fp, "Number of bodies\n");
        fprintf(fp, "%d\n", numBody);
        for (i = 0; i < numBody; i++) {
            fprintf(fp, "Rigid motion design variables for 'Body %d'\n", i+1);
            fprintf(fp, "Var Active         Value               Lower Bound            Upper Bound\n");
            fprintf(fp, "RotRate  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00\n");
            fprintf(fp, "RotFreq  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00\n");
            fprintf(fp, "RotAmpl  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00\n");
            fprintf(fp, "RotOrgx  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00\n");
            fprintf(fp, "RotOrgy  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00\n");
            fprintf(fp, "RotOrgz  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00\n");
            fprintf(fp, "RotVecx  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00\n");
            fprintf(fp, "RotVecy  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00\n");
            fprintf(fp, "RotVecz  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00\n");
            fprintf(fp, "TrnRate  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00\n");
            fprintf(fp, "TrnFreq  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00\n");
            fprintf(fp, "TrnAmpl  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00\n");
            fprintf(fp, "TrnVecx  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00\n");
            fprintf(fp, "TrnVecy  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00\n");
            fprintf(fp, "TrnVecz  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00\n");
            if (fun3dVersion > 12.4) {
                // FUN3D version 13.1 - version 12.4 doesn't have option 5
                fprintf(fp, "Parameterization Scheme (Massoud=1 Bandaids=2 Sculptor=4 User-Defined=5)\n");
                fprintf(fp, "5\n");
            } else {

                fprintf(fp, "Parameterization Scheme (Massoud=1 Bandaids=2 Sculptor=4)\n");
                fprintf(fp, "1\n");
            }

            fprintf(fp, "Number of shape variables for 'Body %d'\n", i+1);
            fprintf(fp, "%d\n", numShapeVar);
            fprintf(fp, "Index Active         Value               Lower Bound            Upper Bound\n");

            m = 1;
            for (j = 0; j < design.numDesignVariable; j++) {

                if (design.designVariable[j].type != DesignVariableGeometry) continue;

                for (k = 0; k < design.designVariable[j].var->length; k++ ) {

                    fprintf(fp, "%d    1   %.15E  %.15E  %.15E\n", m, design.designVariable[j].value[k],
                                                                      design.designVariable[j].lowerBound[k],
                                                                      design.designVariable[j].upperBound[k]);

                    m += 1;
                }
            }
        }

        fprintf(fp, "################################################################################\n");
        fprintf(fp, "############################### Function Information ###########################\n");
        fprintf(fp, "################################################################################\n");
        fprintf(fp, "Number of composite functions for design problem statement\n");
        fprintf(fp, "%d\n", design.numDesignFunctional);

        for (i = 0; i < design.numDesignFunctional; i++) {
            fprintf(fp, "################################################################################\n");
            fprintf(fp, "Cost function (1) or constraint (2)\n");
            fprintf(fp, "1\n");
            fprintf(fp, "If constraint, lower and upper bounds\n");
            fprintf(fp, "0.0 0.0\n");
            fprintf(fp, "Number of components for function   %d\n", i+1);
            fprintf(fp, "%d\n", design.designFunctional[i].numComponent);
            fprintf(fp, "Physical timestep interval where function is defined\n");
            fprintf(fp, "1 1\n");
            fprintf(fp, "Composite function weight, target, and power\n");
            fprintf(fp, "1.0 0.0 1.0\n");
            fprintf(fp, "Components of function   %d: boundary id (0=all)/name/value/weight/target/power\n", i+1);

            for (j = 0; j < design.designFunctional[i].numComponent; j++) {
                //fprintf(fp, "0 clcd          0.000000000000000    1.000   10.00000 2.000\n");
                status = _writeFunctinoalComponent(aimInfo, fp, &design.designFunctional[i].component[j]);
                if (status != CAPS_SUCCESS) goto cleanup;
            }

            fprintf(fp, "Current value of function   %d\n", i+1);
            fprintf(fp, "0.000000000000000\n");
            fprintf(fp, "Current derivatives of function wrt global design variables\n");
            fprintf(fp, "0.000000000000000\n");
            fprintf(fp, "0.000000000000000\n");

            if (fun3dVersion > 12.4) {
                // FUN3D version 13.1 - version 12.4 doesn't have these available
                fprintf(fp, "0.000000000000000\n"); // Yaw
                fprintf(fp, "0.000000000000000\n"); // xrate
                fprintf(fp, "0.000000000000000\n"); // yrate
                fprintf(fp, "0.000000000000000\n"); // zrate
            }

            for (j = 0; j < numBody; j++) {
                fprintf(fp, "Current derivatives of function wrt rigid motion design variables of body %d\n", j+1);
                fprintf(fp, "0.000000000000000\n");
                fprintf(fp, "0.000000000000000\n");
                fprintf(fp, "0.000000000000000\n");
                fprintf(fp, "0.000000000000000\n");
                fprintf(fp, "0.000000000000000\n");
                fprintf(fp, "0.000000000000000\n");
                fprintf(fp, "0.000000000000000\n");
                fprintf(fp, "0.000000000000000\n");
                fprintf(fp, "0.000000000000000\n");
                fprintf(fp, "0.000000000000000\n");
                fprintf(fp, "0.000000000000000\n");
                fprintf(fp, "0.000000000000000\n");
                fprintf(fp, "0.000000000000000\n");
                fprintf(fp, "0.000000000000000\n");
                fprintf(fp, "0.000000000000000\n");
                fprintf(fp, "Current derivatives of function wrt shape design variables of body %d\n", j+1);

                for (k = 0; k < design.numDesignVariable; k++) {

                    if (design.designVariable[k].type != DesignVariableGeometry) continue;

                    for (m = 0; m < design.designVariable[k].var->length; m++ ) fprintf(fp, "0.000000000000000\n");
                }
            }
        }

        status = CAPS_SUCCESS;

    cleanup:
        if (status != CAPS_SUCCESS)
            printf("Error: Premature exit in fun3d_writeRubber status = %d\n",
                   status);

        if (fp != NULL) fclose(fp);

        return status;

        /* Version 13.1 - template
    ################################################################################
    ########################### Design Variable Information ########################
    ################################################################################
    Global design variables (Mach number, AOA, Yaw, Noninertial rates)
    Var Active         Value               Lower Bound            Upper Bound
    Mach    0   0.800000000000000E+00  0.000000000000000E+00  0.900000000000000E+00
      AOA    1   1.000000000000000E+00  0.000000000000000E+00  5.000000000000000E+00
      Yaw    0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    xrate    0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    yrate    0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    zrate    0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    Number of bodies
    2
    Rigid motion design variables for 'wing'
    Var Active         Value               Lower Bound            Upper Bound
    RotRate  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    RotFreq  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    RotAmpl  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    RotOrgx  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    RotOrgy  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    RotOrgz  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    RotVecx  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    RotVecy  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    RotVecz  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    TrnRate  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    TrnFreq  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    TrnAmpl  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    TrnVecx  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    TrnVecy  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    TrnVecz  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    Parameterization Scheme (Massoud=1 Bandaids=2 Sculptor=4 User-Defined=5)
    1
    Number of shape variables for 'wing'
    3
    Index Active         Value               Lower Bound            Upper Bound
    1    1   1.000000000000000E+00  0.000000000000000E+00  2.000000000000000E+00
    2    1   1.000000000000000E+00  0.000000000000000E+00  2.000000000000000E+00
    3    1   1.000000000000000E+00  0.000000000000000E+00  2.000000000000000E+00
    Rigid motion design variables for 'tail'
    Var Active         Value               Lower Bound            Upper Bound
    RotRate  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    RotFreq  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    RotAmpl  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    RotOrgx  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    RotOrgy  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    RotOrgz  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    RotVecx  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    RotVecy  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    RotVecz  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    TrnRate  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    TrnFreq  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    TrnAmpl  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    TrnVecx  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    TrnVecy  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    TrnVecz  0   0.000000000000000E+00  0.000000000000000E+00  0.000000000000000E+00
    Parameterization Scheme (Massoud=1 Bandaids=2 Sculptor=4 User-Defined=5)
    2
    Number of shape variables for 'tail'
    2
    Index Active         Value               Lower Bound            Upper Bound
    1    1   2.000000000000000E+00 -1.000000000000000E+00  5.000000000000000E+00
    2    1   2.000000000000000E+00 -1.000000000000000E+00  5.000000000000000E+00
    ################################################################################
    ############################### Function Information ###########################
    ################################################################################
    Number of composite functions for design problem statement
    2
    ################################################################################
    Cost function (1) or constraint (2)
    1
    If constraint, lower and upper bounds
    0.0 0.0
    Number of components for function   1
    1
    Physical timestep interval where function is defined
    1 1
    Composite function weight, target, and power
    1.0 0.0 1.0
    Components of function   1: boundary id (0=all)/name/value/weight/target/power
    0 cl            0.000000000000000    1.000   10.00000 2.000
    Current value of function   1
    0.000000000000000
    Current derivatives of function wrt global design variables
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    Current derivatives of function wrt rigid motion design variables of body   1
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    Current derivatives of function wrt shape design variables of body   1
    0.000000000000000
    0.000000000000000
    0.000000000000000
    Current derivatives of function wrt rigid motion design variables of body   2
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    Current derivatives of function wrt shape design variables of body   2
    0.000000000000000
    0.000000000000000
    ################################################################################
    Cost function (1) or constraint (2)
    2
    If constraint, lower and upper bounds
    -0.03 -0.01
    Number of components for function   2
    1
    Physical timestep interval where function is defined
    1 1
    Composite function weight, target, and power
    1.0 0.0 1.0
    Components of function   2: boundary id (0=all)/name/value/weight/target/power
    0 cmy           0.000000000000000    1.000    0.00000 1.000
    Current value of function   2
    0.000000000000000
    Current derivatives of function wrt global design variables
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    Current derivatives of function wrt rigid motion design variables of body   1
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    Current derivatives of function wrt shape design variables of body   1
    0.000000000000000
    0.000000000000000
    0.000000000000000
    Current derivatives of function wrt rigid motion design variables of body   2
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    0.000000000000000
    Current derivatives of function wrt shape design variables of body   2
    0.000000000000000
    0.000000000000000
         */
    }

    // Finds a line in rubber.data that matches header
    static int findHeader(const char *header, char **line, size_t *nline, int *iline, FILE *fp)
    {
        int i;
        while (getline(line, nline, fp) != -1) {
            (*iline)++;
            if (nline == 0) continue;
            if ((*line)[0] == '#') continue;
            i = 0;
            while ((*line)[i] == ' ' && (*line)[i] != '\0') i++;
            if (strncasecmp(header, &(*line)[i], strlen(header)) == 0)
              break;
        }
        if (feof(fp)) return CAPS_IOERR;

        return CAPS_SUCCESS;
    }


    // Read objective value and derivatives from FUN3D rubber.data file
    // Will not read shape entries unless explicitly told to check if they are needed
    int fun3d_readRubber(void *aimInfo,
                         cfdDesignStruct design,
                         int checkGeomShape,
                         double fun3dVersion,
                         aimMeshRef *meshRef)
    {
        int status; // Function return status

        int i, j, k, ibody; // Indexing

        int numShapeVar = 0;

        size_t nline;
        int iline=0;
        char *line=NULL;
        char file[] = "rubber.data";
        double dfun_dvar = 0;
        FILE *fp = NULL;

        int numBody;

        // Get AIM bodies
        numBody = meshRef->nmap;

        printf("Reading %s \n", file);

        // Open file
        fp = aim_fopen(aimInfo, file, "r");
        if (fp == NULL) {
            AIM_ERROR(aimInfo, "Unable to open file: %s\n", file);
            status = CAPS_IOERR;
            goto cleanup;
        }

        for (i = 0; i < design.numDesignFunctional; i++) {

            status = findHeader("Current value of function", &line, &nline, &iline, fp);
            AIM_STATUS(aimInfo, status, "rubber.data line %d", iline);

            // Get the line with the objective value
            status = getline(&line, &nline, fp); iline++;
            if (status == -1) status = CAPS_IOERR; else status = CAPS_SUCCESS;
            AIM_STATUS(aimInfo, status, "rubber.data line %d", iline);
            AIM_NOTNULL(line, aimInfo, status);

            status = sscanf(line, "%lf", &design.designFunctional[i].value);
            if (status != 1) status = CAPS_IOERR; else status = CAPS_SUCCESS;
            AIM_STATUS(aimInfo, status, "rubber.data line %d", iline);

            // skip: Current derivatives of function wrt rigid motion design variables of
            status = getline(&line, &nline, fp); iline++;
            if (status == -1) status = CAPS_IOERR; else status = CAPS_SUCCESS;
            AIM_STATUS(aimInfo, status, "rubber.data line %d", iline);
            AIM_NOTNULL(line, aimInfo, status);

            // read the dFun/dMach
            status = getline(&line, &nline, fp); iline++;
            if (status == -1) status = CAPS_IOERR; else status = CAPS_SUCCESS;
            AIM_STATUS(aimInfo, status, "rubber.data line %d", iline);
            AIM_NOTNULL(line, aimInfo, status);

            for (j = 0; j < design.designFunctional[i].numDesignVariable; j++) {

                if (strcasecmp(design.designFunctional[i].dvar[j].name, "Mach") != 0) continue;

                status = sscanf(line, "%lf", &design.designFunctional[i].dvar[j].value[0]);
                if (status != 1) status = CAPS_IOERR; else status = CAPS_SUCCESS;
                AIM_STATUS(aimInfo, status, "rubber.data line %d", iline);
                break;
            }

            // read the dFun/dAOA
            status = getline(&line, &nline, fp); iline++;
            if (status == -1) status = CAPS_IOERR; else status = CAPS_SUCCESS;
            AIM_STATUS(aimInfo, status, "rubber.data line %d", iline);
            AIM_NOTNULL(line, aimInfo, status);

            for (j = 0; j < design.numDesignVariable; j++) {

                if (strcasecmp(design.designFunctional[i].dvar[j].name, "Alpha") != 0) continue;

                status = sscanf(line, "%lf", &design.designFunctional[i].dvar[j].value[0]);
                if (status != 1) status = CAPS_IOERR; else status = CAPS_SUCCESS;
                AIM_STATUS(aimInfo, status, "rubber.data line %d", iline);
                break;
            }

            if (fun3dVersion > 12.4) {
                // FUN3D version 13.1 - version 12.4 doesn't have these available

                // read the dFun/dBeta
                status = getline(&line, &nline, fp); iline++;
                if (status == -1) status = CAPS_IOERR; else status = CAPS_SUCCESS;
                AIM_STATUS(aimInfo, status, "rubber.data line %d", iline);
                AIM_NOTNULL(line, aimInfo, status);

                for (j = 0; j < design.designFunctional[i].numDesignVariable; j++) {

                    if (strcasecmp(design.designFunctional[i].dvar[j].name, "Beta") != 0) continue;

                    status = sscanf(line, "%lf", &design.designFunctional[i].dvar[j].value[0]);
                    if (status != 1) status = CAPS_IOERR; else status = CAPS_SUCCESS;
                    AIM_STATUS(aimInfo, status, "rubber.data line %d", iline);
                    break;
                }

                // read the dFun/dNonInertial_Rotation_Rate[0]
                status = getline(&line, &nline, fp); iline++;
                if (status == -1) status = CAPS_IOERR; else status = CAPS_SUCCESS;
                AIM_STATUS(aimInfo, status, "rubber.data line %d", iline);
                AIM_NOTNULL(line, aimInfo, status);

                for (j = 0; j < design.designFunctional[i].numDesignVariable; j++) {

                    if (strcasecmp(design.designFunctional[i].dvar[j].name, "NonInertial_Rotation_Rate") != 0) continue;

                    status = sscanf(line, "%lf", &design.designFunctional[i].dvar[j].value[0]);
                    if (status != 1) status = CAPS_IOERR; else status = CAPS_SUCCESS;
                    AIM_STATUS(aimInfo, status, "rubber.data line %d", iline);
                    break;
                }

                // read the dFun/dNonInertial_Rotation_Rate[1]
                status = getline(&line, &nline, fp); iline++;
                if (status == -1) status = CAPS_IOERR; else status = CAPS_SUCCESS;
                AIM_STATUS(aimInfo, status, "rubber.data line %d", iline);
                AIM_NOTNULL(line, aimInfo, status);

                for (j = 0; j < design.designFunctional[i].numDesignVariable; j++) {

                    if (strcasecmp(design.designFunctional[i].dvar[j].name, "NonInertial_Rotation_Rate") != 0) continue;

                    status = sscanf(line, "%lf", &design.designFunctional[i].dvar[j].value[1]);
                    if (status != 1) status = CAPS_IOERR; else status = CAPS_SUCCESS;
                    AIM_STATUS(aimInfo, status, "rubber.data line %d", iline);
                    break;
                }

                // read the dFun/dNonInertial_Rotation_Rate[2]
                status = getline(&line, &nline, fp); iline++;
                if (status == -1) status = CAPS_IOERR; else status = CAPS_SUCCESS;
                AIM_STATUS(aimInfo, status, "rubber.data line %d", iline);
                AIM_NOTNULL(line, aimInfo, status);

                for (j = 0; j < design.designFunctional[i].numDesignVariable; j++) {

                    if (strcasecmp(design.designFunctional[i].dvar[j].name, "NonInertial_Rotation_Rate") != 0) continue;

                    status = sscanf(line, "%lf", &design.designFunctional[i].dvar[j].value[2]);
                    if (status != 1) status = CAPS_IOERR; else status = CAPS_SUCCESS;
                    AIM_STATUS(aimInfo, status, "rubber.data line %d", iline);
                    break;
                }
            }


            numShapeVar = 0; // Need to re-zero every functional loop

            if (checkGeomShape == (int) true) {
                // zero out previous derivatives
                for (j = 0; j < design.designFunctional[i].numDesignVariable; j++) {
                    if (design.designFunctional[i].dvar[j].type != DesignVariableGeometry) continue;
                    numShapeVar++;
                    for (k = 0; k < design.designVariable[j].var->length; k++ ) {
                      design.designFunctional[i].dvar[j].value[k] = 0;
                    }
                }
            }

            // No body information if there are no Geometry derivatives
            if (numShapeVar == 0) numBody = 0;

            for (ibody = 0; ibody < numBody; ibody++) {

              // Rigid motion design variables
              status = findHeader("Current derivatives of", &line, &nline, &iline, fp);
              AIM_STATUS(aimInfo, status, "rubber.data line %d", iline);

              // Skip reading rigid motion design variables for now
             
              // Read shape design variables
              status = findHeader("Current derivatives of", &line, &nline, &iline, fp);
              AIM_STATUS(aimInfo, status, "rubber.data line %d", iline);

              for (j = 0; j < design.designFunctional[i].numDesignVariable; j++) {

                  if (design.designFunctional[i].dvar[j].type != DesignVariableGeometry) continue;

                  for (k = 0; k < design.designVariable[j].var->length; k++ ) {

                      // read the dFun/dDesignVar[k]
                      status = getline(&line, &nline, fp); iline++;
                      if (status == -1) status = CAPS_IOERR; else status = CAPS_SUCCESS;
                      AIM_STATUS(aimInfo, status, "rubber.data line %d", iline);
                      AIM_NOTNULL(line, aimInfo, status);

                      status = sscanf(line, "%lf", &dfun_dvar);
                      if (status != 1) status = CAPS_IOERR; else status = CAPS_SUCCESS;
                      AIM_STATUS(aimInfo, status, "rubber.data line %d", iline);

                      // accumulate derivatives across bodies
                      design.designFunctional[i].dvar[j].value[k] += dfun_dvar;
                  }
              }
            }
        }

        status = CAPS_SUCCESS;

    cleanup:

        if (line != NULL) free(line); // Must use free!
        if (fp   != NULL) fclose(fp);

        return status;
    }

    // Make FUN3D directory structure/tree
    int fun3d_makeDirectory(void *aimInfo)
    {
        int status; // Function return status

        char filename[PATH_MAX];

        char flow[] = "Flow";
        char datafile[] = "datafiles";
        char adjoint[] = "Adjoint";
        char rubber[] = "Rubberize";

        printf("Creating FUN3D directory tree\n");

        // Flow
        status = aim_mkDir(aimInfo, flow);
        AIM_STATUS(aimInfo, status);

        // Datafiles
        sprintf(filename, "%s/%s", flow, datafile);
        status = aim_mkDir(aimInfo, filename);
        AIM_STATUS(aimInfo, status);

        // Adjoint
        status = aim_mkDir(aimInfo, adjoint);
        AIM_STATUS(aimInfo, status);

        // Rubber
        status = aim_mkDir(aimInfo, rubber);
        AIM_STATUS(aimInfo, status);

        status = CAPS_SUCCESS;

    cleanup:
        return status;
    }



    // Morphing input
    // First time through and morphing
    // -  Write bodies and mesh to disc only if mesh is linked a
    // -  AIM_CP_FIle the meshfile
    //


    // aim_storeMeshRef(aimInfo, meshRef, MESHEXTENSION);
    //aim_storeMeshRef(aimInfo, meshRef, MESHEXTENSION);
    //
    //aimMeshRef meshRef2;
    //aim_initMeshRef(&meshRef2);
    //aim_loadMeshRef(aimInfo, &meshRef2);
    //aim_freeMeshRef(&meshRef2);

    // Logic for selecting/setting mesh pointer if morphing the mesh
    int fun3d_morphMeshLogic(void *aimInfo, int storeAllowed, aimMeshRef **meshRef, aimMeshRef *meshRefObj)
    {
        int status = CAPS_SUCCESS;

        AIM_NOTNULL(meshRefObj, aimInfo, status);

        if (*meshRef != NULL && storeAllowed == (int) true) { // If we are NOT null lets store the mesh - if allowed

            status = aim_storeMeshRef(aimInfo, *meshRef, MESHEXTENSION);
            AIM_STATUS(aimInfo, status);

            *meshRef =  meshRefObj;

           // Lets "reload" the meshRef now - A little awkward but.....
           status = aim_loadMeshRef(aimInfo, *meshRef);
           AIM_STATUS(aimInfo, status);

        } else if (*meshRef != NULL  && storeAllowed == (int) false) { // If we are NOT null

             // Do nothing

        } else { // *meshRef == NULL

            *meshRef =  meshRefObj;

            // Lets "load" the meshRef now
            status = aim_loadMeshRef(aimInfo, *meshRef);
            AIM_STATUS(aimInfo, status);
        }

        cleanup:
            return status;
    }


    //{
    //    int status = CAPS_SUCCESS;
    //
    //    AIM_NOTNULL(meshRefObj, aimInfo, status);
    //
    //    if (*meshRef != NULL && storeAllowed == (int) true) { // If we are NOT null lets store the mesh
    //        status = aim_storeMeshRef(aimInfo, *meshRef, MESHEXTENSION);
    //        AIM_STATUS(aimInfo, status);
    //
    //        meshRef = NULL;
    //    }
    //
    //    *meshRef =  meshRefObj;
    //
    //    // Lets "reload" the meshRef now - A little awkward but.....
    //    status = aim_loadMeshRef(aimInfo, *meshRef);
    //    AIM_STATUS(aimInfo, status);
    //
    //    cleanup:
    //        return status;
    //}

    // Map the old surface tessellation on to the new bodies
    int fun3d_morpMeshUpdate(void *aimInfo,  aimMeshRef *meshRef, int numBody, ego *bodies)
    {
        int status = CAPS_SUCCESS;
        int i=0;
        int state, numVert;
        ego bodyMapping;
        ego body;
        ego tess;

        // Have the number of bodies changed?
        if (meshRef->nmap != numBody) {
            printf("The number of original surface meshes does NOT equal the number of current bodies!\n");
            AIM_STATUS(aimInfo, CAPS_MISMATCH);
        }

        // Are the bodies topological equivalent?
        for (i = 0; i < meshRef->nmap; i++) {

            status = EG_statusTessBody(meshRef->maps[i].tess, &body, &state, &numVert);
            AIM_STATUS(aimInfo, status);

            status = EG_mapBody(body, bodies[i], "_faceID",  &bodyMapping); // "_faceID" - same as in OpenCSM
            if (status != EGADS_SUCCESS || bodyMapping != NULL) {

                printf("New and old body %d (of %d) do not appear to be topologically equivalent!\n", i+1, meshRef->nmap);

                if (bodyMapping != NULL) {
                    printf("Body mapping isn't NULL!\n");
                }

                AIM_STATUS(aimInfo, CAPS_MISMATCH);
            }
        }

        // Now lets "tweak" the surface tessellation - map the old tessellation to the new bodies
        for (i = 0; i < meshRef->nmap; i++) {

            status = EG_copyObject(bodies[i], NULL, &body);
            AIM_STATUS(aimInfo, status);

            status = EG_mapTessBody(meshRef->maps[i].tess,
                                    body,
                                    &tess);
            AIM_STATUS(aimInfo, status);

            status = EG_deleteObject(meshRef->maps[i].tess);
            meshRef->maps[i].tess = tess;

        }

        cleanup:
            return status;
    }
    ryan / detail
  50. Remove old egads file in aim_storeMeshRef. Set AIMptr to NULL when unlinking. — galbramc / detail
  51. clarify description of SELECT x -1 ibody1 — jfdannen / detail
  52. update to account for latest EGADS fix to EG_generalBoolean — jfdannen / detail
  53. add ocsmRegBcstCB to ocsm.def — jfdannen / detail
  54. update Building button to show progress of build in ESP; fix memory leak when ELEVATE command fails (like it often does in the combine.udc); extend SELECT command to SELECT all Faces/Edges/Nodes in current Body that match a Face in another Body; add SELECT NOT to SELECT the Faces/Edges/Nodes not selected and vice versa; fix bug associated with determining if an Edge is a subset of another; update to account for change in EG_isEquivalent; fix bug that did not copy attributes from an Edge in a WireBody that was INTERSECTed with a SolidBody to the new Edge; add PRINT_LAST_EGO compile-time flag to print the last ego after every step in an ocsmBuild; fix bug in sensCSM when trying to remove .csm or .cpc extension from the casename when it did not exist; extend SELECT command to SELECT all Faces/Edges/Nodes in current Body that are within another Body — jfdannen / detail
  55. fix DUMP x.sens for degenerate Edges; add combine.udc so that old .csm scripts that used the (now deprecated) COMBINE command still work; add combine* test cases; add makeErep journal entry so that ereped can be tested; modify at-itype, at-nface, and at-nedge if there is an erep; add POWLE and POWTE arguments to udfDroop; add new droop* test cases; add udfNaca6mc to warp the input thickness distribution with the sum of multiple NACA 6-series meanlines; add naca6mc* test cases — jfdannen / detail
  56. protect against string overflows in strncpy and strncat; fix mesgCallbackFromOpenCSM so that it does not try to broadcast text when in batch mode; initial version of udpNaca6mc — jfdannen / detail
  57. Remove dependencies from file copies — haimes / detail
  58. clean up usage of include files; update dependencies in Makefile and NMakefile — jfdannen / detail
  59. remove pyscript from lint in serveESP.make — jfdannen / detail
  60. update serveCSM.make to fix stanilizer error for pyscript — jfdannen / detail
  61. another attempt to move building of pyscript from Makefile to serveESP.make — jfdannen / detail
  62. move building of pyscript from Makefile to serveESP.make — jfdannen / detail
  63. fix recently added bug in ereped that only allowed 9 colors; remove timHold from viewBodys.py; add batch flag to ESP structure; add assert command in browserToServer; add method of verifying ereped and plugs if running from journal and in -batch mode; next versions of udfDroop2 and udfDroop3; migrate testing in Makefile.DARWIN64 and Makefile.LINUX64 from serveCSM to serveESP; fix bug in esp.SetCaps() in pyESP — jfdannen / detail
  64. next version of udfDroop2; first version of udfDroop3 — jfdannen / detail
  65. next (still incomplete) version of udfDroop2 — jfdannen / detail
  66. initial version of udfDroop2 — jfdannen / detail
  67. fix test_pyOCSM to account for updated MESSAGE statement — jfdannen / detail
  68. extend MESSAGE to optionally write to a file; modify message test case; update ESP-help and ESP_QuickReference — jfdannen / detail
  69. fix memory leak when a UDF returns an error; remove extra prints when udfFlend or udfOffset detect an error; update sanitize targets in Makefile.LINUX — jfdannen / detail
  70. allow OFFSET to function properly when the number of Edges in the offset is fewer than the original Body; fix escape-close-parenthesis in ESP_QuickReference; add ability to create periodic SPLINE in sketcher by using SSLOPE(0); add sslope6* test cases — jfdannen / detail
  71. add assert in udfOffset to eliminate warnings on gcc — jfdannen / detail
  72. fix bug that caused ereped to not work for a named Body; fix bug that added bar at end of file that was being editted in ESP; extend ereped so that it can be applied if there are multiple Bodys on the stack; add ereped8 test case; fix ESP hint associated with ELEVATE command; fix bug when SELECT EGDE with a bounding box was applied to a Body with degenerate Edges; add extractBodys test case; fix array bounds error in udpEqn2body; fix bug in ereped initializations for a single Body — jfdannen / detail
  73. remove use of _scaleuv Attribute — jfdannen / detail
  74. allow very long message to be sent from browser to server; add setCsmFile instead of setScmFileBeg/Mid/End; if error occurs when processing a timMesg, put the tim back into READY mode; post message when erep is successfully generated; implement udpEqn2body for SheetBodys; update frustrum3 test case; use __usedBodys__ attribute when computing velocity of Faces or Edges for RULE and BLEND; fix bug in sensCSM that did not reset maximum errors for Bodys when in -tess mode; fix big that did not copy .hasdots during RESTORE . (dot) — jfdannen / detail
  75. Fix the new isEquivalent for periodic curves split into 2 Edges — haimes / detail
  76. modified Ints to Int32 and added egadslite tests — docampo / detail
  77. recompiled pdf references were with ?? — docampo / detail
  78. removed EL flag from some of the effectiveTopology functions — docampo / detail
  79. fixed typos — docampo / detail
  80. Remove bug that would error out if there were no Faces during generalBoolean operations — haimes / detail
  81. Update isEquivalent to fix previous problem — haimes / detail
  82. Allow for senses to be flipped when applying EG_isEquivalent — haimes / detail
  83. added the ego. wrap to the getfield operation — docampo / detail
  84. update udpTire for new include scheme — jfdannen / detail
  85. Fix Edge mapping problem introduced since Rev 1.21 — haimes / detail
  86. changed .obj -> .ego — docampo / detail
  87. Add TOML to egadslite dep — galbramc / detail
  88. compiling for Julia 1.6 and 1.8. test fails for 1.8 — docampo / detail
  89. Update makeEnv to reflect newer MAC OS revs — haimes / detail
  90. Add the Julia interface to the EGADS API documentation — haimes / detail
  91. Try to fix stanalyzer errors again — haimes / detail
  92. Fix jlEGADS makefiles — galbramc / detail
  93. Try to fix stanalyzer errors — haimes / detail
  94. Suppress gcc-12 warning — galbramc / detail
  95. Effective Topology: apply last change only for EFaces with more than 1 Face — haimes / detail
  96. Effective Topology: check result of getEdgeUV against invEval and pick the closer — haimes / detail
  97. A different possible fix for the MemcheckOcsm on the rule28 cases — haimes / detail
  98. A different possible fix for the MemcheckOcsm on the rule28 cases — haimes / detail
  99. A possible fix for the MemcheckOcsm on the rule28 cases — haimes / detail
  100. Effective Topology: Mark Planar/Periodic EFaces not to tessellate Edges based on Face Curvature — haimes / detail
  101. Fix scan-build warning from last commit again — haimes / detail
  102. Fix scan-build warning from last commit — haimes / detail
  103. Effective Topology: Fix NM Sheetbody Node removal — haimes / detail
  104. Fix NULL bays at end with multiNode Rule — haimes / detail
  105. Updated ruled sensitvities for multiNode — galbramc / detail
  106. A minor update of the documentation of EG_approximate — haimes / detail
  107. Update the docs to specify the additions to EG_approximate and EG_ruled — haimes / detail
  108. Suppress scipy memory leak — galbramc / detail
  109. Another python valgrind suppression — galbramc / detail
  110. Add Ubuntu 22.04 valgrind suppressions — galbramc / detail
  111. Use system valgrind — galbramc / detail
  112. Small cleanup of email template — galbramc / detail
  113. Remove torch for now — galbramc / detail
  114. Log parser fixes — galbramc / detail
  115. Suppress torch memory leaks — galbramc / detail
  116. Upgrade sanitizer path to 14.0.0 — galbramc / detail
  117. Try dissabling CUDA for pytorch — galbramc / detail
  118. Look for fatal msvc errors — galbramc / detail
  119. Look for fatal gnu/clang errors — galbramc / detail
  120. Revert to serveCSM for plugs — galbramc / detail
  121. Update lingering serveCSM to serveESP — galbramc / detail
  122. Very strange... — galbramc / detail
  123. Allow multiple csm patterns — galbramc / detail
  124. Only check error handling when running all csm files — galbramc / detail
  125. Support changing outLevel and running select files — galbramc / detail
  126. Don't do jlEGADS testing just yet — galbramc / detail
  127. Revert ESP_venv for ESP_Beta — galbramc / detail
  128. Move virtualenv into ESP directory — galbramc / detail
  129. Only run minimal CAPS examples in beta — galbramc / detail
  130. Install torch for beta — galbramc / detail
  131. pip install torch for corsair — galbramc / detail
  132. Smartsr error log parse — galbramc / detail

#166 (Aug 13, 2022, 6:19:59 AM)

  1. Suppress lint warnings — galbramc / detail
  2. Fix for last commit — galbramc / detail
  3. Add aim_storeMeshRef and aim_loadMeshRef — galbramc / detail
  4. Doc updates — haimes / detail
  5. Fix typo in doc, and short Python Units Appendix — galbramc / detail
  6. Minor documentation (and coding) changes for the Python 1-to-1 CAPS interface — haimes / detail
  7. Document 1-to-1 pyCAPS API. Change pyCAPS setLimits/getLimits to setLimitsSize/getLimitsSize. Add support for limit arrays to pyCAPS. — galbramc / detail
  8. Addition of a CAPS API function that anly returns a Value Object's size — haimes / detail
  9. Try again! — haimes / detail
  10. Only set limits on OpenCSM DESPMTRs — haimes / detail
  11. Completed the first cut for per element support of limits — haimes / detail
  12. Adjust Value functions to support per element limits — haimes / detail
  13. Update bfr instructions — galbramc / detail
  14. Fix lint error — haimes / detail
  15. readd multifidelity phasing examples — dongjoon / detail
  16. Read/write and update setLimits/getLimits for per element limits — haimes / detail
  17. Skeleton additions for adding per element limits — haimes / detail
  18. Support for higher precssion xfoil executables. Add BoundaryLayerIntersect meshProp for AFLR4 and AFLR3. — galbramc / detail
  19. Compile needed aflr43 files — galbramc / detail
  20. Upgrade to AFLR 10.22.23 — galbramc / detail
  21. Uprev CAPS — haimes / detail
  22. Update corasairlite README to indicate hoburg_blackbox.py requies scipy. Remove hoburg_lsqp.py from testing. — galbramc / detail
  23. delete duplicate file — dongjoon / detail
  24. update Makefile.LINUX64 — jfdannen / detail
  25. add udfDumpPmtrs and associated test case; implement RESTORE .. (dot dot) to duplicate all Bodys (and Mark) back to the Mark and RESTORE ... (dot dot dot) to duplicate entire stack; add restore3 test case; allow UDF to not return a Body or Bodys; update udfOffset so that it can handle concave corners (but not case where Edges disappear) — jfdannen / detail
  26. update udfOffset for negative distance — jfdannen / detail
  27. fix hint associated with JOIN; change JOIN so that it can make non-manifold WireBodys; modify numWireBody* test cases; create ELEVATE command, which only elevates WireBody to SheetBody and SheetBody to SolidBody;  remove COMBINE command and instruct user to use JOIN and ELEVATE instead; update UNION and JOIN such that if toMark is set but the operation fails, the stack is left unaltered; add unionSolids* test cases; fix typos in ESP-help (and OpenCSM.h) — jfdannen / detail
  28. add (unneeded) initializations in udfOffset to resolve incorrect lint complaints — jfdannen / detail
  29. implement udpEqn2body for WireBodys; modify inflatedBody2b to use the replicate.udc; implement udfOffset; create offset* test cases; improve error handling associated with SELECT SORT; modify flend5* test cases to use udfOffset — jfdannen / detail
  30. uprev to v1.22 — jfdannen / detail
  31. fix bug that did not update at-parameters for REORDER command; fix bug that changed the input Body to REORDER; fix sensitivity bug whan a RULE/BLEND xsect is RESTOREd after being used — jfdannen / detail
  32. initialize *info in ocsmBodyDetails; return sensible error message when a SET is applied to a CONPMTR, CFGPMTR, or DESPMTR — jfdannen / detail
  33. Fix _close! function to properly check garbage collection for both egads and egadslite — galbramc / detail
  34. Fix jlEGADS makefiles for zsh — galbramc / detail
  35. jlEGADS makefile adjusments. Other small cleanup. — galbramc / detail
  36. delete manifest from lite — docampo / detail
  37. now ?? — docampo / detail
  38. egads also used from egadscommon.jl dependencies — docampo / detail
  39. remove unnecesary ENVS — docampo / detail
  40. now with actually things inside — docampo / detail
  41. created egads common folders — docampo / detail
  42. Rule with Loops of differing number of Edges passes simple test with null 'bay' — haimes / detail
  43. Can rule with sections of different number of Edges -- still need to check null bays, multiplicity > 1, curved EndCaps & EG_splineGeom_dot — haimes / detail
  44. Close memory leak in EGADSlite associated with full attributes — haimes / detail
  45. Fix scan-build problem — haimes / detail
  46. Fix warning from last commit — haimes / detail
  47. Second (but incomplete) cut at Rule with differring numbers of Edges — haimes / detail
  48. Another minor correction to last commit — haimes / detail
  49. Minor correction to last commit — haimes / detail
  50. First (but incomplete) cut at Rule with differring numbers of Edges — haimes / detail
  51. Don't allow for Figure 8 Nodes in ELoops (for now) — haimes / detail
  52. Another def file fix — haimes / detail
  53. Add setTessParam to lite def file — haimes / detail
  54. Updates to egadslite to support pyEGADS.egadslite — galbramc / detail
  55. First cut at pyEGADS.egadslite — galbramc / detail
  56. Remove manifest files for jlEGADS. Check array size for makeGeometry. Fix makeGeometry tests. Some references fixes. — galbramc / detail
  57. Fix some pyEGADS refrencing. Check data length for pyEGADS makeGeometry. — galbramc / detail
  58. Fix topObj during EGADSlite import — haimes / detail
  59. Add U periodic options to EG_approximate for surfaces — haimes / detail
  60. try — docampo / detail
  61. Allow for full attribution in EGADSlite — haimes / detail
  62. Add reference counting for jlEGADS finalizers. Lots of other jlEGADS cleanup. — galbramc / detail
  63. Fix for last commit — galbramc / detail
  64. Various Julia cleanup and fixes — galbramc / detail
  65. Update julia Makefile — galbramc / detail
  66. Remove jlEGADS environment variable — galbramc / detail
  67. added package install — docampo / detail
  68. added Makefile — docampo / detail
  69. first jlEGADS — docampo / detail
  70. First cut at fitting periodic curves — haimes / detail
  71. Fix a number of bugs in the Effective Topology — haimes / detail
  72. Uprev EGADS — haimes / detail
  73. Forgot DARWIN64.clang and DARWIN64.gcc — galbramc / detail
  74. Move DARWIN64 ifort link flag to support EGADS test builds — galbramc / detail
  75. Fix for ifort on newer versions of macOS — galbramc / detail
  76. numpy suppression — galbramc / detail
  77. Always build egadslite — galbramc / detail
  78. Update julia path on linux — galbramc / detail
  79. Remove jlEGADS from Sanitizer and Valgrind — galbramc / detail
  80. mark helix3 and frustrum3 as expected successes — jfdannen / detail
  81. Add jlEGADS test to Sanitizer and valgrind — galbramc / detail
  82. Correct path — galbramc / detail
  83. Add julia testing — galbramc / detail
  84. Small fix in xfoil executables — galbramc / detail
  85. Higher output precision xfoil executables — galbramc / detail
  86. Upgrade to AFLR 10.22.23 — galbramc / detail
  87. Another regsx fix — galbramc / detail
  88. Fix email template — galbramc / detail
  89. Ignore ifort warnings — galbramc / detail

#165 (Jul 12, 2022, 9:54:17 PM)

  1. Fix sensitivites with AFLR and TRANSP. Fix check for return codes in caps_info/caps_analysisInfo in pyCAPS. — galbramc / detail
  2. Don't run corsairlite scipy exemples without scipy — galbramc / detail
  3. add wingBEM2 test case; add duct7 test case; add replicate.udc and ereped6* and ereped7 test cases; add sharpFin test case; typo fix in documentation — jfdannen / detail
  4. Use updated wget for ESP_PreBuilt.bat — galbramc / detail
  5. Beta script updates — galbramc / detail

#164 (Jul 12, 2022, 9:09:33 AM)

  1. Remove unused variable to fix stanalizer — galbramc / detail
  2. cCAPS example fix — galbramc / detail
  3. Stanalizer fixes — galbramc / detail
  4. Fixes for dealing with TRANSP surfaces. Add AFLR3 TRANSP tests. — galbramc / detail

#163 (Jul 10, 2022, 5:00:04 AM)

  1. Turn off wake for all tetgen examples — galbramc / detail
  2. protect for degenerate Edges in sensCSM; fix ESP bug that dropped user out of CapsMode if a non-existant pyscript is attempted; update ESP_QuickReference — jfdannen / detail

#162 (Jul 9, 2022, 1:36:12 PM)

  1. Fix SU2 to support internal BCs. TetGen does not support internal surfaces. — galbramc / detail
  2. Add TRANSP_UG3_GBC testing. Don't allow TRANSP_SRC_UG3_GBC for AFLR with CAPS. echo exectables when running tests. — galbramc / detail
  3. Codify corsairlite phasing example as commandline scripts — galbramc / detail
  4. Only set signals from the python Main Thread. — galbramc / detail
  5. Catch import signal ValueError. — galbramc / detail
  6. CorsairLite doc update — haimes / detail
  7. Various small cleanups — galbramc / detail
  8. Fix signal handling in pyCAPS — galbramc / detail
  9. Small Makefile cleanup — galbramc / detail
  10. Add intentPhase to corsair capsPhase examples — galbramc / detail
  11. Upgrade AFLR version check to 10.22.22 libraries. Add more AFLR3 tests. — galbramc / detail
  12. Update training solution — galbramc / detail
  13. Fix bug not removing .tessType for Tri EGADS tessellations — galbramc / detail
  14. Updated training box wing to remove capsIgnore — galbramc / detail
  15. Cleanup corsair examples a little — galbramc / detail
  16. Add intentPhrase and history to pyCAPS — galbramc / detail
  17. Replaced viscous wall bc to inviscid for Wing1 — nitin / detail
  18. changed viscous wall bc to inviscid for Wing1 — nitin / detail
  19. delete corsairlite example multifidelity — dongjoon / detail
  20. Bug fix for cart3d with geometric sensitvities — galbramc / detail
  21. Add symmetry inputs to cart3d — galbramc / detail
  22. Updated SU2 AIM for compatibility with the latest (SU2-7.3.1) — nitin / detail
  23. Protect value index from changes in aimCalcOutput — galbramc / detail
  24. lint fixes — galbramc / detail
  25. fix tabs to spasces — dongjoon / detail
  26. Fix warning — galbramc / detail
  27. Fix aimTransfer to use surface to volume maps — galbramc / detail
  28. modify corsair examples to generate pickle — dongjoon / detail
  29. Lint and memory leak fixes — galbramc / detail
  30. delete duplicate examples/corsairlite/phase_design directory — dongjoon / detail
  31. Update multifidelity examples with no hard coded paths and phaseStart — dongjoon / detail
  32. Update corsairlite examples with multifidelity phasing — dongjoon / detail
  33. Fix copy/pasta mistake — galbramc / detail
  34. Fix boundary indexing for aflr3 to mach capsGroup indexing — galbramc / detail
  35. lint fixes — galbramc / detail
  36. Fix AFLR NMakefiles — galbramc / detail
  37. Plug memory leak in aflr3 and run new Symmetry testing with MESH — galbramc / detail
  38. Reorder AFLR libraries — galbramc / detail
  39. Upgrade to AFLR 10.22.21 with support for sensitvities on symmetry planes — galbramc / detail
  40. fix bug in computing tessellation sensitivities for a NodeBody — jfdannen / detail
  41. try invEvaluateGuess for very hard cases in udfCompare; update compare1b test case — jfdannen / detail
  42. fix bug in IDE that did not properly journal updates to the intent phrase — jfdannen / detail
  43. second cut at tutorial6 write-up — jfdannen / detail
  44. more documentation updates — jfdannen / detail
  45. remove erroneous non-manifold warning associated with NodeBodys — jfdannen / detail
  46. fix bug when visualizing new-style quads and .mixed=NULL; do not process _makeQuads on a Face if the Body has a _makeQuads Attribute — jfdannen / detail
  47. fix bug in ESP.js that did not properly clear the display when velocities were cleared — jfdannen / detail
  48. fix bug that did not persist Edge Attributes to new non-manifold Edges that are generated by EXTRACTing Faces; add attributeA test case; add duct6 test case — jfdannen / detail
  49. fix bug in ESP.js associated with clearing velocities for a scalar; update tutorials 1 to 5 for v1.21 — jfdannen / detail
  50. modify way python is initialized and finalized to work around a problem in numpy - take 2 — jfdannen / detail
  51. modify way python is initialized and finalized to work around a problem in numpy — jfdannen / detail
  52. copy CFGPMTRs and DESPMTRs from temporary CAPS MODL back into the ESP MODL in CapsCommit and CapsSuspend; automatically rebuild MODL in CapsCommit and CapsSuspend — jfdannen / detail
  53. fix recently-added typo to timPyscript.c — jfdannen / detail
  54. add sanitize-thread target to Makefile.DARWIN64 — jfdannen / detail
  55. fix error message associated with not enough arguments to SPHERE; fix bug when user double-clicks in MessageWindow (Chrome and Safari) — jfdannen / detail
  56. Add _USE_MATH_DEFINES for MSVC — galbramc / detail
  57. Add error functions erf and erfc to Surreal — galbramc / detail
  58. README update — haimes / detail
  59. Config update — haimes / detail
  60. Fix bug returning wrong parameters in pyEGADS effectiveMap — galbramc / detail
  61. More defensive Surreal defines — galbramc / detail
  62. Some config updates — haimes / detail
  63. Trying to limit pip install conflits. Cleaner batch pip install function. — galbramc / detail
  64. Skip some long running tests with Sanitizer and Valgrind — galbramc / detail
  65. Missed one more — galbramc / detail
  66. More python suppressions — galbramc / detail
  67. Suppress AFLR uninitialized data — galbramc / detail
  68. Yet more numpy suppressions — galbramc / detail
  69. Another numpy suppression — galbramc / detail
  70. Upgrade to AFLR 10.22.22 on Windoze. Suppress numpy leak. — galbramc / detail
  71. Upgrade to AFLR 10.22.22 on Unix systems — galbramc / detail
  72. Upgrade jenkins to SU2 7.3.1 — galbramc / detail
  73. Look for more tessellation errors in the log file — galbramc / detail
  74. Upgrade to AFLR 10.22.21 — galbramc / detail

#161 (Jun 13, 2022, 1:00:02 AM)

  1. os independent file paths for corsair examples — galbramc / detail
  2. Ignore corsairlite for legacy regression tests — galbramc / detail
  3. Use cdef in fun3dNamelist.pyx to avoid segfault during garbage collection — galbramc / detail
  4. Forgot a file — galbramc / detail
  5. embedded python updates — galbramc / detail
  6. Fix divide by zero — galbramc / detail
  7. fix setting of PYTHONHOME for a double-quoted string — haimes / detail
  8. generate WARNING if importing an .egads file that has an Attribute that cannot be evaluated — jfdannen / detail
  9. Fix up EMP_ThreadDestroy for the M1 — haimes / detail
  10. Stanalizer fixes for Lazy Surreal — galbramc / detail
  11. Use std c++11 for Surreal tests. Fix in depricated EG_loadTess. — galbramc / detail
  12. Use Lazy Expressions for Surreal for slightly better performance on windoze — galbramc / detail
  13. Minor config changes — haimes / detail
  14. Add random sleep for Windoze pip install to avoid install conflicts — galbramc / detail

#159 (Jun 11, 2022, 12:14:09 AM)

  1. Relax tolerances for cvxopt — galbramc / detail
  2. Replaced serveCSM with serveESP — nitin / detail
  3. added figures to html version, consistent with pdf — nitin / detail
  4. Update corsairlite kulfan models — dongjoon / detail
  5. Copy hoburg.pl for corsair testing. Ignore corsair in C regressions — galbramc / detail
  6. Cleanup corsairlite capsPhase examples a little — galbramc / detail
  7. Run corsairlite qp for MINIMAL testing — galbramc / detail
  8. Fix examples copy for corsairlite — galbramc / detail
  9. Enable corsairlite test again. — galbramc / detail
  10. Don't test corsairlite yet — galbramc / detail
  11. Defer scipy import error for corsair kulfan datatype — galbramc / detail
  12. Add corsairlite examples — galbramc / detail
  13. Add corsairlite — galbramc / detail
  14. Some dox cleanup — galbramc / detail
  15. Add caps.phaseNewCSM to pyCAPS. Enalbe use of phaseNewCSM in pyCAPS Problem. — galbramc / detail
  16. Updated CAPS overview dox, *_dox.cfg, *_header.tex to fix warnings — nitin / detail
  17. Add caps_phaseNewCSM to the CAPS API — haimes / detail
  18. Fix bug specifuing Mesh_Sizing on Edges for aflr4 — galbramc / detail
  19. Update phasing tests. Fix Multiple_Mesh for TetGen. — galbramc / detail
  20. Make updates fix directory removal under error — haimes / detail
  21. Cleanup phasing in pyCAPS — galbramc / detail
  22. Include IDE's parent info in phaseUtil and fix some minor issues — haimes / detail
  23. Complete phaseUtil — haimes / detail
  24. Add AnalysisIn derivatives in sens file for tacs and fun3d — galbramc / detail
  25. Add the phaseUtil app -- still needs filling out — haimes / detail
  26. Change skeleton AIM so it can be used as part of the ESP IDE tutorial — haimes / detail
  27. Add mses executables to ESPbfr — galbramc / detail
  28. Protect FunID read against a corrupt file — haimes / detail
  29. Revert adding debug in Problem.__init__ — galbramc / detail
  30. Output file position on Journal read during debugging — haimes / detail
  31. Guard pyCAPS.Problem.__init__ caps calls with caps_debug — galbramc / detail
  32. Stanalizer fix — galbramc / detail
  33. Fix memory leak in Cart3D data transfer — galbramc / detail
  34. Fix bug in last commit. Cleanup caps_preAnalysiX. Add forgotten Cart3D FSI example. — galbramc / detail
  35. Put user name and date/time stamp in capsLock — haimes / detail
  36. fix memory leak in msesAIM — galbramc / detail
  37. Lint fixes — galbramc / detail
  38. another scan-build fix — galbramc / detail
  39. scan-build fixes — galbramc / detail
  40. Add environment variable to output Journal entries as they process — haimes / detail
  41. Fix warning — galbramc / detail
  42. Dedicated function to call aim_UpdateState. Re-formulate phasing/continuation in pyCAPS Problem constructor. Indent AIM error messages. Add Displacement data transfer to Cart3D for FSI. Only allow cals to caps_getDataSet in aimPreAnalysis. Do not allow disk manipulation in aimUpdateState. — galbramc / detail
  43. Fix childByName when in debug mode — haimes / detail
  44. Finally the correct Windoze executable — galbramc / detail
  45. Update avl exec on Linux and Windoze — galbramc / detail
  46. Fix both the documentation and the code so they are consistent with adding const to the AIM functions Pre and Exec — haimes / detail
  47. update jfd3AIM to adjust for recent API changes — jfdannen / detail
  48. Reinstate no intent phrase flag — haimes / detail
  49. Try that again — galbramc / detail
  50. Update macOS avl executables without X11 dylib dependency — galbramc / detail
  51. Don't journal on initPhrase calls in caps_open — haimes / detail
  52. Revert last commit — haimes / detail
  53. Turn off journaling of intentPhrase in caps_open — haimes / detail
  54. lint fixes — galbramc / detail
  55. So many def files... — galbramc / detail
  56. Fixes for delaundo — galbramc / detail
  57. Another forgotten def file — galbramc / detail
  58. Missed one — galbramc / detail
  59. aimUpdateState is no longer optional — galbramc / detail
  60. Forgot def file — galbramc / detail
  61. Add aimUpdateState to su2 AIM, and fix warnings — galbramc / detail
  62. Add aimUpdateState to nastran and poinwise AIMs — galbramc / detail
  63. Add aimUpdateState to hsm, interference, masstran, mystran, and aimUtil.h — galbramc / detail
  64. Fix writing a Integer Value Object that is NULL — haimes / detail
  65. Add aimUpdateState to fun3dAIM. Call aim_UpdateState in post analysis during reload. — galbramc / detail
  66. remove debug print statement — jfdannen / detail
  67. remove side effects from assert() statements, which work differently on Windoze — jfdannen / detail
  68. allow pyscripts to end without a newline — jfdannen / detail
  69. fix bug associated with finishing non-manifold WireBodys — jfdannen / detail
  70. allow overlay to execute a mesg even when in TIM_EXECUTING; update tutorial6 files — jfdannen / detail
  71. remove debug printing of long cast — jfdannen / detail
  72. suppress profile printing if outLevel==0; fix filename updates in ESP after saving new file and after exporting the feature tree — jfdannen / detail
  73. numerous updates — jfdannen / detail
  74. remove try/catch from IDE demo pyscript files; use skeletonAIM instead of jfd3AIM in IDE demo; only show Display table entries (Node/Edge/Face/Csystem/Pips) if it contains any child entities; fix sensitivity bug associated with UDP/UDFs that requires finite differences coupled with an OUTPMTR that is tied to a mass property; add udpTest; initial commit of new timGloves and associated files — jfdannen / detail
  75. add ocsmPrintProfile to explicitly print the profile of CPU usage (instead of printing automatically during ocsmBuild); automagically run pyscript.py when running in continuation mode; change Cvals via script when in CAPS mode; change DESPMTRs via script when in CAPS mode; offer to steal capsLock in CAPS mode; add ogive* test cases — jfdannen / detail
  76. protect non-script caps_ calls with caps_debug; fix bug in timCapsMode that sometimes reverted to original .csm file; update demo_tim.txt; remove duplicate entries from ListHistory; fix typo in ESP-help associated with EVALUATE command — jfdannen / detail
  77. remove -caps command-line option (replaced by timCapsMode); add cleanup process in pyscript.py only between scripts (not at end); fix bug that did not profile MESSAGE statements; reset profile info at beginning of ocsmBuild; fix bug that did not print OUTPMTRs during call to ocsmPrintPmtrs — jfdannen / detail
  78. compile serveESP in Makefile and NMakefile — jfdannen / detail
  79. inform user when CAPS overrides .csm value(s); first attemp at re-running pyscript in continuation mode — jfdannen / detail
  80. add indication for Phases that are a predecessor of current Phase in ListPhases — jfdannen / detail
  81. add calls to caps_debug around List* messages in timCapsMode; update data/ide/*.py to directly rebuild the geometry; reorder Caps menu options — jfdannen / detail
  82. fix ListHistory in timCapsMode — jfdannen / detail
  83. modify demo_caps.txt and demo_tim.txt; make changes needed to make both demos operate properly — jfdannen / detail
  84. fix stanalyzer error — jfdannen / detail
  85. more modifications to support timCapsMode — jfdannen / detail
  86. temporary commit to be used for debugging — jfdannen / detail
  87. remove debug print statement — jfdannen / detail
  88. additional bugs fixed in timCapsMode; add unequalRules test case to show how to RULE between xsects with different nnumber of Edges — jfdannen / detail
  89. remove temporary printf statements — jfdannen / detail
  90. rename capsIde to capsMode; restrict SPHEREs to have a positive radius; allow multiple ESP-help tabs in browser; partial implementation of timCapsMode — jfdannen / detail
  91. Config change for Windows to always start browsers in the background — haimes / detail
  92. Windows configuration change to promote Edge — haimes / detail
  93. update README.txt with note about Safari block-ups — jfdannen / detail
  94. Add makeNmWireBody to pyEGADS — galbramc / detail
  95. Don't clobber LD_LIBRARY_PATH with ESPenv — galbramc / detail
  96. Configuration changes and signature change for EG_makeNmWireBody — haimes / detail
  97. Documentation updates — haimes / detail
  98. More Doc updates — haimes / detail
  99. Some Documentation updates — haimes / detail
  100. Fix the fix for memory issue with getting units from imports — haimes / detail
  101. Fix memory issue with getting units from imports — haimes / detail
  102. Document some of the new updates to SETP/IGES handling — haimes / detail
  103. Support units on both STEP and IGES exports — haimes / detail
  104. Add an option to do less checking with makeTopology on a Loop — haimes / detail
  105. Another touch of splitPeriodics — haimes / detail
  106. Try that again -- revert and get the if statement correct! — haimes / detail
  107. Try that again — haimes / detail
  108. Don't update the shape in splitPeriodics if it fails curve simplification — haimes / detail
  109. Install corsairlite dependencies for ESP_Beta testing — galbramc / detail
  110. Python suppression — galbramc / detail
  111. More error checking in ESPbfr.bat — galbramc / detail
  112. Attempt Python 3.9.13 again — galbramc / detail
  113. macM1 uses 3.8.11 — galbramc / detail
  114. Revert back to ESP_Python 3.8.9 for now — galbramc / detail
  115. install matplotlib and scipy for corsairlite — galbramc / detail
  116. install packaging for corsairlite — galbramc / detail
  117. Upgade to ESP_Python-3.9.13 — galbramc / detail
  118. Install dill and cvxopt — galbramc / detail
  119. Fix funny permissions on DARWIN64 mastros — galbramc / detail
  120. Update macOS intel/M1 mastros with old macOS target — galbramc / detail
  121. astros has array out-of-bound problems — galbramc / detail
  122. Don't define ASTROS_ROOT on machines without astros — galbramc / detail
  123. Restore mastros on macOS, macM1, and restore astros on Linux — galbramc / detail
  124. Uses msys64 wget on Windoze — galbramc / detail
  125. fix formatting error associated with expected errors for frustrum3 and helix3 — jfdannen / detail
  126. fix formatting error associated with expected errors for frustrum3 and helix3 — jfdannen / detail
  127. add frustrum3 and helix3 to expected failures — jfdannen / detail
  128. Windoze requires exe extension — galbramc / detail
  129. Add mses to ESPbfr.tgz with tighter convergence tolerance — galbramc / detail
  130. Another python suppression — galbramc / detail
  131. change testing from serveCSM to serveESP — jfdannen / detail
  132. Third time is the charm — galbramc / detail
  133. Try updating Windoze avl again... — galbramc / detail
  134. Update avl exec on Linux and Windoze — galbramc / detail
  135. Try that again — galbramc / detail
  136. Update macOS avl executables without X11 dylib dependency — galbramc / detail

#158 (Apr 8, 2022, 11:47:49 AM)

  1. Fix lint issue on skeletonAIM — haimes / detail
  2. Fix warnings from last commit — haimes / detail
  3. Protect against writing NULL double value objects; Update skeleton AIM — haimes / detail
  4. Fix applying offset to NULL pointer — galbramc / detail

#157 (Apr 7, 2022, 3:28:24 PM)

  1. Update avl reference data. Plug memory leak in egadsTess — galbramc / detail
  2. Fix aimUpdateState in egadsTess AIM — galbramc / detail
  3. Don't run mystran_PyTest on Windoze — galbramc / detail
  4. Forgot def file — galbramc / detail
  5. Remove on too many lines — galbramc / detail
  6. Fix warning from last commit — galbramc / detail
  7. More aimUpdateState updates. Remove some redundant tests — galbramc / detail
  8. Add aimUpdateState to a number of AIMs. Turn on CLink between Phases. — galbramc / detail
  9. Misseed one astros in training — galbramc / detail
  10. Allow aim_relPath to cross Phases — haimes / detail
  11. Guard againts running mastros test when mastros does not exist — galbramc / detail
  12. Fixes for calling aimUpdateState in caps_preAnalysis — galbramc / detail
  13. Plug memory leak — galbramc / detail
  14. More mses tests. Add aimUpdateStatus to avlAIM. Other small fixes. — galbramc / detail
  15. Allow for debug mode coming into caps_outputObjects — haimes / detail
  16. Allow for debug mode coming into caps_outputObjects — haimes / detail
  17. in pyCAPS.Problem(), use running Caps Problem if it already exists — jfdannen / detail
  18. Add phasing xfoil tests — galbramc / detail
  19. Add Phasing test to pyCAPS unit testing. Fix probemName misspelling. — galbramc / detail
  20. Turn off Analysis links for now — haimes / detail
  21. Fix the warnings from the last commit — haimes / detail
  22. First cut at Analysis disk links for new Phases — haimes / detail
  23. Reduce execution time in legacy xfoil example — galbramc / detail
  24. Faster geometry construction for ESP_Commit — galbramc / detail
  25. Add some more timing — galbramc / detail
  26. Remove incorrect CL derivatives from msesAIM. Run one test per aim for MINIMAL tests. — galbramc / detail
  27. Fix last commit — haimes / detail
  28. Prepare the AIM utilities for Analysis directory links — haimes / detail
  29. Missed one free — galbramc / detail
  30. lint fixes — galbramc / detail
  31. Use aimUpdateState in msesAIM. Cache sensitvities in msesAIM. Expose more inputs in msesAIM. — galbramc / detail
  32. Implement undocumented caps_modifiedDesPmtrs so that the IDE can know which DESPMTRs have been modified by CAPS — haimes / detail
  33. Some more mses updates and regression script fix — galbramc / detail
  34. Don't check alpha for TargetCL Cart3D test — galbramc / detail
  35. Add TargetCL to Cart3D and enable Cart3D 2D mode — galbramc / detail
  36. Try that again — galbramc / detail
  37. Fix infinite loop in convert_doubleToString — galbramc / detail
  38. Fix calling updateState when an AnalysisIn Object has changed — haimes / detail
  39. Add the AIM function aimUpdateState which will allays be the first function executed in the sequence — haimes / detail
  40. A few more mses corrections — galbramc / detail
  41. Adjust mses data — galbramc / detail
  42. mses updates — galbramc / detail
  43. Remove caps_resetAnalysis from the API — haimes / detail
  44. Use curvature weighted point distribution for avl/xfoil/mses — galbramc / detail
  45. Remove caps_resetAnalysis from pyCAPS — galbramc / detail
  46. Don't write 'empty' rubberize files — galbramc / detail
  47. "fix" to fun3dAIM when to write/read shape sensitivities — ryan / detail
  48. Fix intent phrase sync issue on initialization — haimes / detail
  49. Add initial history for GeometryIn, Analysis and AnalysisIn Objects — haimes / detail
  50. Fix journaling so that changing DESPMTR (GeometryIn) dimensioning is handled — haimes / detail
  51. Remove old caps_open CSM reload option and replace it with the new — haimes / detail
  52. Fix history addition for reloading CSM and new GeometryIn Value Object — haimes / detail
  53. Relax tolerance and try to get consistent xfoil results on all machines — galbramc / detail
  54. Update mses data. — galbramc / detail
  55. Update xfoil data. Don't set OCS_CONPMTR when journaling — galbramc / detail
  56. More robust xfoil execution — galbramc / detail
  57. Update legacy test xfoil data — galbramc / detail
  58. Adjust tolerances — galbramc / detail
  59. Add mses OpenMDAO example — galbramc / detail
  60. Missed some test that needed updated reference data — galbramc / detail
  61. Use EGADS tessellator for airfoil files — galbramc / detail
  62. Fix warning clang 13.1.6 warning — galbramc / detail
  63. Fix full path on CSM filename for oPNewCSM option with caps_open — haimes / detail
  64. Forgot to close a file — galbramc / detail
  65. Fix msesAIM sensitivities — galbramc / detail
  66. Support a directory of CSM/UDC files at init or CSM reload for caps_open; test for all analyses past post for completing a Phase — haimes / detail
  67. Don't write fun3d Rubberize files when not needed — galbramc / detail
  68. Don't write out body information to rubber.data for only AnalysisIn sensitvities — galbramc / detail
  69. Stanalizer fix — galbramc / detail
  70. Add avl 3.40 min version error message — galbramc / detail
  71. Remove libquadmath dependency on DARWIN64 — galbramc / detail
  72. Add avl 3.40 to ESPbfr — galbramc / detail
  73. Plug one more memory leak — galbramc / detail
  74. Fix array bound issue in avlmrf — galbramc / detail
  75. Move variable declarations to top of functions — galbramc / detail
  76. Fix avl NMakefile. Fix uninitialized data. — galbramc / detail
  77. Lint and memory leak fixes — galbramc / detail
  78. aimBackdoor has been removed from the avl AIM — galbramc / detail
  79. avl fixes — galbramc / detail
  80. Upgrade to avl 3.41 — galbramc / detail
  81. Change caps_delete to caps_markForDelete in pyCAPS — galbramc / detail
  82. Change caps_delete to caps_markForDelete — haimes / detail
  83. Preserve pyCAPS objects when calling caps_delete. caps_delete and pyCAPS delete will be renamed. — galbramc / detail
  84. revert previous commit — jfdannen / detail
  85. fix previous commit — jfdannen / detail
  86. remove unneeded windoze headers from tim.c — jfdannen / detail
  87. Fix continuation marker in serveESP.mak — haimes / detail
  88. Fix headers for Windows in timPyscript — haimes / detail
  89. fix windoze problem with includes in pervious commit — jfdannen / detail
  90. only allow one open ESP-help tab at a time; implement timFlowchart; modify interpolateUpper.py and conserveUpper.py to demonstrate flowcharting tool — jfdannen / detail
  91. allow for .py scripts to run both in serveESP -caps and in python; add ocsmSetAuxPtr and ocsmGetAuxPtr in OpenCSM API — jfdannen / detail
  92. fix stanalyzer problem in timCapsIde.c — jfdannen / detail
  93. update serveESP.mak in previous commit to include timCapsIde.obj — jfdannen / detail
  94. change display scaling so that if there is only one NodeBody it still gets displayed; initial implementation of timCapsIde — jfdannen / detail
  95. allow nesting of TIMs to MAX_TIM_NESTING=10 — jfdannen / detail
  96. fix bugs in File->Open and File->Edit for cases with and without UDCs — jfdannen / detail
  97. fix bug in ocsmSave that did not write DIMENSION statements correctly when writing either a .csm or .cpc file with a variable nrow and/or ncol; update calls to caps_open to handle UDCs — jfdannen / detail
  98. change -caps mode to accommodate both fore- and back-slashed in filenames; use ocsmTessellate in serveESP; update serveESP.make and serveESP.mak to define USING_CAPS (so that serveESP can be built without CAPS); remove unneeded code that generates a warning with new clang version; use updated caps_open so that cases with UDCs work properly — jfdannen / detail
  99. remove unneeded code that generates a clang warning — jfdannen / detail
  100. add ocsmAdjustUDCs to OpenCSM API to copy all UDCs to the same directory as the .csm file; implement File->New, File->Open, and File-Edit->OK in caps mode — jfdannen / detail
  101. Fix scan-build warning from last commit — haimes / detail
  102. Implement EG_makeNMWireBody — haimes / detail
  103. Straighten out non-manifold WireBodies — haimes / detail
  104. Supress informational message on initial triangulation problem with Planar Faces — haimes / detail
  105. Guard against infinite loop on planar surfaces with bad initial triangulation — haimes / detail
  106. More support for tessellations without face tessellations — galbramc / detail
  107. Enable EG_finishTess to fill in missing face tessellation. — galbramc / detail
  108. Remove warnings generated by a new compiler — haimes / detail
  109. Remove redundant testing — galbramc / detail
  110. Remove mastros.exe on macys and viggen because it does not work — galbramc / detail
  111. Don't look for python deprecation warnings — galbramc / detail
  112. Look for python deprecation warnings — galbramc / detail
  113. Limit threading with OMP_NUM_THREADS — galbramc / detail
  114. Don't pipe sanitize pyOCSM testing — galbramc / detail
  115. add numpy valgrind suppression — galbramc / detail
  116. Look for python unit test timing in log parse — galbramc / detail
  117. Reduce CAPS ESP_Commit testing — galbramc / detail
  118. scipy is now available on Monterey — galbramc / detail
  119. Look for planar meshing errors — galbramc / detail
  120. Update avl and DAWIN64/mastros.exe executables — galbramc / detail
  121. astros.exe does not work on Monterey — galbramc / detail
  122. Upgrade macOS santizer library path — galbramc / detail
  123. Remove libquadmath dependency on DARWIN64 — galbramc / detail
  124. Correct avl version number to 3.40 — galbramc / detail
  125. Reduce avl NVMAX for DARWINM1 — galbramc / detail
  126. Upgrade to avl 3.41 — galbramc / detail

#155 (Mar 12, 2022, 6:41:34 AM)

  1. Revert astros example — galbramc / detail
  2. Fix last scan-build error — haimes / detail
  3. Add history to geometry objects that are either new inputs or all outputs when reloading a CSM file — haimes / detail
  4. Fix scan-build complaint — haimes / detail
  5. lint fixes — galbramc / detail
  6. Restote attributes and history on geometry objects with a CSM reload — haimes / detail
  7. Update legacy su2 test — galbramc / detail
  8. Add some egadsTess NodeBody tests — galbramc / detail
  9. Fix BC indexing for meshRef — galbramc / detail
  10. Add history to geometry outs and other minor updates — haimes / detail
  11. Use enum instead of numerical value for checking stFlag — haimes / detail
  12. Also return row vector in caps_getValue — galbramc / detail
  13. Patch up some potential problems with restart/journaling and changing OpenCSM dimensions — haimes / detail
  14. Add read-only option on caps_open -- not well tested — haimes / detail
  15. remove jfd3AIM from aim/Makefile — jfdannen / detail
  16. add hfd3AIM for testing Bounds and DataSets — jfdannen / detail
  17. Add new caps_delete to pyCAPS — galbramc / detail
  18. Reconcile addHistory with writeSerialNum — haimes / detail
  19. Add history to caps_makeValue — haimes / detail
  20. Upgrade to mystran 13.3a — galbramc / detail
  21. Take away double CSM build ifdef in caps_open — haimes / detail
  22. Add CAPS open option to ignore object delete markers; return delete status in caps_info — haimes / detail
  23. Actually write the serial number in caps_writeSerialNum — haimes / detail
  24. Update legacy example — galbramc / detail
  25. lint fix and update astros Design example — galbramc / detail
  26. Stanalizer fix — galbramc / detail
  27. Fix warning — galbramc / detail
  28. Fix memory leaks. Journal caps_setValue(NULL). — galbramc / detail
  29. Only write sNum to the Problem file if that is all that is needed — haimes / detail
  30. Minor fixes & move sNum to beginning of problem file — haimes / detail
  31. Remove Bound directory during deletion and fix up internal Parameter index when there are deletions — haimes / detail
  32. re-expose su2 variables for aeroelastic cases — nitin / detail
  33. New phase deletion fleshed out — haimes / detail
  34. Remove capsLock beforey renmaing a phase on caps_close — galbramc / detail
  35. Fix logic for tracking pyCPAS ProblemState — galbramc / detail
  36. Update makefiles to run skeletonAIM example — galbramc / detail
  37. Fix skeleton AIM — galbramc / detail
  38. Hold on to and delete User Value Objects at a complete CAPS close; change the meaning of caps_delete (still need to do the deletion at start of a new phase) — haimes / detail
  39. SU2 AIM updated for SU2-7.3.0, Changed boundary markers from numbers to strings — nitin / detail
  40. SU2 AIM updated for SU2-7.3.0, Changed boundary markers from numbers to strings — nitin / detail
  41. Fix mses memory leak — galbramc / detail
  42. Fix legacy script. Add caps_intentPhrase to pyCAPS. — galbramc / detail
  43. scan-build fix — galbramc / detail
  44. Add Cheby_Mode output and sensitvities to msesAIM — galbramc / detail
  45. lint fixes — galbramc / detail
  46. Fix warnings — galbramc / detail
  47. Make relaoding AIMs lazy. Remove unused problem->exec. Enable Journaling/Phasing for aflr3 and fun3d. — galbramc / detail
  48. Fix caps_size in debug mode for Dynamic Outputs — haimes / detail
  49. Implementation of new History tracking — haimes / detail
  50. dox fix — ryan / detail
  51. lint fixes for aflr4 — galbramc / detail
  52. Enable phasing and journaling for aflr4AIM — galbramc / detail
  53. lint fixes — galbramc / detail
  54. Don't call analsysisInfo with problem object — galbramc / detail
  55. scan-build fixes — galbramc / detail
  56. Close ugrid file after reading header — galbramc / detail
  57. Use correct CAPSLONG in pyCAPS on Windoze — galbramc / detail
  58. Update caps_jrnlRead with funID — galbramc / detail
  59. Add funID to caps_jrnlWrite. Fix call to caps_ownerInfo. Fix memory leak in egadsTess. — galbramc / detail
  60. Update pyCAPS for new history tracking. Fixes in su2. Test Phasing with egads/tetgen/su2. — galbramc / detail
  61. Move CAPS towards the new history tracking — haimes / detail
  62. Revert changes in CAPS/src from previous accidental commit — galbramc / detail
  63. Fix training files from last commit. Remove bodyTessMappingStruct in lieu of .mixed attribute. Fixes for Journaling. Fixes for autoExec when the upstream AIM unputs are dirty. Store mesh on disk for egadsTessAIM to work with Journaling and Phasing. Update tetgenAIM to work with Journaling and Phasing. Add tetgen pyCAPS unit tests. — galbramc / detail
  64. Remove caps_freeOwner from pyCAPS. Add commented out caps_getHistory and caps_intentPhrase for the future. — galbramc / detail
  65. Put caps_freeOwner back in the def file — haimes / detail
  66. remove caps_setOwner from caps.py — jfdannen / detail
  67. Prepare for history changes — haimes / detail
  68. Add aim_readBinaryUgridHeader — galbramc / detail
  69. Remove addHistory and getHistory from pyCAPS — galbramc / detail
  70. Allow for removing the entire phase at caps_close — haimes / detail
  71. Don't call caps_delete in pyCAPS for now — galbramc / detail
  72. Add aim_rmFile — galbramc / detail
  73. Change pyCAPS caps_triangulate to caps_getTriangles — galbramc / detail
  74. Add caps_printErrors; change caps_triangulate to caps_getTriangles and other fixes — haimes / detail
  75. Some dox cleanup for tacsAIM — galbramc / detail
  76. tacs dox fix — ryan / detail
  77. fixes for 'Flow' directory switching for fun3dAIM — ryan / detail
  78. Update flowchart viewing and check for parameter value object uniqueness — haimes / detail
  79. Track closure of problem object to avoid calling caps_delete. This will leak attributes for now if the problem object is explicitly closed. — galbramc / detail
  80. Track closure of problem object to avoid calling caps_delete. This will leak attributes for now if the problem object is explicitly closed. — galbramc / detail
  81. Training fix — galbramc / detail
  82. Undo accidental change — galbramc / detail
  83. Some more training fixes — galbramc / detail
  84. Fix comments in training examples — galbramc / detail
  85. dox fixes — galbramc / detail
  86. Require using Design_Variable_Relation for bdf files. — galbramc / detail
  87. dox fix to clarify symmetry boundary condition requirements for fun3d — ryan / detail
  88. Disable autoExecution for functions in capsBase when in debug mode — haimes / detail
  89. Add newline to printing of tuples that are NULL — haimes / detail
  90. Fix memory leak from last commits — haimes / detail
  91. Remove debug prints — haimes / detail
  92. Update outputObjects and add Body and Tess count to AIM info for printObjects — haimes / detail
  93. more improvements for ListHistory; use caps_open in read-only mode in ListPhases and ListHistory — jfdannen / detail
  94. do not let esp.UpdateESP rebuild (force rebuild in .py script) — jfdannen / detail
  95. extend ListHistory to GEOMETRYIN and GEOMETRYOUT — jfdannen / detail
  96. update Makefiles for cases without python; add many features to -caps mode; add data/ide/* to test -caps mode; fix many bugs in -caps mode; expose ocsmTessellate — jfdannen / detail
  97. Makefile updates to supress Python builds — haimes / detail
  98. first implementation of CAPS menu in ESP — jfdannen / detail
  99. allow branching of Phases when running serveESP in -caps mode — jfdannen / detail
  100. caps_PrintObjects signature changed — haimes / detail
  101. allow wild-carding (with 0) in SELECT ADD and SELECT SUB; update connect9a and connect9b; update Makefile.DARWIN64 for 3-way testing — jfdannen / detail
  102. in timViewer, use smoothing shading when looking at BOUNDS; allow wild-carding (with 0) in SELECT ADD and SELECT SUB; update Makefile.DARWIN64 for 3-way testing — jfdannen / detail
  103. add viewing of BOUNDS in timViewer; allow GetToken in tim.c to work if the text string does not end in a delimiter; improve error message generated when SKVAR does not contain triplets; improve CAPS summary at end if running in -caps mode — jfdannen / detail
  104. in timViewer, use a tessellation associated with the AIM if available, otherwise share the tessellation with OpenCSM — jfdannen / detail
  105. fix (long) cast in previous commit — jfdannen / detail
  106. extend STORE such that if index=-1, the first available index is used; fix bug that generated two copies of Bodys on stack when RESTOREing with index<0; if running serveESP in -caps mode, get MODL from CAPS (instead of building it in OpenCSM); extend timViewer to view Bodys associated with an AIM — jfdannen / detail
  107. Config change — haimes / detail
  108. Explicitly open files as bindary in EG_saveTess/EG_loadTess. — galbramc / detail
  109. Upgrade to MYSTRAN 13.3A on Windoze — galbramc / detail
  110. Upgrade to Mystran 13.3a — galbramc / detail
  111. Another python suppression — galbramc / detail
  112. Update jenkins to run skeletonAIM example — galbramc / detail
  113. Upgrade to SU2-v7.3.0 for Jenkins — galbramc / detail

#154 (Feb 12, 2022, 12:08:11 PM)

  1. Remove the regeneration of FieldIn DataSets during AIM restart — haimes / detail
  2. Revamp the way postAnalysis gets invoked during recycling — haimes / detail
  3. Fix jouranling bug — galbramc / detail
  4. Remove first OCSM build when starting a new phase updating the CSM script; don't write out value objects with the link nullified — haimes / detail
  5. Add None arguments to outputObjects — galbramc / detail
  6. Add stream output option to caps_outputObjects — haimes / detail
  7. Add outputObjects to capsProblem — galbramc / detail
  8. Fix bug found by Stanalyzer output — haimes / detail
  9. First cut at DongJoon's CAPS flowchart (debugging) — haimes / detail
  10. Fix caps_rmWild to support spaces in path — haimes / detail
  11. Try to get by Cart3D tests again — haimes / detail
  12. Make sure there is enough room for commands with quoted paths — haimes / detail
  13. Allow for spaces in Problem paths — haimes / detail
  14. Documentation updates and raise an error when caps_open is given a path with an included space — haimes / detail
  15. Fix last commit — galbramc / detail
  16. Try to plug unitialized values from last commit — haimes / detail
  17. Implement a debug mode that does not journal or autoExec/autoFill — haimes / detail
  18. Update .sens file for tacsAIM and fun3dAIM. Add Phasing test for avl. — galbramc / detail
  19. add functional to fun3d options (undocumented) — ryan / detail
  20. Fix writeTecplot — galbramc / detail
  21. Fix for updated pyCAPS constructors — galbramc / detail
  22. Fix last commit — galbramc / detail
  23. Another try at fixing the memory leak. Update pyCAPS for the new caps_triangulate. Add caps_phaseStatus to pyCAPS. — galbramc / detail
  24. Lets try that again — galbramc / detail
  25. Another memory leak fix — galbramc / detail
  26. Update caps.def — galbramc / detail
  27. Update aimTransferUtilTest to test caps_ functions — galbramc / detail
  28. Fix memory leaks from last commit — haimes / detail
  29. Add mesh information to the data returned from caps_triangulate — haimes / detail
  30. Adjust the calling sequence of caps_open to be consistent with the documentation — haimes / detail
  31. Fix deriv bug in caps_writeValue — galbramc / detail
  32. Fix memory corruption from cut & paste — haimes / detail
  33. Add file handling functions to the def file — haimes / detail
  34. Fix SegFault in caps_open CSM reload — haimes / detail
  35. Regression test fixes — galbramc / detail
  36. Update legacy test for PointwiseV18.5R1 — galbramc / detail
  37. Try to fix Stanalyzer error from last commit — haimes / detail
  38. Some corrections to Phase start with CSM reload — haimes / detail
  39. Try again to fix Stanalyzer (scan-build) errors — haimes / detail
  40. Try to fix Stanalyzer (scan-build) errors — haimes / detail
  41. Implementation of new phase startup with CSM script reload -- needs more tessting — haimes / detail
  42. Suppress lint — galbramc / detail
  43. Ignore zero lenth units string in caps_makeValue — galbramc / detail
  44. Fix warning and relax tolerance — galbramc / detail
  45. Fix for caps_getDeriv. Lots of updates to msesAIM — galbramc / detail
  46. Fix for caps_hasDeriv — galbramc / detail
  47. Allow pyCAPS.Problem construction with a capsObj — galbramc / detail
  48. Add missing reference quanteties — galbramc / detail
  49. Updated legacy 'truth' avl values — galbramc / detail
  50. Upgrade pointwimse AIM to support Pointwise 18.5. Rotate airfoils for astros. Add error checking for capsReference quantities in avl. — galbramc / detail
  51. fix bug that either hangs or seg faults serveESP if user closes browser while timPyscript is executing.  this is done by not calling Py_Finalize, which will result in a memory leak; fix bug that occasionally does not show pyscript file when serveESP is starterd with a .py file — jfdannen / detail
  52. improve error message when project.csm exists before the project directorry when running serveESP in -caps mode — jfdannen / detail
  53. improve the creation, editting, and removal of Caps Values in -caps mode; allow multi-valued Caps Values in ESP; add udfShadow and shadow* test cases — jfdannen / detail
  54. added numpts argument to udpKulfan, which helps avoid small dimples near the leading edge of a very thin section — jfdannen / detail
  55. add debug GRAFICs to udpKulfan — jfdannen / detail
  56. add SELECT LOOP to get Edges associated with a Loop in a Face; allow unary + or - in front of (multi-valued) Parameter when used alone in an expression; add sectionsWithHoles* test cases — jfdannen / detail
  57. create hollowSheet and hollowSolid test cases to show all possibile uses of HOLLOW command; improve description of HOLLOW in OpenCSM.h; if HOLLOW fails, input Body is left on the stack; add hollowE* test cases to show some of the limitations of HOLLOW; change timPython to timPyscript to get around problems with linking with python on Windoze; have OpenCSM return SUCCESS even when all the Bodys have been recycled; put note in ESP-help and OpenCSM.h that warns users that results from HOLLOW and SWEEP may be unpredicable due to OpenCASCADE issues; add timPlotter and plotter* test cases — jfdannen / detail
  58. fix bug in ocsmPrintEgo when called on BSPLINE Curve or Surface; require positive scale factors in udfNuscale when applied to a SolidBody; add initial version of timPlotter and test cases; change GetToken in tim.c so that tokens can have unlimited length — jfdannen / detail
  59. fix errors that occur because Py_Finalize cannot be used more than once (because of a numpy bug); fix bug that did not show Cvals for -caps runs in serveESP and Brchs when running without -caps — jfdannen / detail
  60. update Caps Value handling in ESP — jfdannen / detail
  61. do not show Branches in serveESP when in -caps mode; allow serveESP to change .csm file when in -caps mode; fix OpenCSM.c to remove stanalizer error — jfdannen / detail
  62. for a UDP/UDF, only use finite differences if one of its arguements is ATTRREAL and has non-zero dots — jfdannen / detail
  63. modify timPython to get around numPy problems associated with Py_Finalize; modify tim interface to allow a tim to perform cleanup at the end of serveESP — jfdannen / detail
  64. fix tessellation sensitivity bug associated with SheetBodys that are INTERSECTed with a SolidBody; add designP5* test cases; fix Node sensitivity for two parallel Edges; move timViewer from Makefile to serveESP.make; add tufts for Nodes when displaying tessellation sesnitivities in serveESP; ensure that tessellation sensitivities on the boundaries of Face and ends of Edges agree with Edge and Node sensitivities; update .tsen files — jfdannen / detail
  65. remove serveESP from Makefile and NMakefile; create serveESP.make and serveESP.mak — jfdannen / detail
  66. index.lt.0 in RESTORE is now a wildcard, meaning that Bodys in all storages that match name are put onto the stack; initial integration of CAPS into serveESP — jfdannen / detail
  67. update ESP_QuickReference; update ESP so that double-clicking in MessageWindow takes you to filename:linenum that is clicked on, else takes you to last filename:linenum in MessageWindow; update sensCSM to use RMS errors; reduce ratio of errors for alerts from 5 to 2; update .gsen and .tsen files for 7.4.1 and 7.6.0 — jfdannen / detail
  68. Effective Topology: fix memory problem with interior Edge (in Loop twice) — haimes / detail
  69. Remove Degenerate EEdges once Faces are merged in an EFace — haimes / detail
  70. Documentation updates — haimes / detail
  71. Only consider Python 3.8 or newer in makeEnv/winEnv — galbramc / detail
  72. Some config updates — haimes / detail
  73. Export Body CSystems that ar free-standing (aLen = 9) — haimes / detail
  74. Update egads_dot.h. Explicit makefile in scan-build/scan-view. — galbramc / detail
  75. change startup messages in wvServer from stderr to stdout — jfdannen / detail
  76. Add missing DEGENERATE Edge check to EG_exportModel — galbramc / detail
  77. Properly reflect the move of OCC copies in this repository — haimes / detail
  78. Fix array out-of-bounds in EG_extrude_dot — galbramc / detail
  79. Update error messages — galbramc / detail
  80. Upgrade to PointwiseV18.5R1 on macOS — galbramc / detail
  81. Another python suppression — galbramc / detail
  82. Don't build CAPS for OpenCSM with lint — galbramc / detail
  83. Upgrade to Pointwise V18.5R1 on reynolds. Compile serveESP for stanalizer — galbramc / detail
  84. Build serveESP — galbramc / detail

#153 (Jan 8, 2022, 5:30:53 AM)

  1. add stl2tess to Slugs/src — jfdannen / detail
  2. Fix for last commit — galbramc / detail
  3. Add more checks when reading sens file — galbramc / detail
  4. Update the copyright year to 2022 — haimes / detail
  5. Add SensFile reader to fun3dAIM — galbramc / detail
  6. Fix memory leak — galbramc / detail
  7. fix ocsm.def for last commit — jfdannen / detail
  8. implement Debug in .csm editor; fix a few bugs in parent/child Body pointers — jfdannen / detail
  9. fix bugs associated with running a .py file — jfdannen / detail
  10. improve tessellation sensitivities for Edges and Faces that are not trimmed by a boolean operation; update .tsen files — jfdannen / detail
  11. update to special plotting associated with Face velocities for INTERP_VEL=6 — jfdannen / detail
  12. update CAPS/myPlane for intent 4 (avl) and 5 (bem) — jfdannen / detail
  13. .py files are automatically edited in serveESP — jfdannen / detail
  14. fix bug if user with ball loaded a new .csm file and then editted (where user without ball saw wrong file); adjust closeBody2 test case to REORDER the Edges so that all the Faces point outward; inform users if someone has joined or left a collaboration; modify visualization of grids on Faces in serveESP * new-style quads if _makeQuads attribute on Body (.tessType==Quad or .tessType==Mixed) * old-style quads if _makeQuads attribute on Faces (npatch>0) * triangles oherwise; fix valgrind error in sensCSM — jfdannen / detail
  15. check for proper result types from SWEEP; check mass props before topology in verification files; add verification data for OpenCASCADE 7.6.0; change copyright to 2022; give @version a value before the first Body is built; allow DIMENSIONed variable to be returned from a UDC; update expressions test case to demonstrate multi-value returns — jfdannen / detail
  16. more improvements to udfCompare (and fix sanitizer error); update GRAFIC plots of trim curves to show the points on the Pcurves; increase size of plot points in serveESP; add fixture* test cases — jfdannen / detail
  17. update Makefile.LINUX64; fix typo in mitten2.py — jfdannen / detail
  18. add ability of specifying color of items in a plotfile; change udfCompare to group bad points based upon the magnitude of the error; update udfCompare to check distances from Nodes and Edges as well as distances from Faces — jfdannen / detail
  19. update documentation for DUMP in OpenCSM.h — jfdannen / detail
  20. create ESP-viewer.js; pass callbacks to overlay if they exist; reorganize python test cases; change verification tests to include serveESP — jfdannen / detail
  21. update timViewer to be able to make changes to the scene graph; add ocsmLoadFromModel to pyOCSM; add shapes test case; fix initialization bug in ocsm.py — jfdannen / detail
  22. remove timViewModl from timPython.def — jfdannen / detail
  23. add ability to hold python execution while an overlay (viewer) is executing; update viewer.py — jfdannen / detail
  24. change locking mechanism for tim_mesg to use ThreadWait instead of a mutex — jfdannen / detail
  25. fix bug that did not release the lock if tim_mesg is in holding mode — jfdannen / detail
  26. fix mutex problems associated with tim_mesg — jfdannen / detail
  27. fix lint error in tim.c — jfdannen / detail
  28. add timViewer and associated test cases; create mechanism to allow a TIM to call a TIM; fix a few typos in ESP-help — jfdannen / detail
  29. updates to Windoze .def files — jfdannen / detail
  30. Fix up ocsm.def to include GetToken (whose name should probably change) — haimes / detail
  31. When splitting periodic Faces convert the order 1 BSpline to a Line; allow for calling wv_cleanupServers when there are active server threads — haimes / detail
  32. Don't fix Faces for EG_makeTopology(SOLIDBODY) — galbramc / detail
  33. Fix warning — galbramc / detail
  34. Makefile fixes — galbramc / detail
  35. Add EG_rotate_dot and fix EG_extrude_dot to work with OCC. Fix for topological changes in EG_rotate (not enabled yet). — galbramc / detail
  36. Do not use Edge tangents for vertex insertion when the magnitude is deminishing (near C0s) — haimes / detail
  37. Though premature, update the copyright year — haimes / detail
  38. Attempt to mitigate the fillTris error when presented with bad PCurves — haimes / detail
  39. Use the define OCC_VERSION_DEVELOPMENT to tell us we built the OpenCASCADE — haimes / detail
  40. Some trianglulation changes needed to deal with the GMGW3 'foreign' STEP files — haimes / detail
  41. Another config update and another update for Windows threading — haimes / detail
  42. Start FireFox in the background & try to fix stanalizer again — haimes / detail
  43. Try to fix a scan-build error that spuriously popped up! — haimes / detail
  44. Protect EMP against NULL Locks and revert Safari's startup in makeEnv — haimes / detail
  45. Windows config update and fix for EMP_LockTest for Windows — haimes / detail
  46. Fix arcLength (EGADSlite) for BSplines and update some configure scripts — haimes / detail
  47. Fix bugs in EG_arcLenSeg — galbramc / detail
  48. Add MCG_ARCLEN for suggested arclen approach. Add some more quadrature options for arclen. — galbramc / detail
  49. Temp brute-force fix for BSpline arcLength in EGADSlite; fix tessellation progress warning; config updates -- list XQuartz min rev — haimes / detail
  50. Add csmOCCFailure_7.6.txt file — galbramc / detail
  51. Remove csmOCCFailure_7.6.txt directory — galbramc / detail
  52. Should have kept last numpy suppression — galbramc / detail
  53. Another try at suppressing numpy leaks — galbramc / detail
  54. Suppress numpy leak — galbramc / detail
  55. mark closeBody2 as expected success — jfdannen / detail
  56. Update sens parsers — galbramc / detail
  57. Jenkins updates for 7.6.0 — galbramc / detail
  58. update verification data for 7.6.0 — jfdannen / detail
  59. Add Warning: LockRelease to jenkins parser — galbramc / detail
  60. Don't set UDUNITS for beta testing — galbramc / detail
  61. Remove last alias — galbramc / detail
  62. Alias python3 on centOS8 — galbramc / detail

#152 (Dec 11, 2021, 6:02:41 AM)

  1. Fix memory leaks — galbramc / detail
  2. Stanilizer fixes — galbramc / detail
  3. Turn on mses unit test again — galbramc / detail
  4. Fix warning and NMakefile — galbramc / detail
  5. Fixes for last commit. Dissable mses test until we get new executables. — galbramc / detail
  6. MSES seems to be working with sensitvities. Many other small fixes. — galbramc / detail
  7. Fix relative path for geometry.view and analysis.view functions — galbramc / detail
  8. Fix TACS test — galbramc / detail
  9. Fix TACS sensitivity bug — galbramc / detail
  10. First cut at MSES IO functions — haimes / detail
  11. outLevel fixes — galbramc / detail
  12. mses, xfoil and tacs fixes — galbramc / detail
  13. Added initial msesAIM. Very much work in progress — galbramc / detail
  14. Update awave error message — galbramc / detail
  15. Lint fix — galbramc / detail
  16. Dynamic Web Server library cleanup — haimes / detail
  17. Revert back to testing just one wing with cart3d — galbramc / detail
  18. Skiop inconsistency — galbramc / detail
  19. Tighter tolerance — galbramc / detail
  20. Enable multi-body for cart3d — galbramc / detail
  21. Changes due to WV server library going dynamic — haimes / detail
  22. Update getDeriv for GeometryOut for fattened Jacobian — galbramc / detail
  23. Allow for passing a user pointer into browserMessage -- breaks compatibility! — haimes / detail
  24. Fix sanitizer — galbramc / detail
  25. doc fix — galbramc / detail
  26. Warning fix — galbramc / detail
  27. lint fixes — galbramc / detail
  28. Forgot Windoze NMake — galbramc / detail
  29. Add TACS AIM — galbramc / detail
  30. Update pyCAPS examples to use caps_system call — galbramc / detail
  31. Cart3D fixes. Remove os.chdir from doc. Add some Analysis.system calls. — galbramc / detail
  32. Windoze tim.def file fix — jfdannen / detail
  33. Windoze tim.def file fix — jfdannen / detail
  34. Windoze ocsm.def file fix — jfdannen / detail
  35. fix bug in udfNuscale when zscale=0; first cut at calling a TIM from the python TIM — jfdannen / detail
  36. do not change outLevel if ocsmSetOutLevel is called with negative argument; add printout for outLevel.ge.2 when TIM broadcasts a message to the browsers; add option to plot trim curves associated with all Faces — jfdannen / detail
  37. fix Makefile and NMakefile associated with previous commit — jfdannen / detail
  38. reorganize way TIMs broadcast messages out to the browsers — jfdannen / detail
  39. use new BLEND C0 treatment, including in sensitivities; fix bug in timMitten associated with cleanup — jfdannen / detail
  40. remove special sensitivity treatment for BLENDs with C0; update sensCSM verification data; add blend22* and rule22 test cases — jfdannen / detail
  41. add closeBody* test cases; add collab* test cases; add documentation about what happens for a 0 in either SELECT ADD or SELECT SUB statement — jfdannen / detail
  42. fix bug associated with MALLOCs in udpWaffle — jfdannen / detail
  43. fix bug associated with MALLOCs in udpWaffle — jfdannen / detail
  44. fix typo in error message associated with COMBINE command; make error messages from udpWaffle more explicit — jfdannen / detail
  45. remove Windoze warnings from previous commit — jfdannen / detail
  46. fix bug so that wv_broadcastText is not called by a TIM when running in batch mode — jfdannen / detail
  47. improve sensitivities associated with tip treatments in BLEND;  update designG5 and designG6 test cases — jfdannen / detail
  48. update timMesg so that it does not return a response; allow non-*host* to initiate file editting; fix bug that does not give appropriate error messages when having trouble opening a .py file from the command line; fix bug that does not give appropriate error messages when timPython has trouble opening a file from Tools->Python — jfdannen / detail
  49. add IFTHEN/ENDIF to udpWaffle; add waffleA* test cases; update Makefile.LINUX64; fix bug that lost ball if user with ball sends a message — jfdannen / detail
  50. update to latest version of Tblade3; allow default xaxis=0, yaxis=0, and zaxis=0 in ROTATEX, ROTATEY, and ROTATEZ; add rotate1 test case; add LAYOUT to udpWaffle; add waffle3* test case; fix memory leak associated with mprp arrays — jfdannen / detail
  51. remove tim_data from TIM interface (no longer needed); change NXSECT argument in udpGuide to ATTRREAL so that finite difference sensitivities are triggered; remove RAD argument to udpBox in designI2 and designI4 so that finite differences are not triggered; fix bug that re-computed velocities when re-computation was not necessary; fix bug associated with sensitivities of mass properties — jfdannen / detail
  52. update timPython to post stdout to the MessageWindow while python is executing; fix postMessage so that blanks lines between messages are eliminated — jfdannen / detail
  53. fix bug introduced in previous commit — jfdannen / detail
  54. fix bug in previous commit; add ocsmLoadFromModel function; add withTess=0 argument to DUMP command; add dumpTest* test cases; update Makefile.DARWIN64.gcc — jfdannen / detail
  55. updates to way PYTHONHOME is found in timPython; add designV1 test case; set sgFocus in timMitten so that it can start with an empty MODL — jfdannen / detail
  56. fix hint for SCALE in ESP; shorten hint time-out in ESP to 15 sec; fix bug in timPython that did not register MesgCB and SizeCB after a timLoad (which caused a problem if started with a .py file); add (not-yet-implemented) Debug option in ESP editor; update ESP editor so that users without the ball can see the edits as they happen; clean up unneeded routines in TIMs; add debug and debugged testcases; add designU* test cases — jfdannen / detail
  57. For timPython, put Python lib location in rpath — haimes / detail
  58. Build timPython on Windows — haimes / detail
  59. Allow for building TIMs outside of OpenCSM/src — haimes / detail
  60. fix bug that did not remove UDP cache entries when a Body was deleted (such as during recycling); do not consider DEGENERATE Edges when SELECTing either by bounding box of nearest-to — jfdannen / detail
  61. Dynamic Web Server library cleanup — haimes / detail
  62. update Makefile.DARWIN64 to get rpath correct when running a script directly in python (and not serveESP); allow CONPMTR statement if at top-level scope; update scope test case; fix bug that did not compute sensitivity correctly when a DESPMTR was used in an include-type UDC — jfdannen / detail
  63. change button to yellow and RunningPython while python is running; update Makefile.DARWIN64 to get rpath correct when running a script directly in python (and not serveESP) — jfdannen / detail
  64. start getting ready for editing .py files; separate loading and running in timPython; run python in separate thread; add @scope at-parameter; add scope test case — jfdannen / detail
  65. add esp.h — jfdannen / detail
  66. fix stdout sequencing in timPython; get rid of global data in serveESP; update data/python/test* files — jfdannen / detail
  67. fix bug associated with (FD) sensitivities for a MODL with a OUTPMTR being set before a FILLET or CHAMFER; remove code that protects against calls to EG_replaceFaces in SolidBodys with more than one Shell — jfdannen / detail
  68. enable esp.SetModl() to see MODL in serveESP after python finishes — jfdannen / detail
  69. stub routines added to account for change in tim.def — jfdannen / detail
  70. fix bug associated with expressions that contain mass properties; add designT1 test case; fix bug that causes an error to be raised if the mass property sensitivities are computed by FDs before the last Body is built; add designT2 test case; add python mode to CodeMirror; add esp.py — jfdannen / detail
  71. Fix recent update to EGADSlite arcLength code — haimes / detail
  72. Fix the flattening of periodic BSplines when exporting to EGADSlite; Fix the BSPline range in EGADSlite when flattened — haimes / detail
  73. Fix STEP writing with Name attributes — haimes / detail
  74. Mark invalid Bodies during STEP/IGES reads & limit number of inserted verts when tessellating — haimes / detail
  75. Add MSVC 2022 config file — galbramc / detail
  76. C0 Blend fix so that the surface derivatives (normals) do not vanish at the C0 — haimes / detail
  77. Scan-build fix from last commit — haimes / detail
  78. Functions that insert verts based on surface arc-length use underlying Face when all positions are on that Face (Effective Topology). Used in Quadding & HO tessellations. — haimes / detail
  79. Improved error message. Use brep V2 for OCC 7.6. Minor cleanup. — galbramc / detail
  80. Minor cleanup — galbramc / detail
  81. Fix blend tip sensitvities — galbramc / detail
  82. clang analyzer fixes — galbramc / detail
  83. Fix so DEBUG code compiles — galbramc / detail
  84. Allow degenerate Edge in blend/ruled sections — galbramc / detail
  85. Fix Blend with non-planar sections and Degerenate Edges in the Sections — haimes / detail
  86. Target arm64 darwin, not specifically macOS 11 — galbramc / detail
  87. Reduce testing output — galbramc / detail
  88. Cleanup old config files — haimes / detail
  89. No longer support OCC less than Rev 7.0 — haimes / detail
  90. Changes to support OCC 7.6 — haimes / detail
  91. Update Makefiles for OCC 7.6. Improved error message. — galbramc / detail
  92. Other Python config updates — haimes / detail
  93. Dynamic Web Server library cleanup — haimes / detail
  94. Fix the new wvServer dynamic library def file — haimes / detail
  95. NMakefile fix — haimes / detail
  96. Make the WV server library dynamic — haimes / detail
  97. Allow for passing a user pointer into browserMessage -- breaks compatibility! — haimes / detail
  98. minor updates — haimes / detail
  99. update udpTire.c — jfdannen / detail
  100. Update mses path — galbramc / detail
  101. Add PATHs to MSES executables — galbramc / detail
  102. Upgrade Jenkins env for MSVC 2022 — galbramc / detail
  103. Update CASROOT — galbramc / detail
  104. Add 7.6 files (needs updates) — galbramc / detail
  105. Enable awave again for viggen — galbramc / detail
  106. mark closeBody2 as failure on 7.3.1 — jfdannen / detail
  107. changes to accomodate better sensitivities associated with tip treatment in BLENDs — jfdannen / detail
  108. Dissable astros and awave on viggen — galbramc / detail
  109. Use bash on Windoze as well — galbramc / detail
  110. Suppress output from test_pyOCSM to reduce email generation time. — galbramc / detail
  111. Update python suppression — galbramc / detail
  112. add debug.csm as known failure — jfdannen / detail
  113. Update how ESP_Python is found — galbramc / detail
  114. source .bashrc as part of env — galbramc / detail
  115. Fix scipy install on viggen — galbramc / detail

#148 (Oct 14, 2021, 12:41:32 PM)

  1. Relax tolerances and skip openmdao when not available — galbramc / detail
  2. More Cart3D fixes — galbramc / detail
  3. Suppress ocsmBuild for sensitvities (to avoid confusion with an actual build). Fix memory leak in xddm. — galbramc / detail
  4. Fix uninitialized data — galbramc / detail
  5. Cart3D fixes — galbramc / detail
  6. More xddm cleanup — galbramc / detail
  7. xddm fix — galbramc / detail
  8. scan-build fixes — galbramc / detail
  9. Lint fixes — galbramc / detail
  10. Minor cleanup of xddm — galbramc / detail
  11. Tweak cart makefile — galbramc / detail
  12. Update fun3d OpenMDAO example with dynamic outputs — galbramc / detail
  13. Warning fix — galbramc / detail
  14. Fix warning and properly suppress cart test when cart is not present — galbramc / detail
  15. Update training example — galbramc / detail
  16. Add auto execution for AnalysisDynO. Update Fun3D with AnalysisDynO. Add geometric sensitivties to Cart3D. Allow Cart3D trix file to write multiple bodies. Add Cart3D_OpenMDAO examples. Other smaller cleanup. — galbramc / detail
  17. Remove AFLR3 .tmp.b8.ugrid file — galbramc / detail
  18. Surpress delete message (for wildcards that don't exist) on Windows — haimes / detail
  19. add udpTire testing to Makefile.DARWIN64; compute sensitivity correctly with udfCreateBEM, udfCreatePoly, udfEditAttr, udfPrintBbox, udfPrintBrep, and udoPrintEgo; add designS3 test case — jfdannen / detail
  20. fix bug in udfGuide that did not work properly for an input SheetBody; add designS* test cases; fix bug that did not compute sensitivities correctly for UDFs when one (or more) of the input Bodys had non-zero velocities — jfdannen / detail
  21. Allow for Solids with multiple Shells to work in replaceFaces — haimes / detail
  22. update udpTire.c — jfdannen / detail
  23. update udpTire.c for change in cacheUdp — jfdannen / detail
  24. Another Cart3D suppression — galbramc / detail
  25. Typo — galbramc / detail
  26. Use anaconda scipy on M1 for now — galbramc / detail
  27. Don't use system site packages for virtualenv — galbramc / detail
  28. Don't use system site packages for virtualenv — galbramc / detail
  29. Suppress perl memory leaks — galbramc / detail
  30. Remove cart perl command — galbramc / detail
  31. Missed CART3D_ARCH — galbramc / detail
  32. Update env for Cart3D design — galbramc / detail

#146 (Oct 9, 2021, 4:27:00 PM)

  1. add udpBodyList to udp.def for Windoze — jfdannen / detail
  2. initial code for proper UDF sensitivity calculations --- not enabled yet — jfdannen / detail
  3. modify udfNacelle* test cases to use DESPMTRs; add splitWing* test cases; remove gallery/FJ2 from distribution — jfdannen / detail
  4. Retain all Topology labeled as internal but reject external -- keep Nodes that have a non-zero valence — haimes / detail

#145 (Oct 9, 2021, 6:10:56 AM)

  1. Fix Windows Warning — haimes / detail
  2. First complete implementation of Analysis Dynamic Outputs — haimes / detail
  3. Warning fix — galbramc / detail
  4. Add wildcards to aim_rmDir — galbramc / detail
  5. Initialize Dynamic Outs for Restart — haimes / detail
  6. Begin the implementation of Analysis Dynamic Outputs — haimes / detail
  7. Windoze fix — galbramc / detail
  8. Fix last commit — galbramc / detail
  9. Add aim_rmDir and update doc — galbramc / detail
  10. lint fix — galbramc / detail
  11. Fix for reading fun3d .forces with Design_Functional — galbramc / detail
  12. Fix last commit — galbramc / detail
  13. Lint fixes — galbramc / detail
  14. Fix copy/pasta mistake — galbramc / detail
  15. Add interferenceTest to regressoin testing — galbramc / detail
  16. Minor corrections to interfaceAIM — haimes / detail
  17. scan-build fix and dox fix — galbramc / detail
  18. First cut at the Interference AIM and fix for instance number — haimes / detail
  19. Fix AIM_ERROR that was just supposed to be printf — galbramc / detail
  20. scan-build fixes — galbramc / detail
  21. Missed a file — galbramc / detail
  22. Remove static array size for fun3d Design_Functional output. Return list of list for matrix varaibles. — galbramc / detail
  23. Remove fun3d redundant sensitivity calculations — galbramc / detail
  24. Fix some typos — galbramc / detail
  25. Fix bug in limits check in caps_setValue — galbramc / detail
  26. Forgot a file — galbramc / detail
  27. Fixes to fun3d design — galbramc / detail
  28. Another warning fix — galbramc / detail
  29. Fix warning — galbramc / detail
  30. Forgot a file — galbramc / detail
  31. Re-design fun3d design. — galbramc / detail
  32. fixes for the previous fixes — ryan / detail
  33. fix for fun objectives/functionals — ryan / detail
  34. Fix for cart bc — ryan / detail
  35. Lint fix — galbramc / detail
  36. Fix vlmControl memory leak — galbramc / detail
  37. Make avl smarter about when geometry file needs to be written. — galbramc / detail
  38. UpRev to 1.21 — haimes / detail
  39. Fix documentation links — galbramc / detail
  40. another dox fix — ryan / detail
  41. minor dox fix last one — ryan / detail
  42. minor dox fix — ryan / detail
  43. minor dox fix — ryan / detail
  44. Clear aim_setSensivity flags in caps_build — galbramc / detail
  45. Remove spurious Info printf and make sure we are between pre and post in caps_system — haimes / detail
  46. Document AIM execution modes — haimes / detail
  47. Don't use packaging module — galbramc / detail
  48. Fix last commit — galbramc / detail
  49. Try multiple times to get pointwise license in C example — galbramc / detail
  50. Check OpenMDAO version in example — galbramc / detail
  51. allow pyOCSM to import an externally-generated MODL; enable timMitten; create timPython and data/python/test*; allow serveESP to start with a .py file — jfdannen / detail
  52. uprev to v1.21 (ESP.html, ESP-localhost.html, ESP-help.html, OpenCSM.h, ocsm.rc); create unpublished routine convertToBSplines.  modify so that BSPLINES are not reconverted; modify udpParabaloid and udfNuscale to use convertToBSplines.  add DEBUG prints in udpNuscale; add sculpt* test cases; add attrProp test case to show how Edge Attributes get propagated to Faces by EXTRUDE and REVOLVE but not by RULE and BLEND.  attrProp also shows that Edge Attributes are not carried forward by EXTRUDE, REVOLVE, RULE, or BLEND; make INTERFACE DIM obsolete; issue WARNING to use DIMENSION statement instead; add udfNacelle and udfNacelle* test cases — jfdannen / detail
  53. fix bug in StepThru that was introduced when computation of _tParams was delayed — jfdannen / detail
  54. Revert Internal Node change -- needs more work! — haimes / detail
  55. An update to the last commit -- remove only Nodes that do not touch active Edges — haimes / detail
  56. Retain Internal Nodes during parsing -- THIS HAS THE POTENTAL OF BREAKING EGADS FILE READS FOR THE VERY FEW CASES THIS COMES UP (the result is that attributes are not reassigned) — haimes / detail
  57. scan-build suppression — galbramc / detail
  58. Use jenkins setup for scan-build instead — galbramc / detail
  59. Removed one too many — galbramc / detail
  60. Another M1 fix — galbramc / detail
  61. Fix last commit — galbramc / detail
  62. Remove depreacted OCC defines. Fix for M1 scan-build. — galbramc / detail
  63. Enable scan-build on M1 — galbramc / detail
  64. Fix analyzer error — haimes / detail
  65. UpRev to 1.21 and deal with Names on STEP IO — haimes / detail
  66. README updates — haimes / detail
  67. Change back the getdelim suppression — galbramc / detail
  68. Fix bug in Valgrind script not running CAPS unit test — galbramc / detail
  69. Try again to suppress libc problem — galbramc / detail
  70. Adjust sanitizer suppression — galbramc / detail
  71. Upgrade macOS clang 13.0.0 — galbramc / detail
  72. remove testing of plugs7a and plugs8a — jfdannen / detail
  73. Another Windoze fix — galbramc / detail
  74. Turn on echo again for Windoze — galbramc / detail
  75. Fix typo — galbramc / detail
  76. Add M1 to PreBuilt testing — galbramc / detail

#144 (Sep 15, 2021, 10:36:31 PM)

  1. Add caps_system to pyCAPS — galbramc / detail
  2. Add caps_system and add relative path to aim_system — haimes / detail
  3. Upgrade to SU2 7.2.0. Minor avl fixes. — galbramc / detail
  4. Try again and again for pointwise license — galbramc / detail
  5. Adjust regression script — galbramc / detail
  6. Suppress lint warnings — galbramc / detail
  7. Fix warning — galbramc / detail
  8. Cleanup hasDeriv/getDeriv and add derivatives outputs to avlAIM — galbramc / detail
  9. Update the documentation and make the AIM file function naming consistent — haimes / detail
  10. Add caps_journalState to the API — haimes / detail
  11. Upgrade to SU2 7.2.0 — galbramc / detail
  12. Fix for memcheck python script on darwin — galbramc / detail
  13. scipy and curses suppressions — galbramc / detail
  14. Use system site-packages — galbramc / detail
  15. Add openmdao to pip installs — galbramc / detail

#143 (Sep 11, 2021, 5:53:12 AM)

  1. Lint fix — galbramc / detail
  2. Fix warning — galbramc / detail
  3. Forgot the actual journal test — galbramc / detail
  4. First test of journaling. Works with autoExect, manual exec needs work. — galbramc / detail
  5. Replace verbosity with outLevel — galbramc / detail
  6. Use environment variable to set outLevel=0 for CAPS training (not documented). — galbramc / detail
  7. fix minor documentation typos — jfdannen / detail
  8. fix bug in JOIN associated with Faces that share Edges but which have bounding boxes that differ; add halfSpheres test case; remove spurious printout associated with SKEND command; disable timMitten for v1.20 release — jfdannen / detail
  9. fix bug that did not allow a user to EXTRACT multiple Faces from a SolidBody that had an internal hole; add extractShells* test cases; fix memory leak in udpParabaloid; update ESP_QuickReference for v1.20 — jfdannen / detail
  10. Remark on a fix to pyEGADS — galbramc / detail
  11. Fix memory problem in patch tester — haimes / detail
  12. Fix last warnings — haimes / detail
  13. First cut at the Tessellation support for the new AVL; some configuration updates — haimes / detail
  14. Globaly have minimal CAPS output for Jenkins — galbramc / detail
  15. Upgrade pointwise on macys — galbramc / detail
  16. Generalize valgrind suppression — galbramc / detail
  17. Another python suppression — galbramc / detail
  18. remove plugs6a from testing — jfdannen / detail
  19. Add keep-debuginfo=yes to valgrind flags — galbramc / detail
  20. Another attempt at suppression Pointwise memory problems — galbramc / detail
  21. Always suppress alloc dealloc missmatch — galbramc / detail
  22. Undo mystran exclusion from valgrind — galbramc / detail

#142 (Sep 5, 2021, 6:11:55 PM)

  1. Add awave cCAPS example. Cleanup avl AutoSpan test. Return CAPS_EXECERR error on system call failure. — galbramc / detail
  2. Comment out matplotlib import — galbramc / detail
  3. update config/README.txt for v1.20 — jfdannen / detail
  4. Suppress Pointwise alloc dealloc missmatch — galbramc / detail
  5. Exclude mystran from valgrind due to false positives — galbramc / detail
  6. Upgrade Pointwise on Reynolds — galbramc / detail
  7. Undo last change — galbramc / detail
  8. Fix virtualenv on macOS — galbramc / detail
  9. Don't install matplotlib for testing — galbramc / detail
  10. Upgrade valgrind version — galbramc / detail

#140 (Sep 3, 2021, 5:41:19 PM)

  1. add attrScribe test case; add udpParabaloid and test case — jfdannen / detail
  2. implement outLevel in timPlugs; improve -plugs option in serveCSM so that more cases work; mark plugs4c as successful in Jenkins; modify attribute and getattr1 test cases to account for fact that _tParams only exists when Body is recycled — jfdannen / detail
  3. Dissable undefined sanitizer vptr option — galbramc / detail
  4. mark plugs4c as successful — jfdannen / detail

#139 (Sep 3, 2021, 10:50:43 AM)

  1. fund3dAIM fixes — galbramc / detail
  2. Append ESPFLIBS on M1 — galbramc / detail
  3. Update mastros rpath on viggen — galbramc / detail
  4. Update ESPFLIBS for ubsan — galbramc / detail
  5. Another viggen update — galbramc / detail