testgC                 package:ccems                 R Documentation

_T_e_s_t _g _O_b_j_e_c_t _T_C_C _C _C_o_m_p_i_l_a_t_i_o_n

_D_e_s_c_r_i_p_t_i_o_n:

     This function runs 'lsoda' using the compiled C code of the total
     concentration constraint (TCC)  system equations that should have
     been created by 'mkg'.

_U_s_a_g_e:

      testgC(g) 

_A_r_g_u_m_e_n_t_s:

       g: The output of 'mkg'. Since 'mkg' already calls this function,
          it has no independent use. 

_D_e_t_a_i_l_s:

     Components of 'g' used by this function are: 'g$wDir' which is the
     location of the 'models' directory that holds the compiled code; 
     'g$parmsTCC' which includes the total concentrations, the fraction
     of active protein p, and the spur graph complete dissociation
     constant parameter values Kj; and the 'lsoda' ODE integration
     tolerance parameters 'g$rtol' and 'g$atol'.

_V_a_l_u_e:

     No value returned. Ten time points are printed to the screen. The
     free concentrations should show signs of convergence toward a
     steady state.

_A_u_t_h_o_r(_s):

     Tom Radivoyevitch

_S_e_e _A_l_s_o:

     'lsoda'

_E_x_a_m_p_l_e_s:

     library(ccems)
     topology <- list(  
             heads=c("R1t0","R2t0"),  
             sites=list(       
                     s=list(                     # s-site    thread #
                             m=c("R1t1"),        # monomer      1
                             d=c("R2t1","R2t2")  # dimer        2
                     )
             )
     ) 
     g=mkg(topology,TCC=TRUE) 
     testgC(g) # redundantly, this checks it again  ... there should never be a reason to call this function.

