%--------------------------
 % @Author: Jingqiao Hu
 % @Date: 2021-12-23 12:27:07
 % @LastEditTime: 2021-12-23 12:27:10
%--------------------------
function [xval, xold1, xold2, xmin, xmax, low, upp, m_mma, c_mma, d_mma, a0_mma, a_mma] = MMA_paras_mmc(t, t0, t1, nc_local, nd)

    nc_global = sum(nc_local) * nd;
    
    xval = repmat(t, nc_global, 1);
    xmin = repmat(t0,nc_global, 1);
    xmax = repmat(t1,nc_global, 1);

    xold1= xval; xold2 = xval; 
    low  = xmin; upp = xmax; 

    m_mma = 1; % number of constrainm 
    c_mma = 1000*ones(m_mma,1); d_mma = zeros(m_mma,1); a0_mma = 1; a_mma = zeros(m_mma,1); 
end