The C++ function grad_HMM_stat_dir_re() updates the gradient function when the number of cluster is 1 in a static binary directed network. (The above description may need further improvement.)
grad_HMM_stat_dir_re_K1(theta, network, N)
Vector of current theta
(network canonical parameters) values.
Array of a network adjacency matrices.
Number of nodes.
Returns the gradient vector of ELBO with respect to theta parameters for Newton-Raphson update in the M-step.