Commit 29bd6af8 authored by Daniel Brown's avatar Daniel Brown

Fixing missing refractive index factor

parent c49ef0ac
......@@ -1577,7 +1577,7 @@ void fill_mirror_signal_rhs(signal_t *signal, ifo_matrix_vars_t* M_car, ifo_matr
if (!node1->gnd_node) {
// Get output at node 1
sidx = M_sig->node_rhs_idx_2[node1->list_index] + get_port_rhs_idx(f, field_index);
srhs[sidx] = z_by_xph(a11, turn * factor_r, phase_r);
srhs[sidx] = z_by_xph(a11, turn * factor_r * nr1, phase_r);
if (!node2->gnd_node) {
z_inc_z(&srhs[sidx], z_by_xph(a21, (nr1-nr2)*factor_t, phase_t));
......@@ -1589,7 +1589,7 @@ void fill_mirror_signal_rhs(signal_t *signal, ifo_matrix_vars_t* M_car, ifo_matr
if (!node2->gnd_node) {
sidx = M_sig->node_rhs_idx_1[node2->list_index] + get_port_rhs_idx(f, field_index);
srhs[sidx] = z_by_xph(a22, turn * factor_r, phase_r + ph2);
srhs[sidx] = z_by_xph(a22, turn * factor_r * nr2, phase_r + ph2);
if (!node1->gnd_node) {
z_inc_z(&srhs[sidx], z_by_xph(a12, -(nr1-nr2)*factor_t, phase_t));
......
Markdown is supported
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