setupequations.h

Code for L2-optimal three view triangulation based on calculation of stationary points

  eqy[201] = ZERO;
  eqy[202] = ZERO;
  eqy[203] = ZERO;
  eqy[204] = ZERO;
  eqy[205] = ZERO;
  eqy[206] = ZERO;
  eqy[207] = ZERO;
  eqy[208] = ZERO;

  //EQZ
  eqz[0] = ZERO;
  eqz[1] = ZERO;
  eqz[2] = ZERO;
  eqz[3] = ZERO;
  eqz[4] = ZERO;
  eqz[5] = ZERO;
  eqz[6] = ZERO;
  eqz[7] = ZERO;
  eqz[8] = ZERO;
  eqz[9] = ZERO;
  eqz[10] = ZERO;
  eqz[11] = ZERO;
  eqz[12] = ZERO;
  eqz[13] = ZERO;
  eqz[14] = ZERO;
  eqz[15] = ZERO;
  eqz[16] = ZERO;
  eqz[17] = ZERO;
  eqz[18] = ZERO;
  eqz[19] = ZERO;
  eqz[20] = ZERO;
  eqz[21] = ZERO;
  eqz[22] = ZERO;
  eqz[23] = ZERO;
  eqz[24] = ZERO;
  eqz[25] = ZERO;
  eqz[26] = ZERO;
  eqz[27] = ZERO;
  eqz[28] = ZERO;
  eqz[29] = ZERO;
  eqz[30] = ZERO;
  eqz[31] = ZERO;
  eqz[32] = ZERO;
  eqz[33] = ZERO;
  eqz[34] = ZERO;
  eqz[35] = ZERO;
  eqz[36] = ZERO;
  eqz[37] = ZERO;
  eqz[38] = ZERO;
  eqz[39] = ZERO;
  eqz[40] = ZERO;
  eqz[41] = ZERO;
  eqz[42] = ZERO;
  eqz[43] = ZERO;
  eqz[44] = ZERO;
  eqz[45] = ZERO;
  eqz[46] = ZERO;
  eqz[47] = ZERO;
  eqz[48] = ZERO;
  eqz[49] = ZERO;
  eqz[50] = ZERO;
  eqz[51] = ZERO;
  eqz[52] = ZERO;
  eqz[53] = ZERO;
  eqz[54] = ZERO;
  eqz[55] = ZERO;
  eqz[56] = ZERO;
  eqz[57] = ZERO;
  eqz[58] = ZERO;
  eqz[59] = ZERO;
  eqz[60] = ZERO;
  eqz[61] = ZERO;
  eqz[62] = ZERO;
  eqz[63] = ZERO;
  eqz[64] = ZERO;
  eqz[65] = ZERO;
  eqz[66] = ZERO;
  eqz[67] = ZERO;
  eqz[68] = ZERO;
  eqz[69] = ZERO;
  eqz[70] = ZERO;
  eqz[71] = ZERO;
  eqz[72] = ZERO;
  eqz[73] = ZERO;
  eqz[74] = ZERO;
  eqz[75] = ZERO;
  eqz[76] = ZERO;
  eqz[77] = ZERO;
  eqz[78] = ZERO;
  eqz[79] = ZERO;
  eqz[80] = ZERO;
  eqz[81] = ZERO;
  eqz[82] = ZERO;
  eqz[83] = ZERO;
  eqz[84] = ZERO;
  eqz[85] = ZERO;
  eqz[86] = ZERO;
  eqz[87] = ZERO;
  eqz[88] = ZERO;
  eqz[89] = ZERO;
  eqz[90] = ZERO;
  eqz[91] = ZERO;
  eqz[92] = ZERO;
  eqz[93] = ZERO;
  eqz[94] = ZERO;
  eqz[95] = ZERO;
  eqz[96] = ZERO;
  eqz[97] = ZERO;
  eqz[98] = ZERO;
  eqz[99] = ZERO;
  eqz[100] = ZERO;
  eqz[101] = ZERO;
  eqz[102] = ZERO;
  eqz[103] = ZERO;
  eqz[104] = ZERO;
  eqz[105] = ZERO;
  eqz[106] = ZERO;
  eqz[107] = ZERO;
  eqz[108] = ZERO;
  eqz[109] = ZERO;
  eqz[110] = ZERO;
  eqz[111] = ZERO;
  eqz[112] = ZERO;
  eqz[113] = ZERO;
  eqz[114] = ZERO;
  eqz[115] = ZERO;
  eqz[116] = ZERO;
  eqz[117] = ZERO;
  eqz[118] = ZERO;
  eqz[119] = ZERO;
  eqz[120] = ZERO;
  eqz[121] = ZERO;
  eqz[122] = ZERO;
  eqz[123] = ZERO;
  eqz[124] = ZERO;
  eqz[125] = ZERO;
  eqz[126] = ZERO;
  eqz[127] = ZERO;
  eqz[128] = ZERO;
  eqz[129] = ZERO;
  eqz[130] = ZERO;
  eqz[131] = ZERO;
  eqz[132] = ZERO;
  eqz[133] = ZERO;
  eqz[134] = ZERO;
  eqz[135] = ZERO;
  eqz[136] = ZERO;
  eqz[137] = ZERO;
  eqz[138] = ZERO;
  eqz[139] = -0.2e1 * pow(P3[1][0], 0.2e1) - 0.2e1 * pow(P3[0][0], 0.2e1);
  eqz[140] = ZERO;
  eqz[141] = ZERO;
  eqz[142] = -0.4e1 * P3[1][1] * P3[1][0] - 0.4e1 * P3[0][1] * P3[0][0];
  eqz[143] = -0.2e1 * P3[1][2] * P3[1][0] - 0.2e1 * P3[0][2] * P3[0][0];
  eqz[144] = ZERO;
  eqz[145] = 2 * P2[0][2] * P2[0][0] + 2 * P2[1][2] * P2[1][0];
  eqz[146] = -0.2e1 * pow(P3[1][1], 0.2e1) - 0.2e1 * pow(P3[0][1], 0.2e1);
  eqz[147] = -0.2e1 * P3[0][1] * P3[0][2] - 0.2e1 * P3[1][1] * P3[1][2];
  eqz[148] = ZERO;
  eqz[149] = 2 * P2[0][1] * P2[0][2] + 2 * P2[1][1] * P2[1][2];
  eqz[150] = 2 * pow((double) P2[0][2], (double) 2) + 2 * pow((double) P2[1][2], (double) 2);
  eqz[151] = ZERO;
  eqz[152] = ZERO;
  eqz[153] = ZERO;
  eqz[154] = 2 * P1[0][2] * P1[0][0] + 2 * P1[1][2] * P1[1][0];
  eqz[155] = ZERO;
  eqz[156] = ZERO;
  eqz[157] = ZERO;
  eqz[158] = ZERO;
  eqz[159] = ZERO;
  eqz[160] = 2 * P1[0][1] * P1[0][2] + 2 * P1[1][1] * P1[1][2];
  eqz[161] = 2 * pow((double) P1[1][2], (double) 2) + 2 * pow((double) P1[0][2], (double) 2);
  eqz[162] = ZERO;
  eqz[163] = ZERO;
  eqz[164] = ZERO;
  eqz[165] = ZERO;
  eqz[166] = ZERO;
  eqz[167] = -0.4e1 * P3[1][0] * P3[1][3] - 0.4e1 * P3[0][0] * P3[0][3];
  eqz[168] = ZERO;
  eqz[169] = ZERO;
  eqz[170] = -0.4e1 * P3[0][1] * P3[0][3] - 0.4e1 * P3[1][1] * P3[1][3];
  eqz[171] = -0.2e1 * P3[1][2] * P3[1][3] - 0.2e1 * P3[0][2] * P3[0][3];
  eqz[172] = ZERO;
  eqz[173] = 2 * P2[0][2] * P2[0][3] + 2 * P2[1][2] * P2[1][3];
  eqz[174] = ZERO;
  eqz[175] = ZERO;
  eqz[176] = ZERO;
  eqz[177] = ZERO;
  eqz[178] = ZERO;
  eqz[179] = ZERO;
  eqz[180] = ZERO;
  eqz[181] = ZERO;
  eqz[182] = 2 * P1[0][2] * P1[0][3] + 2 * P1[1][2] * P1[1][3];
  eqz[183] = ZERO;
  eqz[184] = ZERO;
  eqz[185] = ZERO;
  eqz[186] = ZERO;
  eqz[187] = ZERO;
  eqz[188] = -0.2e1 * pow(P3[1][3], 0.2e1) - 0.2e1 * pow(P3[0][3], 0.2e1);
  eqz[189] = ZERO;
  eqz[190] = ZERO;
  eqz[191] = ZERO;
  eqz[192] = ZERO;
  eqz[193] = ZERO;
  eqz[194] = ZERO;
  eqz[195] = ZERO;
  eqz[196] = ZERO;
  eqz[197] = ZERO;
  eqz[198] = ZERO;
  eqz[199] = ZERO;
  eqz[200] = ZERO;
  eqz[201] = ZERO;
  eqz[202] = ZERO;
  eqz[203] = ZERO;
  eqz[204] = ZERO;
  eqz[205] = ZERO;
  eqz[206] = ZERO;
  eqz[207] = ZERO;
  eqz[208] = ZERO;

  //EQT
  eqt[0] = ZERO;
  eqt[1] = ZERO;
  eqt[2] = ZERO;
  eqt[3] = ZERO;
  eqt[4] = ZERO;
  eqt[5] = ZERO;
  eqt[6] = ZERO;
  eqt[7] = ZERO;
  eqt[8] = ZERO;
  eqt[9] = ZERO;
  eqt[10] = ZERO;
  eqt[11] = ZERO;
  eqt[12] = ZERO;
  eqt[13] = ZERO;
  eqt[14] = ZERO;
  eqt[15] = ZERO;
  eqt[16] = ZERO;
  eqt[17] = ZERO;
  eqt[18] = ZERO;
  eqt[19] = ZERO;
  eqt[20] = ZERO;
  eqt[21] = ZERO;
  eqt[22] = ZERO;
  eqt[23] = ZERO;
  eqt[24] = ZERO;
  eqt[25] = ZERO;
  eqt[26] = ZERO;
  eqt[27] = ZERO;
  eqt[28] = ZERO;
  eqt[29] = ZERO;
  eqt[30] = ZERO;
  eqt[31] = ZERO;
  eqt[32] = ZERO;
  eqt[33] = ZERO;
  eqt[34] = ZERO;
  eqt[35] = ZERO;
  eqt[36] = ZERO;
  eqt[37] = ZERO;
  eqt[38] = ZERO;
  eqt[39] = ZERO;
  eqt[40] = ZERO;
  eqt[41] = ZERO;
  eqt[42] = ZERO;
  eqt[43] = ZERO;
  eqt[44] = ZERO;
  eqt[45] = ZERO;
  eqt[46] = ZERO;
  eqt[47] = ZERO;
  eqt[48] = ZERO;
  eqt[49] = ZERO;
  eqt[50] = ZERO;
  eqt[51] = ZERO;
  eqt[52] = ZERO;
  eqt[53] = ZERO;
  eqt[54] = ZERO;
  eqt[55] = ZERO;
  eqt[56] = ZERO;
  eqt[57] = ZERO;
  eqt[58] = ZERO;
  eqt[59] = ZERO;
  eqt[60] = ZERO;
  eqt[61] = ZERO;
  eqt[62] = ZERO;
  eqt[63] = ZERO;
  eqt[64] = ZERO;
  eqt[65] = ZERO;
  eqt[66] = ZERO;
  eqt[67] = ZERO;
  eqt[68] = ZERO;
  eqt[69] = ZERO;
  eqt[70] = ZERO;
  eqt[71] = ZERO;
  eqt[72] = ZERO;
  eqt[73] = ZERO;
  eqt[74] = ZERO;
  eqt[75] = ZERO;
  eqt[76] = ZERO;
  eqt[77] = ZERO;
  eqt[78] = ZERO;
  eqt[79] = ZERO;
  eqt[80] = ZERO;
  eqt[81] = ZERO;
  eqt[82] = ZERO;
  eqt[83] = ZERO;
  eqt[84] = ZERO;
  eqt[85] = ZERO;
  eqt[86] = ZERO;
  eqt[87] = ZERO;
  eqt[88] = ZERO;
  eqt[89] = ZERO;
  eqt[90] = ZERO;
  eqt[91] = ZERO;
  eqt[92] = ZERO;
  eqt[93] = ZERO;
  eqt[94] = ZERO;
  eqt[95] = ZERO;
  eqt[96] = ZERO;
  eqt[97] = ZERO;
  eqt[98] = ZERO;
  eqt[99] = ZERO;
  eqt[100] = ZERO;
  eqt[101] = ZERO;
  eqt[102] = ZERO;
  eqt[103] = ZERO;
  eqt[104] = ZERO;
  eqt[105] = ZERO;
  eqt[106] = ZERO;
  eqt[107] = ZERO;
  eqt[108] = ZERO;
  eqt[109] = ZERO;
  eqt[110] = ZERO;
  eqt[111] = ZERO;
  eqt[112] = ZERO;
  eqt[113] = ZERO;
  eqt[114] = ZERO;
  eqt[115] = ZERO;
  eqt[116] = ZERO;
  eqt[117] = ZERO;
  eqt[118] = ZERO;
  eqt[119] = ZERO;
  eqt[120] = ZERO;
  eqt[121] = ZERO;
  eqt[122] = ZERO;
  eqt[123] = ZERO;
  eqt[124] = ZERO;
  eqt[125] = ZERO;
  eqt[126] = ZERO;
  eqt[127] = ZERO;
  eqt[128] = ZERO;
  eqt[129] = ZERO;
  eqt[130] = ZERO;
  eqt[131] = ZERO;
  eqt[132] = ZERO;
  eqt[133] = ZERO;
  eqt[134] = ZERO;
  eqt[135] = ZERO;
  eqt[136] = ZERO;
  eqt[137] = ZERO;
  eqt[138] = ZERO;
  eqt[139] = ZERO;
  eqt[140] = ZERO;
  eqt[141] = ZERO;
  eqt[142] = ZERO;
  eqt[143] = ZERO;
  eqt[144] = ZERO;
  eqt[145] = ZERO;
  eqt[146] = ZERO;
  eqt[147] = ZERO;
  eqt[148] = ZERO;
  eqt[149] = ZERO;
  eqt[150] = ZERO;
  eqt[151] = ZERO;
  eqt[152] = ZERO;
  eqt[153] = ZERO;
  eqt[154] = ZERO;
  eqt[155] = ZERO;
  eqt[156] = ZERO;
  eqt[157] = ZERO;
  eqt[158] = ZERO;
  eqt[159] = ZERO;
  eqt[160] = ZERO;
  eqt[161] = ZERO;
  eqt[162] = ZERO;
  eqt[163] = ZERO;
  eqt[164] = ZERO;
  eqt[165] = ZERO;
  eqt[166] = ZERO;
  eqt[167] = -2 * P3[1][0] * P3[1][3] - 2 * P3[0][0] * P3[0][3];
  eqt[168] = ZERO;
  eqt[169] = -2 * P2[1][0] * P2[1][3] - 2 * P2[0][0] * P2[0][3];
  eqt[170] = -2 * P3[0][1] * P3[0][3] - 2 * P3[1][1] * P3[1][3];
  eqt[171] = -2 * P3[1][2] * P3[1][3] - 2 * P3[0][2] * P3[0][3];
  eqt[172] = -2 * P2[0][1] * P2[0][3] - 2 * P2[1][1] * P2[1][3];
  eqt[173] = -2 * P2[1][2] * P2[1][3] - 2 * P2[0][2] * P2[0][3];
  eqt[174] = ZERO;
  eqt[175] = ZERO;
  eqt[176] = -2 * P1[1][0] * P1[1][3] - 2 * P1[0][0] * P1[0][3];
  eqt[177] = ZERO;
  eqt[178] = ZERO;
  eqt[179] = ZERO;
  eqt[180] = ZERO;
  eqt[181] = -2 * P1[0][1] * P1[0][3] - 2 * P1[1][1] * P1[1][3];
  eqt[182] = -2 * P1[0][2] * P1[0][3] - 2 * P1[1][2] * P1[1][3];
  eqt[183] = ZERO;
  eqt[184] = ZERO;
  eqt[185] = ZERO;
  eqt[186] = ZERO;
  eqt[187] = ZERO;
  eqt[188] = -2 * pow((double) P3[1][3], (double) 2) - 2 * pow((double) P3[0][3], (double) 2);
  eqt[189] = ZERO;
  eqt[190] = -2 * pow((double) P2[0][3], (double) 2) - 2 * pow((double) P2[1][3], (double) 2);
  eqt[191] = ZERO;
  eqt[192] = ZERO;
  eqt[193] = ZERO;
  eqt[194] = ZERO;
  eqt[195] = ZERO;
  eqt[196] = ZERO;
  eqt[197] = -2 * pow((double) P1[1][3], (double) 2) - 2 * pow((double) P1[0][3], (double) 2);
  eqt[198] = ZERO;
  eqt[199] = ZERO;
  eqt[200] = ZERO;
  eqt[201] = ZERO;
  eqt[202] = ZERO;
  eqt[203] = ZERO;
  eqt[204] = ZERO;
  eqt[205] = ZERO;
  eqt[206] = ZERO;
  eqt[207] = ZERO;
  eqt[208] = ZERO;
}
#endif