splines.h

Go to the documentation of this file.

/***************************************************************************
 *
 * Authors:     Carlos Oscar S. Sorzano (coss@cnb.csic.es)
 *
 * Unidad de  Bioinformatica of Centro Nacional de Biotecnologia , CSIC
 *
 * This program is free software; you can redistribute it and/or modify
 * it under the terms of the GNU General Public License as published by
 * the Free Software Foundation; either version 2 of the License, or
 * (at your option) any later version.
 *
 * This program is distributed in the hope that it will be useful,
 * but WITHOUT ANY WARRANTY; without even the implied warranty of
 * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE.  See the
 * GNU General Public License for more details.
 *
 * You should have received a copy of the GNU General Public License
 * along with this program; if not, write to the Free Software
 * Foundation, Inc., 59 Temple Place, Suite 330, Boston, MA
 * 02111-1307  USA
 *
 *  All comments concerning this program package may be sent to the
 *  e-mail address 'xmipp@cnb.csic.es'
 ***************************************************************************/

#ifndef _SPLINES_HH
#define _SPLINES_HH

#include <core/bilib/kernel.h>

#include "grids.h"


/* ========================================================================= */
/* SPLINES                                                                   */
/* ========================================================================= */
// Spline value ------------------------------------------------------------
inline double spatial_Bspline03(const Matrix1D<double> &r)
{
    if (-2 <= XX(r) && XX(r) < 2 &&
        -2 <= YY(r) && YY(r) < 2 &&
        -2 <= ZZ(r) && ZZ(r) < 2)
        return Bspline03(XX(r))*Bspline03(YY(r))*Bspline03(ZZ(r));
    else return 0.0;
}

const int BSPLINE03_SUBSAMPLING = 2000;
double Bspline03LUT(double x);

inline double spatial_Bspline03LUT(const Matrix1D<double> &r)
{
    if (-2 <= XX(r) && XX(r) < 2 &&
        -2 <= YY(r) && YY(r) < 2 &&
        -2 <= ZZ(r) && ZZ(r) < 2)
        return Bspline03LUT(XX(r))*Bspline03LUT(YY(r))*Bspline03LUT(ZZ(r));
    else return 0.0;
}

double sum_spatial_Bspline03_Grid(const Grid &grid);

// Spline projection -------------------------------------------------------
double spatial_Bspline03_proj(const Matrix1D<double> &r, const Matrix1D<double> &u);

void spatial_Bspline032voxels(const GridVolume &vol_splines,
                              MultidimArray<double> *vol_voxels,
                              int Zdim = 0, int Ydim = 0, int Xdim = 0);
#endif