ncgmat.pas 20 KB

123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544
  1. {
  2. Copyright (c) 1998-2002 by Florian Klaempfl
  3. Generate generic mathematical nodes
  4. This program is free software; you can redistribute it and/or modify
  5. it under the terms of the GNU General Public License as published by
  6. the Free Software Foundation; either version 2 of the License, or
  7. (at your option) any later version.
  8. This program is distributed in the hope that it will be useful,
  9. but WITHOUT ANY WARRANTY; without even the implied warranty of
  10. MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
  11. GNU General Public License for more details.
  12. You should have received a copy of the GNU General Public License
  13. along with this program; if not, write to the Free Software
  14. Foundation, Inc., 675 Mass Ave, Cambridge, MA 02139, USA.
  15. ****************************************************************************
  16. }
  17. unit ncgmat;
  18. {$i fpcdefs.inc}
  19. interface
  20. uses
  21. node,nmat,cpubase,cgbase;
  22. type
  23. tcgunaryminusnode = class(tunaryminusnode)
  24. protected
  25. { This routine is called to change the sign of the
  26. floating point value in the floating point
  27. register r.
  28. This routine should be overridden, since
  29. the generic version is not optimal at all. The
  30. generic version assumes that floating
  31. point values are stored in the register
  32. in IEEE-754 format.
  33. }
  34. procedure emit_float_sign_change(r: tregister; _size : tcgsize);virtual;
  35. {$ifdef SUPPORT_MMX}
  36. procedure second_mmx;virtual;abstract;
  37. {$endif SUPPORT_MMX}
  38. {$ifndef cpu64bitalu}
  39. procedure second_64bit;virtual;
  40. {$endif not cpu64bitalu}
  41. procedure second_integer;virtual;
  42. procedure second_float;virtual;
  43. public
  44. procedure pass_generate_code;override;
  45. end;
  46. tcgmoddivnode = class(tmoddivnode)
  47. procedure pass_generate_code;override;
  48. protected
  49. { This routine must do an actual 32-bit division, be it
  50. signed or unsigned. The result must set into the the
  51. @var(num) register.
  52. @param(signed Indicates if the division must be signed)
  53. @param(denum Register containing the denominator
  54. @param(num Register containing the numerator, will also receive result)
  55. The actual optimizations regarding shifts have already
  56. been done and emitted, so this should really a do a divide.
  57. }
  58. procedure emit_div_reg_reg(signed: boolean;denum,num : tregister);virtual;abstract;
  59. { This routine must do an actual 32-bit modulo, be it
  60. signed or unsigned. The result must set into the the
  61. @var(num) register.
  62. @param(signed Indicates if the modulo must be signed)
  63. @param(denum Register containing the denominator
  64. @param(num Register containing the numerator, will also receive result)
  65. The actual optimizations regarding shifts have already
  66. been done and emitted, so this should really a do a modulo.
  67. }
  68. procedure emit_mod_reg_reg(signed: boolean;denum,num : tregister);virtual;abstract;
  69. {$ifndef cpu64bitalu}
  70. { This routine must do an actual 64-bit division, be it
  71. signed or unsigned. The result must set into the the
  72. @var(num) register.
  73. @param(signed Indicates if the division must be signed)
  74. @param(denum Register containing the denominator
  75. @param(num Register containing the numerator, will also receive result)
  76. The actual optimizations regarding shifts have already
  77. been done and emitted, so this should really a do a divide.
  78. Currently, this routine should only be implemented on
  79. 64-bit systems, otherwise a helper is called in 1st pass.
  80. }
  81. procedure emit64_div_reg_reg(signed: boolean;denum,num : tregister64);virtual;
  82. {$endif not cpu64bitalu}
  83. end;
  84. tcgshlshrnode = class(tshlshrnode)
  85. {$ifndef cpu64bitalu}
  86. procedure second_64bit;virtual;
  87. {$endif not cpu64bitalu}
  88. procedure second_integer;virtual;
  89. procedure pass_generate_code;override;
  90. end;
  91. tcgnotnode = class(tnotnode)
  92. protected
  93. procedure second_boolean;virtual;abstract;
  94. {$ifdef SUPPORT_MMX}
  95. procedure second_mmx;virtual;abstract;
  96. {$endif SUPPORT_MMX}
  97. {$ifndef cpu64bitalu}
  98. procedure second_64bit;virtual;
  99. {$endif not cpu64bitalu}
  100. procedure second_integer;virtual;
  101. public
  102. procedure pass_generate_code;override;
  103. end;
  104. implementation
  105. uses
  106. globtype,systems,
  107. cutils,verbose,globals,
  108. symtable,symconst,symtype,symdef,aasmbase,aasmtai,aasmdata,aasmcpu,defutil,
  109. parabase,
  110. pass_2,
  111. ncon,
  112. tgobj,ncgutil,cgobj,cgutils,paramgr,hlcgobj
  113. {$ifndef cpu64bitalu}
  114. ,cg64f32
  115. {$endif not cpu64bitalu}
  116. ;
  117. {*****************************************************************************
  118. TCGUNARYMINUSNODE
  119. *****************************************************************************}
  120. procedure tcgunaryminusnode.emit_float_sign_change(r: tregister; _size : tcgsize);
  121. var
  122. href,
  123. href2 : treference;
  124. begin
  125. { get a temporary memory reference to store the floating
  126. point value
  127. }
  128. tg.gettemp(current_asmdata.CurrAsmList,tcgsize2size[_size],tcgsize2size[_size],tt_normal,href);
  129. { store the floating point value in the temporary memory area }
  130. cg.a_loadfpu_reg_ref(current_asmdata.CurrAsmList,_size,_size,r,href);
  131. { only single and double ieee are supported, for little endian
  132. the signed bit is in the second dword }
  133. href2:=href;
  134. case _size of
  135. OS_F64 :
  136. if target_info.endian = endian_little then
  137. inc(href2.offset,4);
  138. OS_F32 :
  139. ;
  140. else
  141. internalerror(200406021);
  142. end;
  143. { flip sign-bit (bit 31/63) of single/double }
  144. cg.a_op_const_ref(current_asmdata.CurrAsmList,OP_XOR,OS_32,aint($80000000),href2);
  145. cg.a_loadfpu_ref_reg(current_asmdata.CurrAsmList,_size,_size,href,r);
  146. tg.ungetiftemp(current_asmdata.CurrAsmList,href);
  147. end;
  148. {$ifndef cpu64bitalu}
  149. procedure tcgunaryminusnode.second_64bit;
  150. var
  151. tr: tregister;
  152. hl: tasmlabel;
  153. begin
  154. secondpass(left);
  155. location_reset(location,LOC_REGISTER,left.location.size);
  156. location.register64.reglo:=cg.getintregister(current_asmdata.CurrAsmList,OS_INT);
  157. location.register64.reghi:=cg.getintregister(current_asmdata.CurrAsmList,OS_INT);
  158. cg64.a_op64_loc_reg(current_asmdata.CurrAsmList,OP_NEG,OS_S64,
  159. left.location,joinreg64(location.register64.reglo,location.register64.reghi));
  160. { there's only overflow in case left was low(int64) -> -left = left }
  161. if (cs_check_overflow in current_settings.localswitches) then
  162. begin
  163. tr:=cg.getintregister(current_asmdata.CurrAsmList,OS_INT);
  164. cg.a_op_const_reg_reg(current_asmdata.CurrAsmList,OP_XOR,OS_INT,
  165. aint($80000000),location.register64.reghi,tr);
  166. cg.a_op_reg_reg(current_asmdata.CurrAsmList,OP_OR,OS_INT,
  167. location.register64.reglo,tr);
  168. current_asmdata.getjumplabel(hl);
  169. cg.a_cmp_const_reg_label(current_asmdata.CurrAsmList,OS_INT,OC_NE,0,tr,hl);
  170. cg.a_call_name(current_asmdata.CurrAsmList,'FPC_OVERFLOW',false);
  171. cg.a_label(current_asmdata.CurrAsmList,hl);
  172. end;
  173. end;
  174. {$endif not cpu64bitalu}
  175. procedure tcgunaryminusnode.second_float;
  176. begin
  177. secondpass(left);
  178. location_reset(location,LOC_FPUREGISTER,def_cgsize(resultdef));
  179. case left.location.loc of
  180. LOC_REFERENCE,
  181. LOC_CREFERENCE :
  182. begin
  183. location.register:=cg.getfpuregister(current_asmdata.CurrAsmList,location.size);
  184. cg.a_loadfpu_ref_reg(current_asmdata.CurrAsmList,
  185. left.location.size,location.size,
  186. left.location.reference,location.register);
  187. emit_float_sign_change(location.register,def_cgsize(left.resultdef));
  188. end;
  189. LOC_FPUREGISTER:
  190. begin
  191. location.register:=left.location.register;
  192. emit_float_sign_change(location.register,def_cgsize(left.resultdef));
  193. end;
  194. LOC_CFPUREGISTER:
  195. begin
  196. location.register:=cg.getfpuregister(current_asmdata.CurrAsmList,location.size);
  197. cg.a_loadfpu_reg_reg(current_asmdata.CurrAsmList,left.location.size,location.size,left.location.register,location.register);
  198. emit_float_sign_change(location.register,def_cgsize(left.resultdef));
  199. end;
  200. else
  201. internalerror(200306021);
  202. end;
  203. end;
  204. procedure tcgunaryminusnode.second_integer;
  205. var
  206. hl: tasmlabel;
  207. opsize: tdef;
  208. begin
  209. secondpass(left);
  210. { load left operator in a register }
  211. location_copy(location,left.location);
  212. {$ifdef cpunodefaultint}
  213. opsize:=left.resultdef;
  214. {$else cpunodefaultint}
  215. { in case of a 32 bit system that can natively execute 64 bit operations }
  216. if (left.resultdef.size<=sinttype.size) then
  217. opsize:=sinttype
  218. else
  219. opsize:={$ifdef cpu16bitalu}s32inttype{$else}s64inttype{$endif};
  220. {$endif cpunodefaultint}
  221. hlcg.location_force_reg(current_asmdata.CurrAsmList,location,left.resultdef,opsize,false);
  222. hlcg.a_op_reg_reg(current_asmdata.CurrAsmList,OP_NEG,opsize,location.register,location.register);
  223. if (cs_check_overflow in current_settings.localswitches) then
  224. begin
  225. current_asmdata.getjumplabel(hl);
  226. hlcg.a_cmp_const_reg_label(current_asmdata.CurrAsmList,opsize,OC_NE,torddef(opsize).low.svalue,location.register,hl);
  227. hlcg.g_call_system_proc(current_asmdata.CurrAsmList,'fpc_overflow',nil);
  228. hlcg.a_label(current_asmdata.CurrAsmList,hl);
  229. end;
  230. end;
  231. procedure tcgunaryminusnode.pass_generate_code;
  232. begin
  233. {$ifndef cpu64bitalu}
  234. if is_64bit(left.resultdef) then
  235. second_64bit
  236. else
  237. {$endif not cpu64bitalu}
  238. {$ifdef SUPPORT_MMX}
  239. if (cs_mmx in current_settings.localswitches) and is_mmx_able_array(left.resultdef) then
  240. second_mmx
  241. else
  242. {$endif SUPPORT_MMX}
  243. if (left.resultdef.typ=floatdef) then
  244. second_float
  245. else
  246. second_integer;
  247. end;
  248. {*****************************************************************************
  249. TCGMODDIVNODE
  250. *****************************************************************************}
  251. {$ifndef cpu64bitalu}
  252. procedure tcgmoddivnode.emit64_div_reg_reg(signed: boolean; denum,num:tregister64);
  253. begin
  254. { handled in pass_1 already, unless pass_1 is
  255. overridden
  256. }
  257. { should be handled in pass_1 (JM) }
  258. internalerror(200109052);
  259. end;
  260. {$endif not cpu64bitalu}
  261. procedure tcgmoddivnode.pass_generate_code;
  262. var
  263. hreg1 : tregister;
  264. hdenom : tregister;
  265. power : longint;
  266. hl : tasmlabel;
  267. paraloc1 : tcgpara;
  268. opsize : tcgsize;
  269. opdef : tdef;
  270. begin
  271. secondpass(left);
  272. if codegenerror then
  273. exit;
  274. secondpass(right);
  275. if codegenerror then
  276. exit;
  277. location_copy(location,left.location);
  278. {$ifndef cpu64bitalu}
  279. if is_64bit(resultdef) then
  280. begin
  281. if is_signed(left.resultdef) then
  282. opdef:=s64inttype
  283. else
  284. opdef:=u64inttype;
  285. { this code valid for 64-bit cpu's only ,
  286. otherwise helpers are called in pass_1
  287. }
  288. hlcg.location_force_reg(current_asmdata.CurrAsmList,location,left.resultdef,opdef,false);
  289. location_copy(location,left.location);
  290. hlcg.location_force_reg(current_asmdata.CurrAsmList,right.location,right.resultdef,opdef,false);
  291. emit64_div_reg_reg(is_signed(left.resultdef),
  292. joinreg64(right.location.register64.reglo,right.location.register64.reghi),
  293. joinreg64(location.register64.reglo,location.register64.reghi));
  294. end
  295. else
  296. {$endif not cpu64bitalu}
  297. begin
  298. if is_signed(left.resultdef) then
  299. begin
  300. opsize:=OS_SINT;
  301. opdef:=ossinttype;
  302. end
  303. else
  304. begin
  305. opsize:=OS_INT;
  306. opdef:=osuinttype;
  307. end;
  308. { put numerator in register }
  309. hlcg.location_force_reg(current_asmdata.CurrAsmList,left.location,left.resultdef,opdef,false);
  310. hreg1:=left.location.register;
  311. if (nodetype=divn) and
  312. (right.nodetype=ordconstn) and
  313. ispowerof2(tordconstnode(right).value.svalue,power) then
  314. Begin
  315. { for signed numbers, the numerator must be adjusted before the
  316. shift instruction, but not wih unsigned numbers! Otherwise,
  317. "Cardinal($ffffffff) div 16" overflows! (JM) }
  318. If is_signed(left.resultdef) Then
  319. Begin
  320. current_asmdata.getjumplabel(hl);
  321. cg.a_cmp_const_reg_label(current_asmdata.CurrAsmList,OS_INT,OC_GT,0,hreg1,hl);
  322. if power=1 then
  323. cg.a_op_const_reg(current_asmdata.CurrAsmList,OP_ADD,OS_INT,1,hreg1)
  324. else
  325. cg.a_op_const_reg(current_asmdata.CurrAsmList,OP_ADD,OS_INT,Tordconstnode(right).value.svalue-1,hreg1);
  326. cg.a_label(current_asmdata.CurrAsmList,hl);
  327. cg.a_op_const_reg(current_asmdata.CurrAsmList,OP_SAR,OS_INT,power,hreg1);
  328. End
  329. Else { not signed }
  330. cg.a_op_const_reg(current_asmdata.CurrAsmList,OP_SHR,OS_INT,power,hreg1);
  331. End
  332. else
  333. begin
  334. { bring denominator to hdenom }
  335. { hdenom is always free, it's }
  336. { only used for temporary }
  337. { purposes }
  338. hdenom := cg.getintregister(current_asmdata.CurrAsmList,OS_INT);
  339. hlcg.a_load_loc_reg(current_asmdata.CurrAsmList,right.resultdef,osuinttype,right.location,hdenom);
  340. { verify if the divisor is zero, if so return an error
  341. immediately
  342. }
  343. current_asmdata.getjumplabel(hl);
  344. cg.a_cmp_const_reg_label(current_asmdata.CurrAsmList,OS_INT,OC_NE,0,hdenom,hl);
  345. paraloc1.init;
  346. paramanager.getintparaloc(search_system_proc('fpc_handleerror'),1,paraloc1);
  347. cg.a_load_const_cgpara(current_asmdata.CurrAsmList,OS_S32,aint(200),paraloc1);
  348. paramanager.freecgpara(current_asmdata.CurrAsmList,paraloc1);
  349. cg.a_call_name(current_asmdata.CurrAsmList,'FPC_HANDLEERROR',false);
  350. paraloc1.done;
  351. cg.a_label(current_asmdata.CurrAsmList,hl);
  352. if nodetype = modn then
  353. emit_mod_reg_reg(is_signed(left.resultdef),hdenom,hreg1)
  354. else
  355. emit_div_reg_reg(is_signed(left.resultdef),hdenom,hreg1);
  356. end;
  357. location_reset(location,LOC_REGISTER,opsize);
  358. location.register:=hreg1;
  359. end;
  360. cg.g_overflowcheck(current_asmdata.CurrAsmList,location,resultdef);
  361. end;
  362. {*****************************************************************************
  363. TCGSHLRSHRNODE
  364. *****************************************************************************}
  365. {$ifndef cpu64bitalu}
  366. procedure tcgshlshrnode.second_64bit;
  367. begin
  368. { already hanled in 1st pass }
  369. internalerror(2002081501);
  370. end;
  371. {$endif not cpu64bitalu}
  372. procedure tcgshlshrnode.second_integer;
  373. var
  374. op : topcg;
  375. opdef : tdef;
  376. hcountreg : tregister;
  377. opsize : tcgsize;
  378. begin
  379. { determine operator }
  380. case nodetype of
  381. shln: op:=OP_SHL;
  382. shrn: op:=OP_SHR;
  383. end;
  384. {$ifdef cpunodefaultint}
  385. opsize:=left.location.size;
  386. opdef:=left.resultdef;
  387. {$else cpunodefaultint}
  388. { load left operators in a register }
  389. if is_signed(left.resultdef) then
  390. begin
  391. opsize:=OS_SINT;
  392. opdef:=ossinttype
  393. end
  394. else
  395. begin
  396. opsize:=OS_INT;
  397. opdef:=osuinttype;
  398. end;
  399. {$endif cpunodefaultint}
  400. hlcg.location_force_reg(current_asmdata.CurrAsmList,left.location,left.resultdef,opdef,true);
  401. location_reset(location,LOC_REGISTER,opsize);
  402. location.register:=cg.getintregister(current_asmdata.CurrAsmList,opsize);
  403. { shifting by a constant directly coded: }
  404. if (right.nodetype=ordconstn) then
  405. begin
  406. { l shl 32 should 0 imho, but neither TP nor Delphi do it in this way (FK)
  407. if right.value<=31 then
  408. }
  409. cg.a_op_const_reg_reg(current_asmdata.CurrAsmList,op,location.size,
  410. tordconstnode(right).value.uvalue and 31,left.location.register,location.register);
  411. {
  412. else
  413. emit_reg_reg(A_XOR,S_L,hregister1,
  414. hregister1);
  415. }
  416. end
  417. else
  418. begin
  419. { load right operators in a register - this
  420. is done since most target cpu which will use this
  421. node do not support a shift count in a mem. location (cec)
  422. }
  423. if not(right.location.loc in [LOC_CREGISTER,LOC_REGISTER]) then
  424. begin
  425. hcountreg:=cg.getintregister(current_asmdata.CurrAsmList,opsize);
  426. hlcg.a_load_loc_reg(current_asmdata.CurrAsmList,right.resultdef,opdef,right.location,hcountreg);
  427. end
  428. else
  429. hcountreg:=right.location.register;
  430. cg.a_op_reg_reg_reg(current_asmdata.CurrAsmList,op,opsize,hcountreg,left.location.register,location.register);
  431. end;
  432. end;
  433. procedure tcgshlshrnode.pass_generate_code;
  434. begin
  435. secondpass(left);
  436. secondpass(right);
  437. {$ifndef cpu64bitalu}
  438. if is_64bit(left.resultdef) then
  439. second_64bit
  440. else
  441. {$endif not cpu64bitalu}
  442. second_integer;
  443. end;
  444. {*****************************************************************************
  445. TCGNOTNODE
  446. *****************************************************************************}
  447. {$ifndef cpu64bitalu}
  448. procedure tcgnotnode.second_64bit;
  449. begin
  450. secondpass(left);
  451. hlcg.location_force_reg(current_asmdata.CurrAsmList,left.location,left.resultdef,left.resultdef,false);
  452. location_copy(location,left.location);
  453. { perform the NOT operation }
  454. cg64.a_op64_reg_reg(current_asmdata.CurrAsmList,OP_NOT,location.size,left.location.register64,location.register64);
  455. end;
  456. {$endif not cpu64bitalu}
  457. procedure tcgnotnode.second_integer;
  458. begin
  459. secondpass(left);
  460. hlcg.location_force_reg(current_asmdata.CurrAsmList,left.location,left.resultdef,left.resultdef,false);
  461. location_copy(location,left.location);
  462. { perform the NOT operation }
  463. hlcg.a_op_reg_reg(current_asmdata.CurrAsmList,OP_NOT,left.resultdef,location.register,location.register);
  464. end;
  465. procedure tcgnotnode.pass_generate_code;
  466. begin
  467. if is_boolean(resultdef) then
  468. second_boolean
  469. {$ifdef SUPPORT_MMX}
  470. else if (cs_mmx in current_settings.localswitches) and is_mmx_able_array(left.resultdef) then
  471. second_mmx
  472. {$endif SUPPORT_MMX}
  473. {$ifndef cpu64bitalu}
  474. else if is_64bit(left.resultdef) then
  475. second_64bit
  476. {$endif not cpu64bitalu}
  477. else
  478. second_integer;
  479. end;
  480. begin
  481. cmoddivnode:=tcgmoddivnode;
  482. cunaryminusnode:=tcgunaryminusnode;
  483. cshlshrnode:=tcgshlshrnode;
  484. cnotnode:=tcgnotnode;
  485. end.