<html>
<head>
<meta http-equiv="Content-Type" content="text/html; charset=us-ascii">
</head>
<body style="word-wrap: break-word; -webkit-nbsp-mode: space; -webkit-line-break: after-white-space; color: rgb(0, 0, 0); font-size: 14px; font-family: Calibri, sans-serif;">
<div>
<div>
<div>The implementation of the compute derivatives filter has been around a while so the author probably hasn't seen your email. I'm not super familiar with the code, but I believe the vorticity is computed from the first order derivative of the vector field.</div>
<div><br>
</div>
<div>It is not actually using a central difference scheme to compute the gradients. Instead, it is treating the field as a piecewise linear field (or piecewise quadratic or piecewise cubic depending on the cell type) based on the volumes of the cells. The derivatives
 are then computed in the middle of the piecewise C1 continuous parts of the cells from the known function.</div>
<div><br>
</div>
<div>You can take a look at the source for the algorithm here:</div>
<div><br>
</div>
<div><a href="http://vtk.org/gitweb?p=VTK.git;a=blob;f=Filters/General/vtkCellDerivatives.cxx;h=4debdf45eeb7a11e98062ab2d9efa65e938433b6;hb=ab4dc579c340167b6b3ae6cd82b3fd645642796a">http://vtk.org/gitweb?p=VTK.git;a=blob;f=Filters/General/vtkCellDerivatives.cxx;h=4debdf45eeb7a11e98062ab2d9efa65e938433b6;hb=ab4dc579c340167b6b3ae6cd82b3fd645642796a</a></div>
</div>
</div>
<div><br>
</div>
<div>Vorticity is computed between lines 169 and 216.</div>
<div><br>
</div>
<div>-Ken</div>
<div><br>
</div>
<span id="OLK_SRC_BODY_SECTION">
<div style="font-family:Calibri; font-size:11pt; text-align:left; color:black; BORDER-BOTTOM: medium none; BORDER-LEFT: medium none; PADDING-BOTTOM: 0in; PADDING-LEFT: 0in; PADDING-RIGHT: 0in; BORDER-TOP: #b5c4df 1pt solid; BORDER-RIGHT: medium none; PADDING-TOP: 3pt">
<span style="font-weight:bold">From: </span>Jean Hertzberg <<a href="mailto:Hertzberg@Colorado.EDU">Hertzberg@Colorado.EDU</a>><br>
<span style="font-weight:bold">Date: </span>Wednesday, February 11, 2015 at 12:01 PM<br>
<span style="font-weight:bold">To: </span>"<a href="mailto:paraview@paraview.org">paraview@paraview.org</a>" <<a href="mailto:paraview@paraview.org">paraview@paraview.org</a>><br>
<span style="font-weight:bold">Subject: </span>[EXTERNAL] [Paraview] Vorticity calc algorithm?<br>
</div>
<div><br>
</div>
<div xmlns:v="urn:schemas-microsoft-com:vml" xmlns:o="urn:schemas-microsoft-com:office:office" xmlns:w="urn:schemas-microsoft-com:office:word" xmlns:m="http://schemas.microsoft.com/office/2004/12/omml" xmlns="http://www.w3.org/TR/REC-html40">
<meta name="Generator" content="Microsoft Word 15 (filtered medium)">
<style><!--
/* Font Definitions */
@font-face
        {font-family:"Cambria Math";
        panose-1:2 4 5 3 5 4 6 3 2 4;}
@font-face
        {font-family:Calibri;
        panose-1:2 15 5 2 2 2 4 3 2 4;}
/* Style Definitions */
p.MsoNormal, li.MsoNormal, div.MsoNormal
        {margin:0in;
        margin-bottom:.0001pt;
        font-size:11.0pt;
        font-family:"Calibri",sans-serif;}
a:link, span.MsoHyperlink
        {mso-style-priority:99;
        color:#0563C1;
        text-decoration:underline;}
a:visited, span.MsoHyperlinkFollowed
        {mso-style-priority:99;
        color:#954F72;
        text-decoration:underline;}
span.EmailStyle17
        {mso-style-type:personal-compose;
        font-family:"Calibri",sans-serif;
        color:windowtext;}
.MsoChpDefault
        {mso-style-type:export-only;
        font-family:"Calibri",sans-serif;}
@page WordSection1
        {size:8.5in 11.0in;
        margin:1.0in 1.0in 1.0in 1.0in;}
div.WordSection1
        {page:WordSection1;}
--></style><!--[if gte mso 9]><xml>
<o:shapedefaults v:ext="edit" spidmax="1026" />
</xml><![endif]--><!--[if gte mso 9]><xml>
<o:shapelayout v:ext="edit">
<o:idmap v:ext="edit" data="1" />
</o:shapelayout></xml><![endif]-->
<div lang="EN-US" link="#0563C1" vlink="#954F72">
<div class="WordSection1">
<p class="MsoNormal">Hi <o:p></o:p></p>
<p class="MsoNormal">I asked this last week, but got zero response. Please help me, or point me towards the right documentation:<o:p></o:p></p>
<p class="MsoNormal">Exactly what algorithm is used to calculate vorticity in the Compute Derivatives filter? Is it a first or a second order central difference scheme for gradients, or does it use the circulation method?<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Thanks<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Jean Hertzberg<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
<p class="MsoNormal">Dept. Mechanical Engineering<o:p></o:p></p>
<p class="MsoNormal">University of Colorado, Boulder<o:p></o:p></p>
<p class="MsoNormal"><o:p> </o:p></p>
</div>
</div>
</div>
</span>
</body>
</html>