Commit 6177c818 authored by Renaud Pacalet's avatar Renaud Pacalet
Browse files

Add alpha and beta initialization for scmsfp3

parent d3511190
......@@ -1000,6 +1000,13 @@ void scmsfp3_init(mod2sparse *H, double *lratio, char *dblk, double *bprb, scmsf
sd->g[col] = scms_saturate(round(dtmp), sd->size);
sd->h[col] = sd->g[col];
}
for(row = 0; row < M; row++) {
for (e = mod2sparse_first_in_row(H,row), i = 0; !mod2sparse_at_end(e); e = mod2sparse_next_in_row(e), i++) { // for all active columns in row
col = e->col; // column index
sd->b[row][i] = 0; // CN
sd->a[row][i] = 0; // VN
}
}
}
/* compute one iteration */
......
Supports Markdown
0% or .
You are about to add 0 people to the discussion. Proceed with caution.
Finish editing this message first!
Please register or to comment