active_rgn.m

software to find the threshold of an LDPC code in a Gaussian channel

function [A_r,map_r,K_r] = active_rgn(fp,A,delta,lpad,K_rmax)

foo = sum(cumsum(fp) < (sum(fp)-1e-6));
K_r = foo;

map_r = 1:K_r;

A_r = A(1:(foo+2*lpad+1),map_r);