Home
last modified time | relevance | path

Searched refs:Bnum (Results 1 – 25 of 76) sorted by relevance

1234

/dports/science/openmx/openmx3.8/source/
H A DTRAN_Set_SurfOverlap_NC.c55 int Anum, Bnum, wanA, tnoA; in TRAN_Set_SurfOverlap_NC() local
119 Bnum = MP[GB_AN]; in TRAN_Set_SurfOverlap_NC()
149 H00_nc_e[iside][S00_ref(Anum+i,Bnum+j)].r += co*a00 - si*b00; in TRAN_Set_SurfOverlap_NC()
150 H00_nc_e[iside][S00_ref(Anum+i,Bnum+j)].i += co*b00 + si*a00; in TRAN_Set_SurfOverlap_NC()
152 H00_nc_e[iside][S01_ref(Anum+i,Bnum+j)].r += co*a01r - si*(a01i+b01); in TRAN_Set_SurfOverlap_NC()
155 H00_nc_e[iside][S11_ref(Anum+i,Bnum+j)].r += co*a11 - si*b11; in TRAN_Set_SurfOverlap_NC()
156 H00_nc_e[iside][S11_ref(Anum+i,Bnum+j)].i += co*b11 + si*a11; in TRAN_Set_SurfOverlap_NC()
193 H01_nc_e[iside][S00_ref(Anum+i,Bnum+j)].r += co*a00 - si*b00; in TRAN_Set_SurfOverlap_NC()
194 H01_nc_e[iside][S00_ref(Anum+i,Bnum+j)].i += co*b00 + si*a00; in TRAN_Set_SurfOverlap_NC()
199 H01_nc_e[iside][S11_ref(Anum+i,Bnum+j)].r += co*a11 - si*b11; in TRAN_Set_SurfOverlap_NC()
[all …]
H A DTRAN_Set_CentOverlap_NC.c144 int LB_AN, GB_AN, wanB, tnoB, l1,l2,l3, Bnum; in TRAN_Set_CentOverlap_NC() local
176 Bnum = MP[GB_AN]; in TRAN_Set_CentOverlap_NC()
194 SCC_nc[SCC00_ref(Anum+i,Bnum+j)].r += co*S1[q]; in TRAN_Set_CentOverlap_NC()
195 SCC_nc[SCC00_ref(Anum+i,Bnum+j)].i += si*S1[q]; in TRAN_Set_CentOverlap_NC()
197 SCC_nc[SCC11_ref(Anum+i,Bnum+j)].r += co*S1[q]; in TRAN_Set_CentOverlap_NC()
198 SCC_nc[SCC11_ref(Anum+i,Bnum+j)].i += si*S1[q]; in TRAN_Set_CentOverlap_NC()
210 HCC_nc[SCC00_ref(Anum+i,Bnum+j)].r += co*a00 - si*b00; in TRAN_Set_CentOverlap_NC()
255 int GB_AN, wanB, tnoB, Bnum; in TRAN_Set_CentOverlap_NC() local
281 Bnum = MP[GB_AN]; in TRAN_Set_CentOverlap_NC()
334 int GB_AN, wanB, tnoB, Bnum; in TRAN_Set_CentOverlap_NC() local
[all …]
H A DHamiltonian_Cluster_NC.c279 Bnum = MP[GB_AN]; in Hamiltonian_Cluster_NC()
281 H[Anum+i][Bnum+j].r += H1[k]; in Hamiltonian_Cluster_NC()
303 Bnum = MP[GB_AN]; in Hamiltonian_Cluster_NC()
327 Bnum = MP[GB_AN]; in Hamiltonian_Cluster_NC()
352 Bnum = MP[GB_AN]; in Hamiltonian_Cluster_NC()
386 Bnum = MP[GB_AN]; in Hamiltonian_Cluster_NC()
410 Bnum = MP[GB_AN]; in Hamiltonian_Cluster_NC()
434 Bnum = MP[GB_AN]; in Hamiltonian_Cluster_NC()
459 Bnum = MP[GB_AN]; in Hamiltonian_Cluster_NC()
483 Bnum = MP[GB_AN]; in Hamiltonian_Cluster_NC()
[all …]
H A DTRAN_Set_SurfOverlap.c51 int Anum, Bnum, wanA, tnoA; in TRAN_Set_SurfOverlap() local
115 Bnum = MP[GB_AN]; in TRAN_Set_SurfOverlap()
129 S00_e[iside][S00_ref(Anum+i,Bnum+j)].r += co*OLP_e[iside][0][GA_AN][LB_AN][i][j]; in TRAN_Set_SurfOverlap()
130 S00_e[iside][S00_ref(Anum+i,Bnum+j)].i += si*OLP_e[iside][0][GA_AN][LB_AN][i][j]; in TRAN_Set_SurfOverlap()
133 H00_e[iside][k][S00_ref(Anum+i,Bnum+j)].r += co*H_e[iside][k][GA_AN][LB_AN][i][j]; in TRAN_Set_SurfOverlap()
134 H00_e[iside][k][S00_ref(Anum+i,Bnum+j)].i += si*H_e[iside][k][GA_AN][LB_AN][i][j]; in TRAN_Set_SurfOverlap()
145 S01_e[iside][S00_ref(Anum+i,Bnum+j)].r += co*OLP_e[iside][0][GA_AN][LB_AN][i][j]; in TRAN_Set_SurfOverlap()
146 S01_e[iside][S00_ref(Anum+i,Bnum+j)].i += si*OLP_e[iside][0][GA_AN][LB_AN][i][j]; in TRAN_Set_SurfOverlap()
149 H01_e[iside][k][S00_ref(Anum+i,Bnum+j)].r += co*H_e[iside][k][GA_AN][LB_AN][i][j]; in TRAN_Set_SurfOverlap()
150 H01_e[iside][k][S00_ref(Anum+i,Bnum+j)].i += si*H_e[iside][k][GA_AN][LB_AN][i][j]; in TRAN_Set_SurfOverlap()
H A DTRAN_Set_CentOverlap.c111 int LB_AN, GB_AN, wanB, tnoB, l1,l2,l3, Bnum; in TRAN_Set_CentOverlap() local
144 Bnum = MP[GB_AN]; in TRAN_Set_CentOverlap()
161 SCC[SCC_ref(Anum+i,Bnum+j)].r += co*S1[q]; in TRAN_Set_CentOverlap()
162 SCC[SCC_ref(Anum+i,Bnum+j)].i += si*S1[q]; in TRAN_Set_CentOverlap()
165 HCC[k][SCC_ref(Anum+i,Bnum+j)].r += co*H1[k][q]; in TRAN_Set_CentOverlap()
166 HCC[k][SCC_ref(Anum+i,Bnum+j)].i += si*H1[k][q]; in TRAN_Set_CentOverlap()
186 int GB_AN, wanB, tnoB, Bnum; in TRAN_Set_CentOverlap() local
212 Bnum = MP[GB_AN]; in TRAN_Set_CentOverlap()
236 int GB_AN, wanB, tnoB, Bnum; in TRAN_Set_CentOverlap() local
264 Bnum = MP[GB_AN]; in TRAN_Set_CentOverlap()
[all …]
H A DTRAN_Main_Analysis_NC.c3590 int Anum, Bnum, wanA, tnoA; in MTRAN_Set_SurfOverlap() local
3658 Bnum = MP[GB_AN]; in MTRAN_Set_SurfOverlap()
3922 Bnum = MP[GB_AN]; in MTRAN_Set_CentOverlap()
4033 Bnum = MP[GB_AN]; in MTRAN_Set_CentOverlap()
4047 VCC[SCC00_ref(Anum + i, Bnum + j)].r = VCC[SCC00_ref(Anum + i, Bnum + j)].r in MTRAN_Set_CentOverlap()
4050 VCC[SCC01_ref(Anum + i, Bnum + j)].r = VCC[SCC01_ref(Anum + i, Bnum + j)].r in MTRAN_Set_CentOverlap()
4053 VCC[SCC10_ref(Anum + i, Bnum + j)].r = VCC[SCC10_ref(Anum + i, Bnum + j)].r in MTRAN_Set_CentOverlap()
4056 VCC[SCC11_ref(Anum + i, Bnum + j)].r = VCC[SCC11_ref(Anum + i, Bnum + j)].r in MTRAN_Set_CentOverlap()
4074 VCC[SCC00_ref(Anum + i, Bnum + j)].i = VCC[SCC00_ref(Anum + i, Bnum + j)].i in MTRAN_Set_CentOverlap()
4077 VCC[SCC01_ref(Anum + i, Bnum + j)].i = VCC[SCC01_ref(Anum + i, Bnum + j)].i in MTRAN_Set_CentOverlap()
[all …]
H A DHamiltonian_Band_NC.c25 int i,j,k,wanA,wanB,tnoA,tnoB,Anum,Bnum; in Hamiltonian_Band_NC() local
147 Bnum = MP[GB_AN]; in Hamiltonian_Band_NC()
155 H11r[Anum+i][Bnum+j] += co*RH[0][MA_AN][LB_AN][i][j]; in Hamiltonian_Band_NC()
156 H11i[Anum+i][Bnum+j] += si*RH[0][MA_AN][LB_AN][i][j]; in Hamiltonian_Band_NC()
157 H22r[Anum+i][Bnum+j] += co*RH[1][MA_AN][LB_AN][i][j]; in Hamiltonian_Band_NC()
158 H22i[Anum+i][Bnum+j] += si*RH[1][MA_AN][LB_AN][i][j]; in Hamiltonian_Band_NC()
159 H12r[Anum+i][Bnum+j] += co*RH[2][MA_AN][LB_AN][i][j] - si*RH[3][MA_AN][LB_AN][i][j]; in Hamiltonian_Band_NC()
160 H12i[Anum+i][Bnum+j] += si*RH[2][MA_AN][LB_AN][i][j] + co*RH[3][MA_AN][LB_AN][i][j]; in Hamiltonian_Band_NC()
169 H11r[Anum+i][Bnum+j] += co*RH[0][MA_AN][LB_AN][i][j] - si*IH[0][MA_AN][LB_AN][i][j]; in Hamiltonian_Band_NC()
170 H11i[Anum+i][Bnum+j] += si*RH[0][MA_AN][LB_AN][i][j] + co*IH[0][MA_AN][LB_AN][i][j]; in Hamiltonian_Band_NC()
[all …]
H A DTRAN_DFT_Dosout.c285 Bnum = MP[GB_AN]; in TRAN_DFT_Kdependent_NC()
308 tmpr =-0.5*(GCR[ (Bnum+j-1)*NUM_c+(Anum+i)-1 ].i - GCA[ (Bnum+j-1)*NUM_c+(Anum+i)-1 ].i); in TRAN_DFT_Kdependent_NC()
309 tmpi = 0.5*(GCR[ (Bnum+j-1)*NUM_c+(Anum+i)-1 ].r - GCA[ (Bnum+j-1)*NUM_c+(Anum+i)-1 ].r); in TRAN_DFT_Kdependent_NC()
356 Bnum = MP_e[iside][GB_AN_e]; in TRAN_DFT_Kdependent_NC()
417 Bnum = MP_e[iside][GB_AN_e]; in TRAN_DFT_Kdependent_NC()
1198 Bnum = MP[GB_AN]; in TRAN_DFT_Kdependent_Col()
1234 tmpr =-0.5*(GCR[ v_idx( Anum+i, Bnum+j ) ].i - GCA[ v_idx( Anum+i, Bnum+j ) ].i); in TRAN_DFT_Kdependent_Col()
1235 tmpi = 0.5*(GCR[ v_idx( Anum+i, Bnum+j ) ].r - GCA[ v_idx( Anum+i, Bnum+j ) ].r); in TRAN_DFT_Kdependent_Col()
1278 Bnum = MP_e[iside][GB_AN_e]; in TRAN_DFT_Kdependent_Col()
1295 tmpr =-0.5*(GCL_R[ v_idx( Anum+i, Bnum+j) ].i - GCL_A[ v_idx( Anum+i, Bnum+j) ].i); in TRAN_DFT_Kdependent_Col()
[all …]
H A DDivide_Conquer_Dosout.c640 Bnum = MP[j]; in DC_Dosout_Col()
661 S_DC[Anum+m][Bnum+n] = 0.0; in DC_Dosout_Col()
668 H_DC[spin][Anum+m][Bnum+n] = 0.0; in DC_Dosout_Col()
784 Bnum = MP[h_AN]; in DC_Dosout_Col()
1590 Bnum = MP[j]; in DC_Dosout_NonCol()
1610 H_DC[Bnum+n+NUM][Anum+m ].r = H_DC[Anum+m ][Bnum+n+NUM].r; in DC_Dosout_NonCol()
1611 H_DC[Bnum+n+NUM][Anum+m ].i = -H_DC[Anum+m ][Bnum+n+NUM].i; in DC_Dosout_NonCol()
1626 H_DC[Bnum+n+NUM][Anum+m ].r = H_DC[Anum+m ][Bnum+n+NUM].r; in DC_Dosout_NonCol()
1627 H_DC[Bnum+n+NUM][Anum+m ].i = -H_DC[Anum+m ][Bnum+n+NUM].i; in DC_Dosout_NonCol()
1638 S_DC[Anum+m][Bnum+n] = 0.0; in DC_Dosout_NonCol()
[all …]
H A DDivide_Conquer.c732 Bnum = MP[j]; in DC_Col()
758 S_DC[Anum+m][Bnum+n] = 0.0; in DC_Col()
766 H_DC[spin][Anum+m][Bnum+n] = 0.0; in DC_Col()
1040 Bnum = MP[h_AN]; in DC_Col()
2274 Bnum = MP[j]; in DC_NonCol()
2297 H_DC[Bnum+n+NUM][Anum+m ].r = H_DC[Anum+m ][Bnum+n+NUM].r; in DC_NonCol()
2298 H_DC[Bnum+n+NUM][Anum+m ].i = -H_DC[Anum+m ][Bnum+n+NUM].i; in DC_NonCol()
2313 H_DC[Bnum+n+NUM][Anum+m ].r = H_DC[Anum+m ][Bnum+n+NUM].r; in DC_NonCol()
2314 H_DC[Bnum+n+NUM][Anum+m ].i = -H_DC[Anum+m ][Bnum+n+NUM].i; in DC_NonCol()
2326 S_DC[Anum+m][Bnum+n] = 0.0; in DC_NonCol()
[all …]
H A DCluster_DFT_ScaLAPACK.c207 int MA_AN,LB_AN,GB_AN,Bnum,MaxN; in Cluster_collinear() local
1013 Bnum = MP[GB_AN]; in Cluster_collinear()
1064 Bnum = MP[GB_AN]; in Cluster_collinear()
2146 Bnum = MP[GB_AN]; in Cluster_non_collinear()
2209 Bnum = MP[GB_AN]; in Cluster_non_collinear()
2678 Bnum = MP[GB_AN]; in Save_DOS_Col()
2931 Bnum = MP[GB_AN]; in Save_DOS_NonCol()
2936 dum = C[k][Anum+i].r*C[k][Bnum+j].r + C[k][Anum+i].i*C[k][Bnum+j].i; in Save_DOS_NonCol()
2940 dum = C[k][Anum+i+n].r*C[k][Bnum+j+n].r + C[k][Anum+i+n].i*C[k][Bnum+j+n].i; in Save_DOS_NonCol()
2944 dum = C[k][Anum+i].r*C[k][Bnum+j+n].r + C[k][Anum+i].i*C[k][Bnum+j+n].i; in Save_DOS_NonCol()
[all …]
H A DCluster_DFT.c143 int MA_AN,LB_AN,GB_AN,Bnum,MaxN; in Cluster_collinear() local
1074 Bnum = MP[GB_AN]; in Cluster_collinear()
1125 Bnum = MP[GB_AN]; in Cluster_collinear()
2215 Bnum = MP[GB_AN]; in Cluster_non_collinear()
2278 Bnum = MP[GB_AN]; in Cluster_non_collinear()
2747 Bnum = MP[GB_AN]; in Save_DOS_Col()
3000 Bnum = MP[GB_AN]; in Save_DOS_NonCol()
3005 dum = C[k][Anum+i].r*C[k][Bnum+j].r + C[k][Anum+i].i*C[k][Bnum+j].i; in Save_DOS_NonCol()
3009 dum = C[k][Anum+i+n].r*C[k][Bnum+j+n].r + C[k][Anum+i+n].i*C[k][Bnum+j+n].i; in Save_DOS_NonCol()
3013 dum = C[k][Anum+i].r*C[k][Bnum+j+n].r + C[k][Anum+i].i*C[k][Bnum+j+n].i; in Save_DOS_NonCol()
[all …]
H A DpolB.c1080 int Anum,Bnum,tnoA,tnoB; in Overlap_k1k2() local
1481 Bnum = MP[Gh_AN]; in Overlap_k1k2()
1573 Bnum = MP[Gh_AN]; in Overlap_k1k2()
1740 Bnum = MP[GB_AN]; in Overlap_Band()
1744 S1[Anum+i][Bnum+j] += s*co; in Overlap_Band()
1745 S2[Anum+i][Bnum+j] += s*si; in Overlap_Band()
1834 Bnum = MP[GB_AN]; in Hamiltonian_Band()
1838 H1[Anum+i][Bnum+j] += h*co; in Hamiltonian_Band()
1839 H2[Anum+i][Bnum+j] += h*si; in Hamiltonian_Band()
1884 static int i,j,k,wanA,wanB,tnoA,tnoB,Anum,Bnum; in Hamiltonian_Band_NC() local
[all …]
/dports/math/e-antic/flint2-ae7ec89/fmpq_mat/
H A Dsolve_fraction_free.c43 fmpz_mat_t Bnum; in fmpq_mat_solve_fraction_free() local
49 fmpz_mat_init(Bnum, B->r, B->c); in fmpq_mat_solve_fraction_free()
53 fmpq_mat_get_fmpz_mat_rowwise_2(Anum, Bnum, NULL, A, B); in fmpq_mat_solve_fraction_free()
56 success = fmpz_mat_solve(Xnum, den, Anum, Bnum); in fmpq_mat_solve_fraction_free()
62 fmpz_mat_clear(Bnum); in fmpq_mat_solve_fraction_free()
H A Dsolve_multi_mod.c185 fmpz_mat_t Bnum; in fmpq_mat_solve_multi_mod() local
189 fmpz_mat_init(Bnum, B->r, B->c); in fmpq_mat_solve_multi_mod()
191 fmpq_mat_get_fmpz_mat_rowwise_2(Anum, Bnum, NULL, A, B); in fmpq_mat_solve_multi_mod()
192 success = fmpq_mat_solve_fmpz_mat_multi_mod(X, Anum, Bnum); in fmpq_mat_solve_multi_mod()
195 fmpz_mat_clear(Bnum); in fmpq_mat_solve_multi_mod()
H A Dsolve_dixon.c196 fmpz_mat_t Bnum; in fmpq_mat_solve_dixon() local
200 fmpz_mat_init(Bnum, B->r, B->c); in fmpq_mat_solve_dixon()
202 fmpq_mat_get_fmpz_mat_rowwise_2(Anum, Bnum, NULL, A, B); in fmpq_mat_solve_dixon()
203 success = fmpq_mat_solve_fmpz_mat_dixon(X, Anum, Bnum); in fmpq_mat_solve_dixon()
206 fmpz_mat_clear(Bnum); in fmpq_mat_solve_dixon()
/dports/math/flint2/flint-2.8.4/fmpq_mat/
H A Dsolve_fraction_free.c43 fmpz_mat_t Bnum; in fmpq_mat_solve_fraction_free() local
49 fmpz_mat_init(Bnum, B->r, B->c); in fmpq_mat_solve_fraction_free()
53 fmpq_mat_get_fmpz_mat_rowwise_2(Anum, Bnum, NULL, A, B); in fmpq_mat_solve_fraction_free()
56 success = fmpz_mat_solve(Xnum, den, Anum, Bnum); in fmpq_mat_solve_fraction_free()
62 fmpz_mat_clear(Bnum); in fmpq_mat_solve_fraction_free()
H A Dsolve_multi_mod.c185 fmpz_mat_t Bnum; in fmpq_mat_solve_multi_mod() local
189 fmpz_mat_init(Bnum, B->r, B->c); in fmpq_mat_solve_multi_mod()
191 fmpq_mat_get_fmpz_mat_rowwise_2(Anum, Bnum, NULL, A, B); in fmpq_mat_solve_multi_mod()
192 success = fmpq_mat_solve_fmpz_mat_multi_mod(X, Anum, Bnum); in fmpq_mat_solve_multi_mod()
195 fmpz_mat_clear(Bnum); in fmpq_mat_solve_multi_mod()
H A Dsolve_dixon.c196 fmpz_mat_t Bnum; in fmpq_mat_solve_dixon() local
200 fmpz_mat_init(Bnum, B->r, B->c); in fmpq_mat_solve_dixon()
202 fmpq_mat_get_fmpz_mat_rowwise_2(Anum, Bnum, NULL, A, B); in fmpq_mat_solve_dixon()
203 success = fmpq_mat_solve_fmpz_mat_dixon(X, Anum, Bnum); in fmpq_mat_solve_dixon()
206 fmpz_mat_clear(Bnum); in fmpq_mat_solve_dixon()
/dports/math/e-antic/e-antic-1.0.0-rc.13/libeantic/upstream/antic/fmpq_mat/
H A Dsolve_fraction_free.c43 fmpz_mat_t Bnum; in fmpq_mat_solve_fraction_free() local
49 fmpz_mat_init(Bnum, B->r, B->c); in fmpq_mat_solve_fraction_free()
53 fmpq_mat_get_fmpz_mat_rowwise_2(Anum, Bnum, NULL, A, B); in fmpq_mat_solve_fraction_free()
56 success = fmpz_mat_solve(Xnum, den, Anum, Bnum); in fmpq_mat_solve_fraction_free()
62 fmpz_mat_clear(Bnum); in fmpq_mat_solve_fraction_free()
H A Dsolve_dixon.c196 fmpz_mat_t Bnum; in fmpq_mat_solve_dixon() local
200 fmpz_mat_init(Bnum, B->r, B->c); in fmpq_mat_solve_dixon()
202 fmpq_mat_get_fmpz_mat_rowwise_2(Anum, Bnum, NULL, A, B); in fmpq_mat_solve_dixon()
203 success = fmpq_mat_solve_fmpz_mat_dixon(X, Anum, Bnum); in fmpq_mat_solve_dixon()
206 fmpz_mat_clear(Bnum); in fmpq_mat_solve_dixon()
H A Dsolve_multi_mod.c185 fmpz_mat_t Bnum; in fmpq_mat_solve_multi_mod() local
189 fmpz_mat_init(Bnum, B->r, B->c); in fmpq_mat_solve_multi_mod()
191 fmpq_mat_get_fmpz_mat_rowwise_2(Anum, Bnum, NULL, A, B); in fmpq_mat_solve_multi_mod()
192 success = fmpq_mat_solve_fmpz_mat_multi_mod(X, Anum, Bnum); in fmpq_mat_solve_multi_mod()
195 fmpz_mat_clear(Bnum); in fmpq_mat_solve_multi_mod()
/dports/japanese/sj3-lib/sj3-2.0.1.20/sj3/
H A Dedit.c89 while (j < cv->Bnum) { in Bbig()
133 if (j < cv->Bnum) in Bbig()
153 if (j == cv->Bnum) { in Bsmall()
154 cv->Bnum ++; in Bsmall()
155 cv->Bkettei[cv->Bnum] = Bend; in Bsmall()
164 if (j == cv->Bnum) in Bsmall()
168 for (k = cv->Bnum ; k > j ; k --) { in Bsmall()
181 cv->Bnum ++; in Bsmall()
189 if (j == cv->Bnum) in Bsmall()
260 if (j == cv->Bnum) { in Bdelete()
[all …]
H A Dhenkan.c128 cv->Bnum = 0; in exec_henkan()
173 cv->Bkettei[cv->Bnum] = Bhenkan; in exec_henkan()
204 cv->Bnum++; in exec_henkan()
207 cv->Bnum++; in exec_henkan()
210 cv->Bkettei[cv->Bnum] = Bend; in exec_henkan()
212 for (i = 0 ; i < cv->Bnum ; i ++) { in exec_henkan()
499 cv->Bnum = 1; in exec_muhenkan()
817 if (cv->CurBun >= cv->Bnum) { in Bright()
846 cv->CurBun = cv->Bnum - 1; in Bleft()
861 cv->CurBun = cv->Bnum - 1; in Blast()
[all …]
/dports/math/giacxcas/giac-1.6.0/check/
H A Dmorley_demo.cas2 local a1,a2,A,Anum,B,Bnum,C0,C,Cnum,P1,P1num,P2,P2num,P3,P3num,P4,P5,j,k,q,r,s,tmp,cpt,trequi;
11 Bnum:=evalf(B);
23 P1num[q] =< affix(single_inter(line(Anum[0],Anum[k]),line(Bnum[0],Bnum[j])));
24 P2num[q] =< affix(single_inter(line(Bnum[0],Bnum[k]),line(Cnum[0],Cnum[j])));

1234