add transform option to Slugs; in Slugs, the KutTriangles option automatically creates a link; add optional GRAFIC to cutTriangles in Slugs; lots of improvements to Slugs, especially for cases with very noisy .stl file
(detail)
by jfdannen
fix bug that did not allow interactive use of Slugs; add Kut triangles option in Slugs; update testScript.txt to include interactive Slugs testing
(detail)
by jfdannen
Fix surface to volume map for MultiDomain tetgen meshes
(detail)
by galbramc
Add geometric sensitvities to masstran AIM
(detail)
by galbramc
Suppress 'could not create compact unwind on M1' for now.
(detail)
by galbramc
Allow for FD sensitivities within the AIM sensitivity functions
(detail)
by haimes
add torch tests for kulfan corsairlite examples
(detail)
by dongjoon
update problem.__init__ to work if running pyscript but not from capsMode
(detail)
by jfdannen
Simplify multifidelity examples. Fig unit bug in caps.py. Support pint dimensionless to caps.py. Support udunits import in corsairlite.
(detail)
by galbramc
// 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; }
*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;
// 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); }
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; }
/* 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;
// 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;
// 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++) {
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++) {
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);
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); }
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; } }
// 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; }
// 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) {
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++) {
#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++) {
// 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; }
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);
// 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;
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);
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];
// 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;
// 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; }
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
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
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
// 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
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;
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);
// 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; } }
// 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);
// 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);
// 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);
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;
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
// 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);
// 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;
Fix typo in doc, and short Python Units Appendix
(detail)
by galbramc
Minor documentation (and coding) changes for the Python 1-to-1 CAPS interface
(detail)
by haimes
Document 1-to-1 pyCAPS API. Change pyCAPS setLimits/getLimits to setLimitsSize/getLimitsSize. Add support for limit arrays to pyCAPS.
(detail)
by galbramc
Addition of a CAPS API function that anly returns a Value Object's size
(detail)
by haimes
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.
(detail)
by galbramc
Fix childByName when in debug mode
(detail)
by haimes
Finally the correct Windoze executable
(detail)
by galbramc
Update avl exec on Linux and Windoze
(detail)
by galbramc
Fix both the documentation and the code so they are consistent with adding const to the AIM functions Pre and Exec
(detail)
by haimes
update jfd3AIM to adjust for recent API changes
(detail)
by jfdannen
Reinstate no intent phrase flag
(detail)
by haimes
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)
(detail)
by haimes
SU2 AIM updated for SU2-7.3.0, Changed boundary markers from numbers to strings
(detail)
by nitin
SU2 AIM updated for SU2-7.3.0, Changed boundary markers from numbers to strings
(detail)
by nitin
Add funID to caps_jrnlWrite. Fix call to caps_ownerInfo. Fix memory leak in egadsTess.
(detail)
by galbramc
Update pyCAPS for new history tracking. Fixes in su2. Test Phasing with egads/tetgen/su2.
(detail)
by galbramc
Move CAPS towards the new history tracking
(detail)
by haimes
Revert changes in CAPS/src from previous accidental commit
(detail)
by galbramc
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.
(detail)
by galbramc
Remove caps_freeOwner from pyCAPS. Add commented out caps_getHistory and caps_intentPhrase for the future.
(detail)
by galbramc
Put caps_freeOwner back in the def file
(detail)
by haimes
fixes for 'Flow' directory switching for fun3dAIM
(detail)
by ryan
Update flowchart viewing and check for parameter value object uniqueness
(detail)
by haimes
Track closure of problem object to avoid calling caps_delete. This will leak attributes for now if the problem object is explicitly closed.
(detail)
by galbramc
Track closure of problem object to avoid calling caps_delete. This will leak attributes for now if the problem object is explicitly closed.
(detail)
by galbramc
Upgrade pointwimse AIM to support Pointwise 18.5. Rotate airfoils for astros. Add error checking for capsReference quantities in avl.
(detail)
by galbramc
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.
(detail)
by galbramc
Remove duplicated CFD meshes on disk. Lots of refactoring...
(detail)
by galbramc
fixed uninitialized variable in Slugs
(detail)
by jfdannen
First cut at documenting mesh writers
(detail)
by haimes
fix addressing bug in udpWarp; initial implementation of timVspSetup
(detail)
by jfdannen
convert backslashes to forward slashes in filenames in udpVsp3 because anglescipt expects forward slashes
(detail)
by jfdannen
properly rename vsp3 verification test cases
(detail)
by jfdannen
ATTRRECYCLE forces UDPRIM to not be recycled; fix printing of SELECT with attrValu a number; update udpVsp3 to automatically run vsp; add data/vsp3 test cases; add propInteg*a test cases; use _plugsIgnore on Faces that should be ignored during cloud point reclassifications in Plugs; add udpWarp; add ablate0* test cases
(detail)
by jfdannen
modify ocsmPrintBrep to show list Edges as degenerate rather than skipping them; add SPECIAL printCPs command to print the control points associated with all Nurbs in the Body on the top of the stack; add first version of udpVsp3
(detail)
by jfdannen
fix bug that did not retain Edge attributes during an ELEVATE operation
(detail)
by jfdannen
add mincp argument to udpNuscale and convertToBSplines functions; add ablate* test cases, including using Plugs; extend Plugs to be able to use DIMENSIONed DESPMTRs; fix bug in velocityOfEdge for Edges with dots that have changing tranges
(detail)
by jfdannen
temporarily remove vspFuse* test cases
(detail)
by jfdannen
analytic sensitivities for sketches; update .gsen and .tsen files; move sketch13* test cases to designY*
(detail)
by jfdannen
start implementing analytic sensitivities for Sketches; add sketch13* test cases
(detail)
by jfdannen
initial, very early, preliminary version of udpVspFuse and associated test cases
(detail)
by jfdannen
fix error in serveESP.mak in previous commit
(detail)
by jfdannen
fix error in serveESP.mak; fix debug code error for windoze
(detail)
by jfdannen
change the SUBTRACTion of two co-planar SheetBodys such that if the Loop associated with the second Body does not cut a clean hole in the first Body, revert to the (default) scribing behavior for SheetBodys, that is scribe the first Body with the Edges from the second Body; update the documentation for SUBTRACT; add scribe11* test cases; fix occasional bug associated with analytical sensitivities when DUMPing a .sens file from with .csm script
(detail)
by jfdannen
add support for BLENDs with .C1side Attributes; add blend23* test cases; change serveESP.mak so that it does not error if PYTHONINC is not set
(detail)
by jfdannen
uprev to v1.24 (ESP.html, ESP-localhost.html, ESP-help.html, OpenCSM.h, ocsm.rc)
(detail)
by jfdannen
added ability to double-click on python error to open pyscript editor at the correct line
(detail)
by jfdannen
add ocsmTraceAttrs, ocsmTraceStors, ocsmGetFiletree, and ocsmShowTblOfContents; add ability to trace Attributes, Storges, the FileTree, and the provided representations in the .csm editor; limit number of files in File menu in ESP so that it does not overflow the display; add SPECIAL provides statement to work with ocsmShowTblOfContents
(detail)
by jfdannen
fix bug that did not properly decrement the pattern counter for ENDIF statements after a PATBREAK was executed; fix bug that returned error after SELECTing a degenerate Edge
(detail)
by jfdannen
take shortcut when creating tessellation sensitivities for Faces that do not have Edges that are the result of Boolean operations; update .tsen files for both 7.6.0 and 7.7.0
(detail)
by jfdannen
fix typo in text string in serveESP
(detail)
by jfdannen
fix bug that copied Csystems into a UDPRIM; add udfMechanism and mechanism* test cases
(detail)
by jfdannen
add udfDeform to change BSpline control points associated with Faces (and Edges and Nodes) in a Body; add deform* test cases; modify some verification data to account for differences between 7.6.0 and 7.7.0 on all operating systems
(detail)
by jfdannen
remove GLOVES from ESP menu (for now)
(detail)
by jfdannen
update data/capsExamples/verify_7.7.0 files; use EG_copyObject since EGADS now allows the copying of PCurves; fix typos in testScript.txt; updates to avoid warnings in latest version of clang
(detail)
by jfdannen
allow cases with Edge-only tessellations to work; add edgeGrid* test cases; update 7.7.0 verification data for imprint5 (should fail); add PLOT_DISPLACE to plot displaced Edges based upon tessellation sensitivities
(detail)
by jfdannen
fix recently-added bug associated with FD tessellation sensitivities for NodeBodys
(detail)
by jfdannen
create verification data for OCC 7.7.0; make sure EVALUATE returns dot=0
(detail)
by jfdannen
compute tessellatiopn sensitivities correctly for a NodeBody; set and report sensitivity range properly when looking at only one component of tessellation sensitivities; add -forceFDs command line option to force finite difference sensitivities in serveESP
(detail)
by jfdannen
further improvements to flends; add flend10* test cases; fix memory leak when udfOffset encounters an error
(detail)
by jfdannen
add code to protect for UDFs that have references to its input Model; add (hidden) option to create timingReport.txt file; improve flends; add flend9* test cases
(detail)
by jfdannen
update flend8* test cases and their verification data
(detail)
by jfdannen
fix lint error in udfMatchBodys; re-add 7.4.1 verification data for flend8* test cases
(detail)
by jfdannen
change flend fits to be better in some cases; track down memory leaks
(detail)
by jfdannen
fix problem erroneously flagged by lint in udpOffset; update errors reported by udfFlend; fix bug where wrong UDP storage was used in udfOffset (which caused problems when udpOffset was called more than once); remove SPECIAL test of EG_otherCurve; document .multiNode use in RULE command; fix bug where UDP storage is incorrectly cleared for an embedded UDP (in udpCsm); fix bug in udfFlend that sometimes did not correctly reverse one of the Loops; fix memory leak in EVALUATE EDGEKT, EDGECP, FACEUKT, FACEVKT, and FACECP when a signal is raised
(detail)
by jfdannen
fix compiler warning in latest commit
(detail)
by jfdannen
account for fact that Py_SetPythonHome was deprecated in Python 3.11; update udfOffset so that it can generate offsets on SolidBodys (for FLENDs); add offset(I-R)* test cases; use zebra stripe coloring for maximum curvature in ESP (displayType=4); update convertToBsplines to reuse Pcurves whenever possible; add udpBspline and bspline* test cases
(detail)
by jfdannen
improve printing of Pcurves in ocsmPrintEgo; improve marking of Faces with _body and _brch Attributes when returned from a UDF; add PRINT_MAPPING_INFO define to OpenCSM; fix bug that did not account for Node and Edge sensitivities when report sensitivity range in serveESP; fix bug that occasionally did not properly compute Node tessellation sensitivities
(detail)
by jfdannen
add EVALUATE EDGEKT, EVALUATE EDGECP, EVALUATE FACEUKT, EVALUATE FACEVKT, and EVALUATE FACECP; add evaluate3* test cases; add multiRegions* test cases; add SELECT FLIP to flip the order of rntities in at-sellist
(detail)
by jfdannen
add copyAttr=0 argument to RULE and BLEND to copy Edge Attributes in xsects onto RULE/BLEND Edges; add blend29 and rule29 test cases
(detail)
by jfdannen
improve test for matching Nodes in udfMatchBodys; adjust matchBodys5 test case to check at-at-parameters
(detail)
by jfdannen
add arcLengths to ocsmPrintEgo; fix bug in udfMatchBodys when tol was large
(detail)
by jfdannen
modify udfMatchBodys to transfer an Attribute from one Body to the other when the Nodes, Edges, or Faces match; add matchBody5 test case
(detail)
by jfdannen
fix bug associated with UDP/UDF outputs that change shape
(detail)
by jfdannen
remove stanalyzer false positive by disabling checks for one loop in udfTester1.c
(detail)
by jfdannen
attempt to fix stanalizer error in udfTester1.c
(detail)
by jfdannen
allow at-at-parameters from a UDP/UDF to have velocities; split udpTest into udfTester1 (to demonstrate matrix outputs from a UDP/UDF) and udpTester2 (to demonstrate output velocities from a UDP/UDF); adjust Makefiles and associated test cases
(detail)
by jfdannen
add test.so to Makefiles; add diagnostic printouts when a Face does not tessellate; make sure all UDP/UDFs return output values to the numUdp instance
(detail)
by jfdannen
allow number of rows and columns to be accessible in UDPs; allow at-at-parameters from a UDP/UDF to be vectors or matricies; add udp_error1 to udp_error9 as possible argument to CATBEG statement; fix bug that did not restore at-at-parameters when a UDP/UDF was recycled; refactor way UDP/UDF arguments are handled
(detail)
by jfdannen
update header files to allow use of C++ main program; add simpleCSM.cpp and modify Makefile.DARWIN64 for testing
(detail)
by jfdannen
change udfSlices to skip a slice if OpenCASCADE fails (rather than return an eror); add SET statement to udpWaffle; add waffleB1 and waffleB2 test cases; fix bug in plotting of normals in serveESP and serveCSM; add PRINT_RULE_KNOTS to print number of knots in BSPLINEs entering a RULE or BLEND; allow -ATTRREALSENS to indicate a UDPRIM parameter that is intended to signal that Parameter is used within an included file; add waffleB3 test case to demonstrate this; add waffleB3.jrnl to demonstrate generation of .sens file for every DESPMTR without rebuilds
(detail)
by jfdannen
remove training/data/session10/version99/*
(detail)
by jfdannen
uprev to v1.23 (ESP.html, ESP-localhost.html, ESP-help.html, OpenCSM.h, ocsm.rc)
(detail)
by jfdannen
put 1 second delay when opening a pyscript to make sure that the browser is ready to display the pyscript in the editor; do not change status of ESP if an error was found while executing a pyscript; update session10 training files to be consistent with new method for creating a large, multi-component, multi-discipline model; update Makefile.DARWIN64; update ESP_QuickReference; update testScript.txt
(detail)
by jfdannen
write and test viewSimplySupport for transport test case
(detail)
by jfdannen
more updates to training/session10 and CAPS/training/python
(detail)
by jfdannen
more updates to transport training files
(detail)
by jfdannen
update error handling in flapz.udc; first attempt at python files for new transport for CAPS training
(detail)
by jfdannen
account for forward- and backward-slashes when opening a .udc
(detail)
by jfdannen
update ereped6a and project verification data again; fix recently-added bug that causes a memory fault for some non-manifold Edges
(detail)
by jfdannen
add Trace button in .csm editor to trace all top-level Parameters; reorganize session10 files to make it easier to build large, multi-component, multi-view models; update ereped6a, hollowC6, hollowC7, hollowCA, hollowCB, and project 3 verification data on 7.6.0 (because of issues highlighted by Xcode 14 on M1); update ereped6a and project3 verification data on 7.4.1 (because of issues highlighted by Xcode 14 on M1); add description of -dxdd flag to ESP-help
(detail)
by jfdannen
update testScript.txt; fix possible buffer overflow in SPECIAL command
(detail)
by jfdannen
update verification data for clearance1 and clerance4 to allow for slight differenes due to chaotic behavior of the EGADS tessellator
(detail)
by jfdannen
convert all (deprecated) sprintf calls to snprintf
(detail)
by jfdannen
update Verification.sh to use serveESP
(detail)
by jfdannen
fix bug where result of SPECIAL/clearance was not recycled
(detail)
by jfdannen
update .tsen files for designG5 and designG6
(detail)
by jfdannen
remove suppressions of sensitivity checks on designH2, designJ1c, designJ1d, and designJ5; add suppressions of sensitivity checks to designG5 and designG6 (for rounded tips); fix valgrind error in sensitivities in udpWaffle
(detail)
by jfdannen
fix some compiler warnings in most recent commit
(detail)
by jfdannen
add ocsmClearance; add clearance* test cases
(detail)
by jfdannen
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
(detail)
by jfdannen
add ocsmAdjoint and testOcsmAdjoint (in serveESP); allow timPlotter to have two y-axes; add plotter3.py test case
(detail)
by jfdannen
update adjoint test to remove (incorrect) compiler warning
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
clarify description of SELECT x -1 ibody1
(detail)
by jfdannen
update to account for latest EGADS fix to EG_generalBoolean
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
Remove dependencies from file copies
(detail)
by haimes
clean up usage of include files; update dependencies in Makefile and NMakefile
(detail)
by jfdannen
remove pyscript from lint in serveESP.make
(detail)
by jfdannen
update serveCSM.make to fix stanilizer error for pyscript
(detail)
by jfdannen
another attempt to move building of pyscript from Makefile to serveESP.make
(detail)
by jfdannen
move building of pyscript from Makefile to serveESP.make
(detail)
by jfdannen
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
(detail)
by jfdannen
next version of udfDroop2; first version of udfDroop3
(detail)
by jfdannen
next (still incomplete) version of udfDroop2
(detail)
by jfdannen
fix test_pyOCSM to account for updated MESSAGE statement
(detail)
by jfdannen
extend MESSAGE to optionally write to a file; modify message test case; update ESP-help and ESP_QuickReference
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
add assert in udfOffset to eliminate warnings on gcc
(detail)
by jfdannen
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
(detail)
by jfdannen
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)
(detail)
by jfdannen
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)
(detail)
by jfdannen
temporarily remove offset9* test cases
(detail)
by jfdannen
update udfOffset for negative distance
(detail)
by jfdannen
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)
(detail)
by jfdannen
add (unneeded) initializations in udfOffset to resolve incorrect lint complaints
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
initialize *info in ocsmBodyDetails; return sensible error message when a SET is applied to a CONPMTR, CFGPMTR, or DESPMTR
(detail)
by jfdannen
return stub.py to expected state at beginning of tutorial6
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
remove erroneous non-manifold warning associated with NodeBodys
(detail)
by jfdannen
fix bug when visualizing new-style quads and .mixed=NULL; do not process _makeQuads on a Face if the Body has a _makeQuads Attribute
(detail)
by jfdannen
fix bug in ESP.js that did not properly clear the display when velocities were cleared
(detail)
by jfdannen
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
(detail)
by jfdannen
fix bug in ESP.js associated with clearing velocities for a scalar; update tutorials 1 to 5 for v1.21
(detail)
by jfdannen
modify way python is initialized and finalized to work around a problem in numpy - take 2
(detail)
by jfdannen
modify way python is initialized and finalized to work around a problem in numpy
(detail)
by jfdannen
copy CFGPMTRs and DESPMTRs from temporary CAPS MODL back into the ESP MODL in CapsCommit and CapsSuspend; automatically rebuild MODL in CapsCommit and CapsSuspend
(detail)
by jfdannen
fix recently-added typo to timPyscript.c
(detail)
by jfdannen
add sanitize-thread target to Makefile.DARWIN64
(detail)
by jfdannen
fix error message associated with not enough arguments to SPHERE; fix bug when user double-clicks in MessageWindow (Chrome and Safari)
(detail)
by jfdannen
fix setting of PYTHONHOME for a double-quoted string
(detail)
by haimes
generate WARNING if importing an .egads file that has an Attribute that cannot be evaluated
(detail)
by jfdannen
add sub03 and f104 to studentGallery
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
compile serveESP in Makefile and NMakefile
(detail)
by jfdannen
inform user when CAPS overrides .csm value(s); first attemp at re-running pyscript in continuation mode
(detail)
by jfdannen
add indication for Phases that are a predecessor of current Phase in ListPhases
(detail)
by jfdannen
add calls to caps_debug around List* messages in timCapsMode; update data/ide/*.py to directly rebuild the geometry; reorder Caps menu options
(detail)
by jfdannen
additional bugs fixed in timCapsMode; add unequalRules test case to show how to RULE between xsects with different nnumber of Edges
(detail)
by jfdannen
rename capsIde to capsMode; restrict SPHEREs to have a positive radius; allow multiple ESP-help tabs in browser; partial implementation of timCapsMode
(detail)
by jfdannen
remove unneeded windoze headers from tim.c
(detail)
by jfdannen
Fix continuation marker in serveESP.mak
(detail)
by haimes
Fix headers for Windows in timPyscript
(detail)
by haimes
fix windoze problem with includes in pervious commit
(detail)
by jfdannen
only allow one open ESP-help tab at a time; implement timFlowchart; modify interpolateUpper.py and conserveUpper.py to demonstrate flowcharting tool
(detail)
by jfdannen
allow for .py scripts to run both in serveESP -caps and in python; add ocsmSetAuxPtr and ocsmGetAuxPtr in OpenCSM API
(detail)
by jfdannen
fix stanalyzer problem in timCapsIde.c
(detail)
by jfdannen
update serveESP.mak in previous commit to include timCapsIde.obj
(detail)
by jfdannen
change display scaling so that if there is only one NodeBody it still gets displayed; initial implementation of timCapsIde
(detail)
by jfdannen
allow nesting of TIMs to MAX_TIM_NESTING=10
(detail)
by jfdannen
fix bugs in File->Open and File->Edit for cases with and without UDCs
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
remove unneeded code that generates a clang warning
(detail)
by jfdannen
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
(detail)
by jfdannen
more improvements for ListHistory; use caps_open in read-only mode in ListPhases and ListHistory
(detail)
by jfdannen
do not let esp.UpdateESP rebuild (force rebuild in .py script)
(detail)
by jfdannen
extend ListHistory to GEOMETRYIN and GEOMETRYOUT
(detail)
by jfdannen
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
(detail)
by jfdannen
Makefile updates to supress Python builds
(detail)
by haimes
first implementation of CAPS menu in ESP
(detail)
by jfdannen
allow branching of Phases when running serveESP in -caps mode
(detail)
by jfdannen
caps_PrintObjects signature changed
(detail)
by haimes
allow wild-carding (with 0) in SELECT ADD and SELECT SUB; update connect9a and connect9b; update Makefile.DARWIN64 for 3-way testing
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
in timViewer, use a tessellation associated with the AIM if available, otherwise share the tessellation with OpenCSM
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
improve error message when project.csm exists before the project directorry when running serveESP in -caps mode
(detail)
by jfdannen
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
(detail)
by jfdannen
added numpts argument to udpKulfan, which helps avoid small dimples near the leading edge of a very thin section
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
for a UDP/UDF, only use finite differences if one of its arguements is ATTRREAL and has non-zero dots
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
remove serveESP from Makefile and NMakefile; create serveESP.make and serveESP.mak
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
update Makefile.LINUX64; fix typo in mitten2.py
(detail)
by jfdannen
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
(detail)
by jfdannen
update documentation for DUMP in OpenCSM.h
(detail)
by jfdannen
create ESP-viewer.js; pass callbacks to overlay if they exist; reorganize python test cases; change verification tests to include serveESP
(detail)
by jfdannen
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
(detail)
by jfdannen
remove timViewModl from timPython.def
(detail)
by jfdannen
update test.py to remove call to removed esp.ViewModl()
(detail)
by jfdannen
add ability to hold python execution while an overlay (viewer) is executing; update viewer.py
(detail)
by jfdannen
change locking mechanism for tim_mesg to use ThreadWait instead of a mutex
(detail)
by jfdannen
fix bug that did not release the lock if tim_mesg is in holding mode
(detail)
by jfdannen
fix mutex problems associated with tim_mesg
(detail)
by jfdannen
fix bug in udfNuscale when zscale=0; first cut at calling a TIM from the python TIM
(detail)
by jfdannen
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
(detail)
by jfdannen
fix Makefile and NMakefile associated with previous commit
(detail)
by jfdannen
reorganize way TIMs broadcast messages out to the browsers
(detail)
by jfdannen
use new BLEND C0 treatment, including in sensitivities; fix bug in timMitten associated with cleanup
(detail)
by jfdannen
remove special sensitivity treatment for BLENDs with C0; update sensCSM verification data; add blend22* and rule22 test cases
(detail)
by jfdannen
add closeBody* test cases; add collab* test cases; add documentation about what happens for a 0 in either SELECT ADD or SELECT SUB statement
(detail)
by jfdannen
fix bug associated with MALLOCs in udpWaffle
(detail)
by jfdannen
fix bug associated with MALLOCs in udpWaffle
(detail)
by jfdannen
fix typo in error message associated with COMBINE command; make error messages from udpWaffle more explicit
(detail)
by jfdannen
remove Windoze warnings from previous commit
(detail)
by jfdannen
fix bug so that wv_broadcastText is not called by a TIM when running in batch mode
(detail)
by jfdannen
improve sensitivities associated with tip treatments in BLEND; update designG5 and designG6 test cases
(detail)
by jfdannen
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
(detail)
by jfdannen
add IFTHEN/ENDIF to udpWaffle; add waffleA* test cases; update Makefile.LINUX64; fix bug that lost ball if user with ball sends a message
(detail)
by jfdannen
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
(detail)
by jfdannen
missed designU1c.tsen in previous commit
(detail)
by jfdannen
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
(detail)
by jfdannen
update timPython to post stdout to the MessageWindow while python is executing; fix postMessage so that blanks lines between messages are eliminated
(detail)
by jfdannen
fix bug introduced in previous commit
(detail)
by jfdannen
fix bug in previous commit; add ocsmLoadFromModel function; add withTess=0 argument to DUMP command; add dumpTest* test cases; update Makefile.DARWIN64.gcc
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
For timPython, put Python lib location in rpath
(detail)
by haimes
Allow for building TIMs outside of OpenCSM/src
(detail)
by haimes
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
(detail)
by jfdannen
Dynamic Web Server library cleanup
(detail)
by haimes
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
(detail)
by jfdannen
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)
(detail)
by jfdannen
fix pyOCSM test errors caused by addition of @scope
(detail)
by jfdannen
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
(detail)
by jfdannen
fix stdout sequencing in timPython; get rid of global data in serveESP; update data/python/test* files
(detail)
by jfdannen
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
(detail)
by jfdannen
enable esp.SetModl() to see MODL in serveESP after python finishes
(detail)
by jfdannen
stub routines added to account for change in tim.def
(detail)
by jfdannen
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
(detail)
by jfdannen
add udpTire testing to Makefile.DARWIN64; compute sensitivity correctly with udfCreateBEM, udfCreatePoly, udfEditAttr, udfPrintBbox, udfPrintBrep, and udoPrintEgo; add designS3 test case
(detail)
by jfdannen
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
(detail)
by jfdannen
add udpBodyList to udp.def for Windoze
(detail)
by jfdannen
initial code for proper UDF sensitivity calculations --- not enabled yet
(detail)
by jfdannen
modify udfNacelle* test cases to use DESPMTRs; add splitWing* test cases; remove gallery/FJ2 from distribution
(detail)
by jfdannen
allow pyOCSM to import an externally-generated MODL; enable timMitten; create timPython and data/python/test*; allow serveESP to start with a .py file
(detail)
by jfdannen
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
(detail)
by jfdannen
fix bug in StepThru that was introduced when computation of _tParams was delayed
(detail)
by jfdannen
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
(detail)
by jfdannen
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
(detail)
by jfdannen
add attrScribe test case; add udpParabaloid and test case
(detail)
by jfdannen
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
(detail)
by jfdannen
fix long delay when browser without ball syncs display
(detail)
by jfdannen
remove inside-out test for SolidBodys to improve performance; move Plugs to a TIM in serveESP; do not add _tParams attribute unless the Body will be tessellated; remove coarseBox test case
(detail)
by jfdannen
remove print in tim.c to eliminate Windoze warning
(detail)
by jfdannen
fix bug that sometimes uses FD sensitivities when analytic derivatives could be used instead (associated with UDPARG statements); change method for setting mass properties so that rebuild is not needed when changing spectrum limits in ESP
(detail)
by jfdannen
apply attributes on LINSEG, ARC, CIRARC, BEZIER, SPLINE, and SSLOPE Branches to the Edge that they create; add attrSketch* test cases; remove Gloves and ErepEd from serveCSM; next version of serveESP that uses TIMs; move ErepEd to a TIM in serveESP; move Gloves to a TIM in serveESP; create Mitten as a TIM example
(detail)
by jfdannen
Update EGADS doc for makeTessBody with negative legnth
(detail)
by galbramc
Try fixing bodies in generalBoolean
(detail)
by galbramc
Fix error message. Do not change topObj for bodies. Why does this not work...
(detail)
by galbramc
Correctly set topObj for Geometry in Bodies. Set topObj for bodies in Models. Use EG_removeCntxtRef in makeFace. Check Body references before deleting a Model. More robust cleanup in EG_close.
(detail)
by galbramc
Rule with Loops of differing number of Edges passes simple test with null 'bay'
(detail)
by haimes
Can rule with sections of different number of Edges -- still need to check null bays, multiplicity > 1, curved EndCaps & EG_splineGeom_dot
(detail)
by haimes
Close memory leak in EGADSlite associated with full attributes
(detail)
by haimes
Export Body CSystems that ar free-standing (aLen = 9)
(detail)
by haimes
Update egads_dot.h. Explicit makefile in scan-build/scan-view.
(detail)
by galbramc
change startup messages in wvServer from stderr to stdout
(detail)
by jfdannen
Add missing DEGENERATE Edge check to EG_exportModel
(detail)
by galbramc
Properly reflect the move of OCC copies in this repository
(detail)
by haimes
Fix array out-of-bounds in EG_extrude_dot
(detail)
by galbramc
When splitting periodic Faces convert the order 1 BSpline to a Line; allow for calling wv_cleanupServers when there are active server threads
(detail)
by haimes
Don't fix Faces for EG_makeTopology(SOLIDBODY)
(detail)
by galbramc
C0 Blend fix so that the surface derivatives (normals) do not vanish at the C0
(detail)
by haimes
Scan-build fix from last commit
(detail)
by haimes
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.
(detail)
by haimes
Improved error message. Use brep V2 for OCC 7.6. Minor cleanup.
(detail)
by galbramc
update udpTire.c for change in cacheUdp
(detail)
by jfdannen
Retain all Topology labeled as internal but reject external -- keep Nodes that have a non-zero valence
(detail)
by haimes
Revert Internal Node change -- needs more work!
(detail)
by haimes
An update to the last commit -- remove only Nodes that do not touch active Edges
(detail)
by haimes
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)
(detail)
by haimes
![[Jenkins]](/csi/static/a404d988/images/svgs/logo.svg){kind=logo}
