gamini.m

function bigger=gamini(data,folding)
% bigger=GAMINI(data,folding)
%
% Replicates the ith entry of a vector as many times as the ith element of
% another vector (or a fixed number of times when given only one value ).
%
% INPUT:
%
% data       A certain data vector
% folding    The folding factor(s)
%            if the same size as the data vector, element-wise replication;
%            if a scalar, it applies to all data elements [default: 3];
%            if zero or negative, no replication occurs
%
% EXAMPLE:
%
% [a,b]=degamini(gamini([1 2 3 1 4 5 ],[1 2 3 2 4 2]))
%
% See also DEGAMINI,GAMINI2
%
% Last modified by fjsimons-at-alum.mit.edu, 05/26/2021

defval('folding',3)

% Copy a single scalar folding for all elements in the array 
if prod(size(folding))==1
  folding=repmat(folding,prod(size(data)),1);
end

% Never had a replication factor of 0 before
% But only if they are of the same length!
if length(data)==length(folding)
  data=data(folding>0);
  folding=folding(folding>0);
end

% Rearrange
data=data(:)';
folding=folding(:)';

if ~all(size(data)==size(folding)) 
  error([ ' Sizes of input and folding must be the same'])
end

gelp=zeros(1,sum(folding));
gelp([1 cumsum(folding(1:end-1))+1])=1; 
if ~isempty(data)
  bigger=data(cumsum(gelp));
else
  bigger=[];
end