diff --git a/GENetLib/__init__.py b/GENetLib/__init__.py
index c545a60..602a4a5 100644
--- a/GENetLib/__init__.py
+++ b/GENetLib/__init__.py
@@ -2,11 +2,9 @@
 from .create_basis import create_bspline_basis,create_expon_basis,create_fourier_basis,create_monomial_basis,create_power_basis,create_constant_basis
 from .basis_mat import bspline_mat,expon_mat,fourier_mat,monomial_mat,polyg_mat,power_mat
 from .get_basis_matrix import get_basis_matrix
-from .pp_func import ppbspline,ppderiv
 from .GE_Net import GE_Net
 from .survival_costfunc_cindex import neg_par_log_likelihood, c_index
 from .inprod import inprod
-from .inprod_bspline import inprod_bspline
 from .fd import fd
 from .pre_data1 import pre_data1
 from .pre_data2 import pre_data2
@@ -16,7 +14,6 @@
 from .fd_chk import fd_chk
 from .knotmultchk import knotmultchk
 from .eval_basis_fd import eval_basis,eval_fd
-from .polyprod import polyprod
 from .dense_to_func import dense_to_func
 from .scalar_l2train import scalar_l2train
 from .scalar_mcp_l2train import scalar_mcp_l2train
@@ -32,8 +29,8 @@
 __all__ = ['basis_fd', 'create_bspline_basis','create_expon_basis','create_fourier_basis',
            'create_monomial_basis','create_power_basis','create_constant_basis',
            'bspline_mat','expon_mat','fourier_mat','monomial_mat','polyg_mat','power_mat',
-           'get_basis_matrix','ppbspline','ppderiv','GE_Net','neg_par_log_likelihood','c_index',
-           'inprod','inprod_bspline','fd','pre_data1','pre_data2','sim_data_scalar','sim_data_func',
-           'spline_design','fd_chk','knotmultchk','eval_basis','eval_fd','polyprod','dense_to_func',
+           'get_basis_matrix','GE_Net','neg_par_log_likelihood','c_index',
+           'inprod','fd','pre_data1','pre_data2','sim_data_scalar','sim_data_func',
+           'spline_design','fd_chk','knotmultchk','eval_basis','eval_fd','dense_to_func',
            'scalar_l2train','scalar_mcp_l2train','scalar_ge','grid_scalar_ge','func_ge','grid_func_ge',
            'plot_fd','plot_rawdata','plot_func']