LCOV - code coverage report
Current view: top level - src/verdict - V_KnifeMetric.cpp (source / functions) Hit Total Coverage
Test: coverage_sk.info Lines: 0 37 0.0 %
Date: 2020-12-16 07:07:30 Functions: 0 2 0.0 %
Branches: 0 50 0.0 %

           Branch data     Line data    Source code
       1                 :            : /*=========================================================================
       2                 :            : 
       3                 :            :   Module:    $RCSfile: V_KnifeMetric.cpp,v $
       4                 :            : 
       5                 :            :   Copyright (c) 2006 Sandia Corporation.
       6                 :            :   All rights reserved.
       7                 :            :   See Copyright.txt or http://www.kitware.com/Copyright.htm for details.
       8                 :            : 
       9                 :            :      This software is distributed WITHOUT ANY WARRANTY; without even
      10                 :            :      the implied warranty of MERCHANTABILITY or FITNESS FOR A PARTICULAR
      11                 :            :      PURPOSE.  See the above copyright notice for more information.
      12                 :            : 
      13                 :            : =========================================================================*/
      14                 :            : 
      15                 :            : /*
      16                 :            :  *
      17                 :            :  * KnifeMetrics.cpp contains quality calculations for knives
      18                 :            :  *
      19                 :            :  * This file is part of VERDICT
      20                 :            :  *
      21                 :            :  */
      22                 :            : 
      23                 :            : #define VERDICT_EXPORTS
      24                 :            : 
      25                 :            : #include "moab/verdict.h"
      26                 :            : #include "VerdictVector.hpp"
      27                 :            : #include <memory.h>
      28                 :            : 
      29                 :            : /*  a knife element
      30                 :            : 
      31                 :            :           3
      32                 :            :          _/\_
      33                 :            :        _/  | \_
      34                 :            :    0 _/        \_ 2
      35                 :            :     |\_    | ___/|
      36                 :            :     |  \  __/    |
      37                 :            :     |  1\/ |     |
      38                 :            :     |    \       |
      39                 :            :     |_____\|_____|
      40                 :            :    4       5      6
      41                 :            : 
      42                 :            : 
      43                 :            :     (edge 3,5  is is a hidden line if you will)
      44                 :            : 
      45                 :            :     if this is hard to visualize, consider a hex
      46                 :            :     with nodes 5 and 7 becoming the same node
      47                 :            : 
      48                 :            : 
      49                 :            : */
      50                 :            : 
      51                 :            : /*!
      52                 :            :   calculates the volume of a knife element
      53                 :            : 
      54                 :            :   this is done by dividing the knife into 4 tets
      55                 :            :   and summing the volumes of each.
      56                 :            : */
      57                 :            : 
      58                 :          0 : C_FUNC_DEF double v_knife_volume( int num_nodes, double coordinates[][3] )
      59                 :            : {
      60                 :          0 :     double volume = 0;
      61 [ #  # ][ #  # ]:          0 :     VerdictVector side1, side2, side3;
                 [ #  # ]
      62                 :            : 
      63         [ #  # ]:          0 :     if( num_nodes == 7 )
      64                 :            :     {
      65                 :            : 
      66                 :            :         // divide the knife into 4 tets and calculate the volume
      67                 :            : 
      68                 :          0 :         side1.set( coordinates[1][0] - coordinates[0][0], coordinates[1][1] - coordinates[0][1],
      69         [ #  # ]:          0 :                    coordinates[1][2] - coordinates[0][2] );
      70                 :          0 :         side2.set( coordinates[3][0] - coordinates[0][0], coordinates[3][1] - coordinates[0][1],
      71         [ #  # ]:          0 :                    coordinates[3][2] - coordinates[0][2] );
      72                 :          0 :         side3.set( coordinates[4][0] - coordinates[0][0], coordinates[4][1] - coordinates[0][1],
      73         [ #  # ]:          0 :                    coordinates[4][2] - coordinates[0][2] );
      74                 :            : 
      75 [ #  # ][ #  # ]:          0 :         volume = side3 % ( side1 * side2 ) / 6;
      76                 :            : 
      77                 :          0 :         side1.set( coordinates[5][0] - coordinates[1][0], coordinates[5][1] - coordinates[1][1],
      78         [ #  # ]:          0 :                    coordinates[5][2] - coordinates[1][2] );
      79                 :          0 :         side2.set( coordinates[3][0] - coordinates[1][0], coordinates[3][1] - coordinates[1][1],
      80         [ #  # ]:          0 :                    coordinates[3][2] - coordinates[1][2] );
      81                 :          0 :         side3.set( coordinates[4][0] - coordinates[1][0], coordinates[4][1] - coordinates[1][1],
      82         [ #  # ]:          0 :                    coordinates[4][2] - coordinates[1][2] );
      83                 :            : 
      84 [ #  # ][ #  # ]:          0 :         volume += side3 % ( side1 * side2 ) / 6;
      85                 :            : 
      86                 :          0 :         side1.set( coordinates[2][0] - coordinates[1][0], coordinates[2][1] - coordinates[1][1],
      87         [ #  # ]:          0 :                    coordinates[2][2] - coordinates[1][2] );
      88                 :          0 :         side2.set( coordinates[3][0] - coordinates[1][0], coordinates[3][1] - coordinates[1][1],
      89         [ #  # ]:          0 :                    coordinates[3][2] - coordinates[1][2] );
      90                 :          0 :         side3.set( coordinates[6][0] - coordinates[1][0], coordinates[6][1] - coordinates[1][1],
      91         [ #  # ]:          0 :                    coordinates[6][2] - coordinates[1][2] );
      92                 :            : 
      93 [ #  # ][ #  # ]:          0 :         volume += side3 % ( side1 * side2 ) / 6;
      94                 :            : 
      95                 :          0 :         side1.set( coordinates[3][0] - coordinates[1][0], coordinates[3][1] - coordinates[1][1],
      96         [ #  # ]:          0 :                    coordinates[3][2] - coordinates[1][2] );
      97                 :          0 :         side2.set( coordinates[5][0] - coordinates[1][0], coordinates[5][1] - coordinates[1][1],
      98         [ #  # ]:          0 :                    coordinates[5][2] - coordinates[1][2] );
      99                 :          0 :         side3.set( coordinates[6][0] - coordinates[1][0], coordinates[6][1] - coordinates[1][1],
     100         [ #  # ]:          0 :                    coordinates[6][2] - coordinates[1][2] );
     101                 :            : 
     102 [ #  # ][ #  # ]:          0 :         volume += side3 % ( side1 * side2 ) / 6;
     103                 :            :     }
     104                 :            : 
     105                 :          0 :     return (double)volume;
     106                 :            : }
     107                 :            : 
     108                 :            : /*!
     109                 :            : 
     110                 :            :   calculate the quality metrics of a knife element.
     111                 :            : 
     112                 :            :   There is only one, but we put this here to be consistent with
     113                 :            :   functions for other element types.  Who knows if we'll add
     114                 :            :   more metrics.
     115                 :            : */
     116                 :            : 
     117                 :          0 : C_FUNC_DEF void v_knife_quality( int num_nodes, double coordinates[][3], unsigned int metrics_request_flag,
     118                 :            :                                  KnifeMetricVals* metric_vals )
     119                 :            : {
     120                 :          0 :     memset( metric_vals, 0, sizeof( KnifeMetricVals ) );
     121                 :            : 
     122         [ #  # ]:          0 :     if( metrics_request_flag & V_KNIFE_VOLUME ) metric_vals->volume = v_knife_volume( num_nodes, coordinates );
     123                 :          0 : }

Generated by: LCOV version 1.11