% get K subsets of gradients from one given gradient set, such that
% gradients in each subset are evenly distributed.
%
% USAGE:
% [gradCell,grb_result,indexMatrix] = OptimalSamplingMultiSubsetsFromSameSet(gradAll, param, grbParam)
%
% INPUT
% gradAll : Nx3 gradient matrix, where each row is a point in sphere.
% param.numSamples : Kx1 or 1xK vector which gives number of samples in K subsets
% param.w : weight for single shell term. 0<w<1. Default value is 0.5.
% param.sos : 1 (default), use sos constraint
% 0, do not use sos
% param.ModelFile : string. If set, write the model into a file
%
% grbParam.start : a warm start if provided
% grbParam.TimeLimit : time limit to terminate the program
% grbParam.MIPFocus : MIPFocus
%
% OUTPUT
% gradCell : Kx1 cell, each element is a gradient matrix.
% grb_result : output from GUROBI
%
%
% Reference:
% 1. "Single- and Multiple-Shell Uniform Sampling Schemes for Diffusion MRI Using Spherical Codes",
% Jian Cheng, Dinggang Shen, Pew-Thian Yap, Peter J. Basser, IEEE Transactions on Medical Imaging, 2017.
% 2. "Designing Single- and Multiple-Shell Sampling Schemes for Diffusion MRI Using Spherical Code",
% Jian Cheng, Dinggang Shen, Pew-Thian Yap, MICCAI 2014.
%
% Copyright (c) 2013, Jian Cheng <jian.cheng.1983@gmail.com>
%
