Oscillatory Shear Index of a Point Cloud

A function to calculate the Oscillatory Shear Index of a Point Cloud

This is a function that calculates the oscillatory shear index (OSI) of a point cloud. It is based on definition of OSI given by Soulis et al. 2011 - OSI = 0.5*(1.0-AWSSV/AWSS), where AWSSV = magnitude of the time-averaged WSS vector, and AWSS = time-averaged WSS magnitude. The code is intended for use as Fluent/CFD-post data post processing. As such, inputs are formatted to coincide with the format that is output by Fluent with each time-step.
Authors: G. Riccardello, K. Thomas, A. Changa.
INPUT
pointcld_WSS_WSScomp = nx7xs array [Xcoordinate Ycoordinate Zcoordinate...
MagnitudeWSS XcomponentWSS YcomponentWSS ZcomponentWSS]
where n = number of points in the point cloud and s = number
of steps in the simulation.

OUTPUT
pointcld_OSI = nx4 array [Xcoordinate Ycoordinate Zcoordinate OSI]
where n = number of points in the point cloud

Soulis, Johannes V., Olga P. Lampri, Dimitrios K. Fytanidis, and George D. Giannoglou. "Relative Residence Time and Oscillatory Shear Index of Non-Newtonian Flow Models in Aorta." 2011 10th International Workshop on Biomedical Engineering (2011): n. pag. Web.