LCOV - code coverage report
Current view: top level - src - Util.cpp (source / functions) Hit Total Coverage
Test: coverage_sk.info Lines: 1 36 2.8 %
Date: 2020-12-16 07:07:30 Functions: 2 4 50.0 %
Branches: 2 64 3.1 %

           Branch data     Line data    Source code
       1                 :            : /**
       2                 :            :  * MOAB, a Mesh-Oriented datABase, is a software component for creating,
       3                 :            :  * storing and accessing finite element mesh data.
       4                 :            :  *
       5                 :            :  * Copyright 2004 Sandia Corporation.  Under the terms of Contract
       6                 :            :  * DE-AC04-94AL85000 with Sandia Corporation, the U.S. Government
       7                 :            :  * retains certain rights in this software.
       8                 :            :  *
       9                 :            :  * This library is free software; you can redistribute it and/or
      10                 :            :  * modify it under the terms of the GNU Lesser General Public
      11                 :            :  * License as published by the Free Software Foundation; either
      12                 :            :  * version 2.1 of the License, or (at your option) any later version.
      13                 :            :  *
      14                 :            :  */
      15                 :            : 
      16                 :            : //-------------------------------------------------------------------------
      17                 :            : // Purpose       : This file contains utility functions for use in MOAB
      18                 :            : //
      19                 :            : // Special Notes : This is a pure virtual class, to prevent instantiation.
      20                 :            : //                 All functions are static, called like this:
      21                 :            : //                 Util::function_name();
      22                 :            : //-------------------------------------------------------------------------
      23                 :            : #include "moab/Util.hpp"
      24                 :            : #include "moab/Interface.hpp"
      25                 :            : #include <assert.h>
      26                 :            : #include <algorithm>
      27                 :            : #include <limits>
      28                 :            : #if defined( _MSC_VER ) || defined( __MINGW32__ )
      29                 :            : #include <float.h>
      30                 :            : #define finite( A ) _finite( A )
      31                 :            : #ifndef MOAB_HAVE_FINITE
      32                 :            : #define MOAB_HAVE_FINITE
      33                 :            : #endif
      34                 :            : #endif
      35                 :            : 
      36                 :            : namespace moab
      37                 :            : {
      38                 :            : 
      39                 :            : //! temporary normal function for MBEntities.  This should be moved to
      40                 :            : //! an appropriate MB algorithms file
      41                 :            : 
      42                 :          0 : void Util::normal( Interface* MB, EntityHandle handle, double& x, double& y, double& z )
      43                 :            : {
      44                 :            :     // get connectivity
      45                 :          0 :     const EntityHandle* connectivity = NULL;
      46                 :          0 :     int number_nodes                 = 0;
      47                 :            :     // TODO make the return value nonvoid
      48 [ #  # ][ #  # ]:          0 :     ErrorCode rval = MB->get_connectivity( handle, connectivity, number_nodes, true );MB_CHK_SET_ERR_RET( rval, "can't get_connectivity" );
         [ #  # ][ #  # ]
         [ #  # ][ #  # ]
                 [ #  # ]
      49         [ #  # ]:          0 :     assert( number_nodes >= 3 );
      50                 :            : 
      51                 :            :     // get_coordinates
      52                 :            :     double coords[3][3];
      53         [ #  # ]:          0 :     MB->get_coords( &( connectivity[0] ), 1, coords[0] );
      54         [ #  # ]:          0 :     MB->get_coords( &( connectivity[1] ), 1, coords[1] );
      55         [ #  # ]:          0 :     MB->get_coords( &( connectivity[2] ), 1, coords[2] );
      56                 :            : 
      57                 :            :     double vecs[2][3];
      58                 :          0 :     vecs[0][0] = coords[1][0] - coords[0][0];
      59                 :          0 :     vecs[0][1] = coords[1][1] - coords[0][1];
      60                 :          0 :     vecs[0][2] = coords[1][2] - coords[0][2];
      61                 :          0 :     vecs[1][0] = coords[2][0] - coords[0][0];
      62                 :          0 :     vecs[1][1] = coords[2][1] - coords[0][1];
      63                 :          0 :     vecs[1][2] = coords[2][2] - coords[0][2];
      64                 :            : 
      65                 :          0 :     x = vecs[0][1] * vecs[1][2] - vecs[0][2] * vecs[1][1];
      66                 :          0 :     y = vecs[0][2] * vecs[1][0] - vecs[0][0] * vecs[1][2];
      67                 :          0 :     z = vecs[0][0] * vecs[1][1] - vecs[0][1] * vecs[1][0];
      68                 :            : 
      69                 :          0 :     double mag = sqrt( x * x + y * y + z * z );
      70         [ #  # ]:          0 :     if( mag > std::numeric_limits< double >::epsilon() )
      71                 :            :     {
      72                 :          0 :         x /= mag;
      73                 :          0 :         y /= mag;
      74                 :          0 :         z /= mag;
      75                 :            :     }
      76                 :            : }
      77                 :            : 
      78                 :          0 : void Util::centroid( Interface* MB, EntityHandle handle, CartVect& coord )
      79                 :            : {
      80                 :          0 :     const EntityHandle* connectivity = NULL;
      81                 :          0 :     int number_nodes                 = 0;
      82                 :            :     // TODO make the return value nonvoid
      83 [ #  # ][ #  # ]:          0 :     ErrorCode rval = MB->get_connectivity( handle, connectivity, number_nodes, true );MB_CHK_SET_ERR_RET( rval, "can't get_connectivity" );
         [ #  # ][ #  # ]
         [ #  # ][ #  # ]
                 [ #  # ]
      84                 :            : 
      85 [ #  # ][ #  # ]:          0 :     coord[0] = coord[1] = coord[2] = 0.0;
                 [ #  # ]
      86                 :            : 
      87         [ #  # ]:          0 :     for( int i = 0; i < number_nodes; i++ )
      88                 :            :     {
      89                 :            :         double node_coords[3];
      90         [ #  # ]:          0 :         MB->get_coords( &( connectivity[i] ), 1, node_coords );
      91                 :            : 
      92         [ #  # ]:          0 :         coord[0] += node_coords[0];
      93         [ #  # ]:          0 :         coord[1] += node_coords[1];
      94         [ #  # ]:          0 :         coord[2] += node_coords[2];
      95                 :            :     }
      96                 :            : 
      97         [ #  # ]:          0 :     coord[0] /= (double)number_nodes;
      98         [ #  # ]:          0 :     coord[1] /= (double)number_nodes;
      99         [ #  # ]:          0 :     coord[2] /= (double)number_nodes;
     100                 :            : }
     101                 :            : 
     102                 :            : /*//This function calculates the coordinates for the centers of each edges of the entity specified
     103                 :            : by handle. The coordinates are returned in the list coords_list void Util::edge_centers(Interface
     104                 :            : *MB, EntityHandle handle, std::vector<Coord> &coords_list)
     105                 :            : {
     106                 :            :   MB canon_tool(MB);
     107                 :            :   EntityType type;
     108                 :            :   int i = 0;
     109                 :            :   int number_nodes = 0;
     110                 :            :   double coords[2][3];
     111                 :            :   const EntityHandle *connectivity;
     112                 :            : 
     113                 :            :   MB->get_connectivity(handle, connectivity, number_nodes,true);
     114                 :            : 
     115                 :            :   MB->type_from_handle(handle,type);
     116                 :            : 
     117                 :            :   const struct MBCN::ConnMap* conn_map = &(canon_tool.mConnectivityMap[type][0]); //get edge
     118                 :            : sub_elements
     119                 :            : 
     120                 :            :   coords_list.resize(conn_map->num_sub_elements);
     121                 :            : 
     122                 :            :   for(i = 0; i<conn_map->num_sub_elements; i++)
     123                 :            :   {
     124                 :            : 
     125                 :            :     MB->get_coords(connectivity[conn_map->conn[i][0]], coords[0]);
     126                 :            :     MB->get_coords(connectivity[conn_map->conn[i][1]], coords[1]);
     127                 :            : 
     128                 :            :     coords_list[i].x = (coords[0][0] + coords[1][0])/2.0;
     129                 :            :     coords_list[i].y = (coords[0][1] + coords[1][1])/2.0;
     130                 :            :     coords_list[i].z = (coords[0][2] + coords[1][2])/2.0;
     131                 :            :   }
     132                 :            : }
     133                 :            : */
     134                 :            : 
     135                 :            : /*
     136                 :            : void Util::face_centers(Interface *MB, EntityHandle handle, std::vector<Coord> &coords_list)
     137                 :            : {
     138                 :            :   MB canon_tool(MB);
     139                 :            :   EntityType type;
     140                 :            :   int i = 0;
     141                 :            :   int number_nodes = 0;
     142                 :            :   double node_coords[3];
     143                 :            :   const EntityHandle *connectivity;
     144                 :            : 
     145                 :            :   MB->get_connectivity(handle, connectivity, number_nodes,true);
     146                 :            : 
     147                 :            :   MB->type_from_handle(handle,type);
     148                 :            : 
     149                 :            :   const struct MBCN::ConnMap* conn_map = &(canon_tool.mConnectivityMap[type][1]); //get face
     150                 :            : sub_elements
     151                 :            : 
     152                 :            :   coords_list.resize(conn_map->num_sub_elements);
     153                 :            : 
     154                 :            :   for(i = 0; i<conn_map->num_sub_elements;i++)
     155                 :            :   {
     156                 :            :     int number_nodes_per_element = conn_map->num_nodes_per_sub_element[i];
     157                 :            : 
     158                 :            :     for(int j = 0; j<number_nodes_per_element; j++)
     159                 :            :     {
     160                 :            :       MB->get_coords(connectivity[conn_map->conn[i][j]], node_coords);
     161                 :            : 
     162                 :            :       coords_list[i].x+=node_coords[0];
     163                 :            :       coords_list[i].y+=node_coords[1];
     164                 :            :       coords_list[i].z+=node_coords[2];
     165                 :            :     }
     166                 :            : 
     167                 :            :     coords_list[i].x/=(double)number_nodes_per_element;
     168                 :            :     coords_list[i].y/=(double)number_nodes_per_element;
     169                 :            :     coords_list[i].z/=(double)number_nodes_per_element;
     170                 :            :   }
     171                 :            : }
     172                 :            : */
     173                 :            : 
     174                 :            : // Explicit template specializations
     175                 :            : template bool Util::is_finite< double >( double value );
     176                 :            : 
     177 [ +  - ][ +  - ]:        228 : }  // namespace moab

Generated by: LCOV version 1.11