macarray_interpolator3.h

Go to the documentation of this file.

/*
**  macarray_interpolator3.h
**
**  This is part of Shiokaze, a research-oriented fluid solver for computer graphics.
**  Created by Ryoichi Ando <rand@nii.ac.jp> on Feb 25, 2018.
**
**  Permission is hereby granted, free of charge, to any person obtaining a copy of
**  this software and associated documentation files (the "Software"), to deal in
**  the Software without restriction, including without limitation the rights to use,
**  copy, modify, merge, publish, distribute, sublicense, and/or sell copies of the
**  Software, and to permit persons to whom the Software is furnished to do so,
**  subject to the following conditions:
**
**  The above copyright notice and this permission notice shall be included in all copies
**  or substantial portions of the Software.
**
**  THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS OR IMPLIED,
**  INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY, FITNESS FOR A
**  PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL THE AUTHORS OR COPYRIGHT
**  HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER LIABILITY, WHETHER IN AN ACTION OF
**  CONTRACT, TORT OR OTHERWISE, ARISING FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE
**  OR THE USE OR OTHER DEALINGS IN THE SOFTWARE.
*/
//
#ifndef SHKZ_MACARRAY_INTERPOLATOR3_H
#define SHKZ_MACARRAY_INTERPOLATOR3_H
//
#include <shiokaze/math/vec.h>
#include "array_interpolator3.h"
//
SHKZ_BEGIN_NAMESPACE
//
namespace macarray_interpolator3 {
    template<class T> static vec3<T> interpolate( const macarray3<T> &array, const vec3d &p, bool only_actives=false ) {
        vec3<T> result;
        for( int dim : DIMS3 ) {
            const vec3d pos = vec3d(p[0]-0.5*(dim!=0),p[1]-0.5*(dim!=1),p[2]-0.5*(dim!=2));
            result[dim] = array_interpolator3::interpolate<T>(array[dim],pos,only_actives);
        }
        return result;
    }
    template<class T> static vec3<T> interpolate( const macarray3<T> &array, const vec3d &origin, double dx, const vec3d &p, bool only_actives=false ) {
        return interpolate<T>(array,(p-origin)/dx,only_actives);
    }
};
//
SHKZ_END_NAMESPACE
//
#endif
//