You can not select more than 25 topics Topics must start with a chinese character,a letter or number, can include dashes ('-') and can be up to 35 characters long.

setparam-ref.c 36 kB

12 years ago
12 years ago
123456789101112131415161718192021222324252627282930313233343536373839404142434445464748495051525354555657585960616263646566676869707172737475767778798081828384858687888990919293949596979899100101102103104105106107108109110111112113114115116117118119120121122123124125126127128129130131132133134135136137138139140141142143144145146147148149150151152153154155156157158159160161162163164165166167168169170171172173174175176177178179180181182183184185186187188189190191192193194195196197198199200201202203204205206207208209210211212213214215216217218219220221222223224225226227228229230231232233234235236237238239240241242243244245246247248249250251252253254255256257258259260261262263264265266267268269270271272273274275276277278279280281282283284285286287288289290291292293294295296297298299300301302303304305306307308309310311312313314315316317318319320321322323324325326327328329330331332333334335336337338339340341342343344345346347348349350351352353354355356357358359360361362363364365366367368369370371372373374375376377378379380381382383384385386387388389390391392393394395396397398399400401402403404405406407408409410411412413414415416417418419420421422423424425426427428429430431432433434435436437438439440441442443444445446447448449450451452453454455456457458459460461462463464465466467468469470471472473474475476477478479480481482483484485486487488489490491492493494495496497498499500501502503504505506507508509510511512513514515516517518519520521522523524525526527528529530531532533534535536537538539540541542543544545546547548549550551552553554555556557558559560561562563564565566567568569570571572573574575576577578579580581582583584585586587588589590591592593594595596597598599600601602603604605606607608609610611612613614615616617618619620621622623624625626627628629630631632633634635636637638639640641642643644645646647648649650651652653654655656657658659660661662663664665666667668669670671672673674675676677678679680681682683684685686687688689690691692693694695696697698699700701702703704705706707708709710711712713714715716717718719720721722723724725726727728729730731732733734735736737738739740741742743744745746747748749750751752753754755756757758759760761762763764765766767768769770771772773774775776777778779780781782783784785786787788789790791792793794795796797798799800801802803804805806807808809810811812813814815816817818819820821822823824825826827828829830831832833834835836837838839840841842843844845846847848849850851852853854855856857858859860861862863864865866867868869870871872873874875876877878879880881882883884885886887888889890891892893894895896897898899900901902903904905906907908909910911912913914915916917918919920921922923924925926927928929930931932933934935936937938939940941942943944945946947948949950951952953954955956957958959960961962963964965966967968969970971972973974975976977978979980981982983984985986987988989990991992993994995996997998999100010011002100310041005100610071008100910101011101210131014101510161017101810191020102110221023102410251026102710281029103010311032103310341035103610371038103910401041104210431044104510461047104810491050105110521053105410551056105710581059106010611062106310641065106610671068106910701071107210731074107510761077107810791080108110821083108410851086108710881089109010911092109310941095109610971098109911001101110211031104110511061107110811091110111111121113111411151116111711181119112011211122112311241125112611271128112911301131113211331134113511361137113811391140114111421143114411451146114711481149115011511152115311541155115611571158115911601161116211631164116511661167116811691170117111721173117411751176117711781179118011811182118311841185118611871188118911901191119211931194119511961197119811991200120112021203120412051206120712081209121012111212121312141215121612171218121912201221122212231224122512261227122812291230123112321233
  1. /*********************************************************************/
  2. /* Copyright 2009, 2010 The University of Texas at Austin. */
  3. /* All rights reserved. */
  4. /* */
  5. /* Redistribution and use in source and binary forms, with or */
  6. /* without modification, are permitted provided that the following */
  7. /* conditions are met: */
  8. /* */
  9. /* 1. Redistributions of source code must retain the above */
  10. /* copyright notice, this list of conditions and the following */
  11. /* disclaimer. */
  12. /* */
  13. /* 2. Redistributions in binary form must reproduce the above */
  14. /* copyright notice, this list of conditions and the following */
  15. /* disclaimer in the documentation and/or other materials */
  16. /* provided with the distribution. */
  17. /* */
  18. /* THIS SOFTWARE IS PROVIDED BY THE UNIVERSITY OF TEXAS AT */
  19. /* AUSTIN ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, */
  20. /* INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF */
  21. /* MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE */
  22. /* DISCLAIMED. IN NO EVENT SHALL THE UNIVERSITY OF TEXAS AT */
  23. /* AUSTIN OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, */
  24. /* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES */
  25. /* (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE */
  26. /* GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR */
  27. /* BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF */
  28. /* LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT */
  29. /* (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT */
  30. /* OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE */
  31. /* POSSIBILITY OF SUCH DAMAGE. */
  32. /* */
  33. /* The views and conclusions contained in the software and */
  34. /* documentation are those of the authors and should not be */
  35. /* interpreted as representing official policies, either expressed */
  36. /* or implied, of The University of Texas at Austin. */
  37. /*********************************************************************/
  38. #include <stdio.h>
  39. #include <string.h>
  40. #include "common.h"
  41. #ifdef BUILD_KERNEL
  42. #include "kernelTS.h"
  43. #endif
  44. #undef DEBUG
  45. static void init_parameter(void);
  46. gotoblas_t TABLE_NAME = {
  47. DTB_DEFAULT_ENTRIES ,
  48. GEMM_DEFAULT_OFFSET_A, GEMM_DEFAULT_OFFSET_B, GEMM_DEFAULT_ALIGN,
  49. 0, 0, 0,
  50. SGEMM_DEFAULT_UNROLL_M, SGEMM_DEFAULT_UNROLL_N,
  51. #ifdef SGEMM_DEFAULT_UNROLL_MN
  52. SGEMM_DEFAULT_UNROLL_MN,
  53. #else
  54. MAX(SGEMM_DEFAULT_UNROLL_M, SGEMM_DEFAULT_UNROLL_N),
  55. #endif
  56. #ifdef HAVE_EXCLUSIVE_CACHE
  57. 1,
  58. #else
  59. 0,
  60. #endif
  61. samax_kTS, samin_kTS, smax_kTS, smin_kTS,
  62. isamax_kTS, isamin_kTS, ismax_kTS, ismin_kTS,
  63. snrm2_kTS, sasum_kTS, scopy_kTS, sdot_kTS,
  64. dsdot_kTS,
  65. srot_kTS, saxpy_kTS, sscal_kTS, sswap_kTS,
  66. sgemv_nTS, sgemv_tTS, sger_kTS,
  67. ssymv_LTS, ssymv_UTS,
  68. sgemm_kernelTS, sgemm_betaTS,
  69. #if SGEMM_DEFAULT_UNROLL_M != SGEMM_DEFAULT_UNROLL_N
  70. sgemm_incopyTS, sgemm_itcopyTS,
  71. #else
  72. sgemm_oncopyTS, sgemm_otcopyTS,
  73. #endif
  74. sgemm_oncopyTS, sgemm_otcopyTS,
  75. strsm_kernel_LNTS, strsm_kernel_LTTS, strsm_kernel_RNTS, strsm_kernel_RTTS,
  76. #if SGEMM_DEFAULT_UNROLL_M != SGEMM_DEFAULT_UNROLL_N
  77. strsm_iunucopyTS, strsm_iunncopyTS, strsm_iutucopyTS, strsm_iutncopyTS,
  78. strsm_ilnucopyTS, strsm_ilnncopyTS, strsm_iltucopyTS, strsm_iltncopyTS,
  79. #else
  80. strsm_ounucopyTS, strsm_ounncopyTS, strsm_outucopyTS, strsm_outncopyTS,
  81. strsm_olnucopyTS, strsm_olnncopyTS, strsm_oltucopyTS, strsm_oltncopyTS,
  82. #endif
  83. strsm_ounucopyTS, strsm_ounncopyTS, strsm_outucopyTS, strsm_outncopyTS,
  84. strsm_olnucopyTS, strsm_olnncopyTS, strsm_oltucopyTS, strsm_oltncopyTS,
  85. strmm_kernel_RNTS, strmm_kernel_RTTS, strmm_kernel_LNTS, strmm_kernel_LTTS,
  86. #if SGEMM_DEFAULT_UNROLL_M != SGEMM_DEFAULT_UNROLL_N
  87. strmm_iunucopyTS, strmm_iunncopyTS, strmm_iutucopyTS, strmm_iutncopyTS,
  88. strmm_ilnucopyTS, strmm_ilnncopyTS, strmm_iltucopyTS, strmm_iltncopyTS,
  89. #else
  90. strmm_ounucopyTS, strmm_ounncopyTS, strmm_outucopyTS, strmm_outncopyTS,
  91. strmm_olnucopyTS, strmm_olnncopyTS, strmm_oltucopyTS, strmm_oltncopyTS,
  92. #endif
  93. strmm_ounucopyTS, strmm_ounncopyTS, strmm_outucopyTS, strmm_outncopyTS,
  94. strmm_olnucopyTS, strmm_olnncopyTS, strmm_oltucopyTS, strmm_oltncopyTS,
  95. #if SGEMM_DEFAULT_UNROLL_M != SGEMM_DEFAULT_UNROLL_N
  96. ssymm_iutcopyTS, ssymm_iltcopyTS,
  97. #else
  98. ssymm_outcopyTS, ssymm_oltcopyTS,
  99. #endif
  100. ssymm_outcopyTS, ssymm_oltcopyTS,
  101. #ifndef NO_LAPACK
  102. sneg_tcopyTS, slaswp_ncopyTS,
  103. #else
  104. NULL,NULL,
  105. #endif
  106. 0, 0, 0,
  107. DGEMM_DEFAULT_UNROLL_M, DGEMM_DEFAULT_UNROLL_N,
  108. #ifdef DGEMM_DEFAULT_UNROLL_MN
  109. DGEMM_DEFAULT_UNROLL_MN,
  110. #else
  111. MAX(DGEMM_DEFAULT_UNROLL_M, DGEMM_DEFAULT_UNROLL_N),
  112. #endif
  113. damax_kTS, damin_kTS, dmax_kTS, dmin_kTS,
  114. idamax_kTS, idamin_kTS, idmax_kTS, idmin_kTS,
  115. dnrm2_kTS, dasum_kTS, dcopy_kTS, ddot_kTS,
  116. drot_kTS, daxpy_kTS, dscal_kTS, dswap_kTS,
  117. dgemv_nTS, dgemv_tTS, dger_kTS,
  118. dsymv_LTS, dsymv_UTS,
  119. dgemm_kernelTS, dgemm_betaTS,
  120. #if DGEMM_DEFAULT_UNROLL_M != DGEMM_DEFAULT_UNROLL_N
  121. dgemm_incopyTS, dgemm_itcopyTS,
  122. #else
  123. dgemm_oncopyTS, dgemm_otcopyTS,
  124. #endif
  125. dgemm_oncopyTS, dgemm_otcopyTS,
  126. dtrsm_kernel_LNTS, dtrsm_kernel_LTTS, dtrsm_kernel_RNTS, dtrsm_kernel_RTTS,
  127. #if DGEMM_DEFAULT_UNROLL_M != DGEMM_DEFAULT_UNROLL_N
  128. dtrsm_iunucopyTS, dtrsm_iunncopyTS, dtrsm_iutucopyTS, dtrsm_iutncopyTS,
  129. dtrsm_ilnucopyTS, dtrsm_ilnncopyTS, dtrsm_iltucopyTS, dtrsm_iltncopyTS,
  130. #else
  131. dtrsm_ounucopyTS, dtrsm_ounncopyTS, dtrsm_outucopyTS, dtrsm_outncopyTS,
  132. dtrsm_olnucopyTS, dtrsm_olnncopyTS, dtrsm_oltucopyTS, dtrsm_oltncopyTS,
  133. #endif
  134. dtrsm_ounucopyTS, dtrsm_ounncopyTS, dtrsm_outucopyTS, dtrsm_outncopyTS,
  135. dtrsm_olnucopyTS, dtrsm_olnncopyTS, dtrsm_oltucopyTS, dtrsm_oltncopyTS,
  136. dtrmm_kernel_RNTS, dtrmm_kernel_RTTS, dtrmm_kernel_LNTS, dtrmm_kernel_LTTS,
  137. #if DGEMM_DEFAULT_UNROLL_M != DGEMM_DEFAULT_UNROLL_N
  138. dtrmm_iunucopyTS, dtrmm_iunncopyTS, dtrmm_iutucopyTS, dtrmm_iutncopyTS,
  139. dtrmm_ilnucopyTS, dtrmm_ilnncopyTS, dtrmm_iltucopyTS, dtrmm_iltncopyTS,
  140. #else
  141. dtrmm_ounucopyTS, dtrmm_ounncopyTS, dtrmm_outucopyTS, dtrmm_outncopyTS,
  142. dtrmm_olnucopyTS, dtrmm_olnncopyTS, dtrmm_oltucopyTS, dtrmm_oltncopyTS,
  143. #endif
  144. dtrmm_ounucopyTS, dtrmm_ounncopyTS, dtrmm_outucopyTS, dtrmm_outncopyTS,
  145. dtrmm_olnucopyTS, dtrmm_olnncopyTS, dtrmm_oltucopyTS, dtrmm_oltncopyTS,
  146. #if DGEMM_DEFAULT_UNROLL_M != DGEMM_DEFAULT_UNROLL_N
  147. dsymm_iutcopyTS, dsymm_iltcopyTS,
  148. #else
  149. dsymm_outcopyTS, dsymm_oltcopyTS,
  150. #endif
  151. dsymm_outcopyTS, dsymm_oltcopyTS,
  152. #ifndef NO_LAPACK
  153. dneg_tcopyTS, dlaswp_ncopyTS,
  154. #else
  155. NULL, NULL,
  156. #endif
  157. #ifdef EXPRECISION
  158. 0, 0, 0,
  159. QGEMM_DEFAULT_UNROLL_M, QGEMM_DEFAULT_UNROLL_N, MAX(QGEMM_DEFAULT_UNROLL_M, QGEMM_DEFAULT_UNROLL_N),
  160. qamax_kTS, qamin_kTS, qmax_kTS, qmin_kTS,
  161. iqamax_kTS, iqamin_kTS, iqmax_kTS, iqmin_kTS,
  162. qnrm2_kTS, qasum_kTS, qcopy_kTS, qdot_kTS,
  163. qrot_kTS, qaxpy_kTS, qscal_kTS, qswap_kTS,
  164. qgemv_nTS, qgemv_tTS, qger_kTS,
  165. qsymv_LTS, qsymv_UTS,
  166. qgemm_kernelTS, qgemm_betaTS,
  167. #if QGEMM_DEFAULT_UNROLL_M != QGEMM_DEFAULT_UNROLL_N
  168. qgemm_incopyTS, qgemm_itcopyTS,
  169. #else
  170. qgemm_oncopyTS, qgemm_otcopyTS,
  171. #endif
  172. qgemm_oncopyTS, qgemm_otcopyTS,
  173. qtrsm_kernel_LNTS, qtrsm_kernel_LTTS, qtrsm_kernel_RNTS, qtrsm_kernel_RTTS,
  174. #if QGEMM_DEFAULT_UNROLL_M != QGEMM_DEFAULT_UNROLL_N
  175. qtrsm_iunucopyTS, qtrsm_iunncopyTS, qtrsm_iutucopyTS, qtrsm_iutncopyTS,
  176. qtrsm_ilnucopyTS, qtrsm_ilnncopyTS, qtrsm_iltucopyTS, qtrsm_iltncopyTS,
  177. #else
  178. qtrsm_ounucopyTS, qtrsm_ounncopyTS, qtrsm_outucopyTS, qtrsm_outncopyTS,
  179. qtrsm_olnucopyTS, qtrsm_olnncopyTS, qtrsm_oltucopyTS, qtrsm_oltncopyTS,
  180. #endif
  181. qtrsm_ounucopyTS, qtrsm_ounncopyTS, qtrsm_outucopyTS, qtrsm_outncopyTS,
  182. qtrsm_olnucopyTS, qtrsm_olnncopyTS, qtrsm_oltucopyTS, qtrsm_oltncopyTS,
  183. qtrmm_kernel_RNTS, qtrmm_kernel_RTTS, qtrmm_kernel_LNTS, qtrmm_kernel_LTTS,
  184. #if QGEMM_DEFAULT_UNROLL_M != QGEMM_DEFAULT_UNROLL_N
  185. qtrmm_iunucopyTS, qtrmm_iunncopyTS, qtrmm_iutucopyTS, qtrmm_iutncopyTS,
  186. qtrmm_ilnucopyTS, qtrmm_ilnncopyTS, qtrmm_iltucopyTS, qtrmm_iltncopyTS,
  187. #else
  188. qtrmm_ounucopyTS, qtrmm_ounncopyTS, qtrmm_outucopyTS, qtrmm_outncopyTS,
  189. qtrmm_olnucopyTS, qtrmm_olnncopyTS, qtrmm_oltucopyTS, qtrmm_oltncopyTS,
  190. #endif
  191. qtrmm_ounucopyTS, qtrmm_ounncopyTS, qtrmm_outucopyTS, qtrmm_outncopyTS,
  192. qtrmm_olnucopyTS, qtrmm_olnncopyTS, qtrmm_oltucopyTS, qtrmm_oltncopyTS,
  193. #if QGEMM_DEFAULT_UNROLL_M != QGEMM_DEFAULT_UNROLL_N
  194. qsymm_iutcopyTS, qsymm_iltcopyTS,
  195. #else
  196. qsymm_outcopyTS, qsymm_oltcopyTS,
  197. #endif
  198. qsymm_outcopyTS, qsymm_oltcopyTS,
  199. #ifndef NO_LAPACK
  200. qneg_tcopyTS, qlaswp_ncopyTS,
  201. #else
  202. NULL, NULL,
  203. #endif
  204. #endif
  205. 0, 0, 0,
  206. CGEMM_DEFAULT_UNROLL_M, CGEMM_DEFAULT_UNROLL_N,
  207. #ifdef CGEMM_DEFAULT_UNROLL_MN
  208. CGEMM_DEFAULT_UNROLL_MN,
  209. #else
  210. MAX(CGEMM_DEFAULT_UNROLL_M, CGEMM_DEFAULT_UNROLL_N),
  211. #endif
  212. camax_kTS, camin_kTS, icamax_kTS, icamin_kTS,
  213. cnrm2_kTS, casum_kTS, ccopy_kTS,
  214. cdotu_kTS, cdotc_kTS, csrot_kTS,
  215. caxpy_kTS, caxpyc_kTS, cscal_kTS, cswap_kTS,
  216. cgemv_nTS, cgemv_tTS, cgemv_rTS, cgemv_cTS,
  217. cgemv_oTS, cgemv_uTS, cgemv_sTS, cgemv_dTS,
  218. cgeru_kTS, cgerc_kTS, cgerv_kTS, cgerd_kTS,
  219. csymv_LTS, csymv_UTS,
  220. chemv_LTS, chemv_UTS, chemv_MTS, chemv_VTS,
  221. cgemm_kernel_nTS, cgemm_kernel_lTS, cgemm_kernel_rTS, cgemm_kernel_bTS,
  222. cgemm_betaTS,
  223. #if CGEMM_DEFAULT_UNROLL_M != CGEMM_DEFAULT_UNROLL_N
  224. cgemm_incopyTS, cgemm_itcopyTS,
  225. #else
  226. cgemm_oncopyTS, cgemm_otcopyTS,
  227. #endif
  228. cgemm_oncopyTS, cgemm_otcopyTS,
  229. ctrsm_kernel_LNTS, ctrsm_kernel_LTTS, ctrsm_kernel_LRTS, ctrsm_kernel_LCTS,
  230. ctrsm_kernel_RNTS, ctrsm_kernel_RTTS, ctrsm_kernel_RRTS, ctrsm_kernel_RCTS,
  231. #if CGEMM_DEFAULT_UNROLL_M != CGEMM_DEFAULT_UNROLL_N
  232. ctrsm_iunucopyTS, ctrsm_iunncopyTS, ctrsm_iutucopyTS, ctrsm_iutncopyTS,
  233. ctrsm_ilnucopyTS, ctrsm_ilnncopyTS, ctrsm_iltucopyTS, ctrsm_iltncopyTS,
  234. #else
  235. ctrsm_ounucopyTS, ctrsm_ounncopyTS, ctrsm_outucopyTS, ctrsm_outncopyTS,
  236. ctrsm_olnucopyTS, ctrsm_olnncopyTS, ctrsm_oltucopyTS, ctrsm_oltncopyTS,
  237. #endif
  238. ctrsm_ounucopyTS, ctrsm_ounncopyTS, ctrsm_outucopyTS, ctrsm_outncopyTS,
  239. ctrsm_olnucopyTS, ctrsm_olnncopyTS, ctrsm_oltucopyTS, ctrsm_oltncopyTS,
  240. ctrmm_kernel_RNTS, ctrmm_kernel_RTTS, ctrmm_kernel_RRTS, ctrmm_kernel_RCTS,
  241. ctrmm_kernel_LNTS, ctrmm_kernel_LTTS, ctrmm_kernel_LRTS, ctrmm_kernel_LCTS,
  242. #if CGEMM_DEFAULT_UNROLL_M != CGEMM_DEFAULT_UNROLL_N
  243. ctrmm_iunucopyTS, ctrmm_iunncopyTS, ctrmm_iutucopyTS, ctrmm_iutncopyTS,
  244. ctrmm_ilnucopyTS, ctrmm_ilnncopyTS, ctrmm_iltucopyTS, ctrmm_iltncopyTS,
  245. #else
  246. ctrmm_ounucopyTS, ctrmm_ounncopyTS, ctrmm_outucopyTS, ctrmm_outncopyTS,
  247. ctrmm_olnucopyTS, ctrmm_olnncopyTS, ctrmm_oltucopyTS, ctrmm_oltncopyTS,
  248. #endif
  249. ctrmm_ounucopyTS, ctrmm_ounncopyTS, ctrmm_outucopyTS, ctrmm_outncopyTS,
  250. ctrmm_olnucopyTS, ctrmm_olnncopyTS, ctrmm_oltucopyTS, ctrmm_oltncopyTS,
  251. #if CGEMM_DEFAULT_UNROLL_M != CGEMM_DEFAULT_UNROLL_N
  252. csymm_iutcopyTS, csymm_iltcopyTS,
  253. #else
  254. csymm_outcopyTS, csymm_oltcopyTS,
  255. #endif
  256. csymm_outcopyTS, csymm_oltcopyTS,
  257. #if CGEMM_DEFAULT_UNROLL_M != CGEMM_DEFAULT_UNROLL_N
  258. chemm_iutcopyTS, chemm_iltcopyTS,
  259. #else
  260. chemm_outcopyTS, chemm_oltcopyTS,
  261. #endif
  262. chemm_outcopyTS, chemm_oltcopyTS,
  263. 0, 0, 0,
  264. #if defined(USE_GEMM3M)
  265. #ifdef CGEMM3M_DEFAULT_UNROLL_M
  266. CGEMM3M_DEFAULT_UNROLL_M, CGEMM3M_DEFAULT_UNROLL_N, MAX(CGEMM3M_DEFAULT_UNROLL_M, CGEMM3M_DEFAULT_UNROLL_N),
  267. #else
  268. SGEMM_DEFAULT_UNROLL_M, SGEMM_DEFAULT_UNROLL_N, MAX(SGEMM_DEFAULT_UNROLL_M, SGEMM_DEFAULT_UNROLL_N),
  269. #endif
  270. cgemm3m_kernelTS,
  271. cgemm3m_incopybTS, cgemm3m_incopyrTS,
  272. cgemm3m_incopyiTS, cgemm3m_itcopybTS,
  273. cgemm3m_itcopyrTS, cgemm3m_itcopyiTS,
  274. cgemm3m_oncopybTS, cgemm3m_oncopyrTS,
  275. cgemm3m_oncopyiTS, cgemm3m_otcopybTS,
  276. cgemm3m_otcopyrTS, cgemm3m_otcopyiTS,
  277. csymm3m_iucopybTS, csymm3m_ilcopybTS,
  278. csymm3m_iucopyrTS, csymm3m_ilcopyrTS,
  279. csymm3m_iucopyiTS, csymm3m_ilcopyiTS,
  280. csymm3m_oucopybTS, csymm3m_olcopybTS,
  281. csymm3m_oucopyrTS, csymm3m_olcopyrTS,
  282. csymm3m_oucopyiTS, csymm3m_olcopyiTS,
  283. chemm3m_iucopybTS, chemm3m_ilcopybTS,
  284. chemm3m_iucopyrTS, chemm3m_ilcopyrTS,
  285. chemm3m_iucopyiTS, chemm3m_ilcopyiTS,
  286. chemm3m_oucopybTS, chemm3m_olcopybTS,
  287. chemm3m_oucopyrTS, chemm3m_olcopyrTS,
  288. chemm3m_oucopyiTS, chemm3m_olcopyiTS,
  289. #else
  290. 0, 0, 0,
  291. NULL,
  292. NULL, NULL,
  293. NULL, NULL,
  294. NULL, NULL,
  295. NULL, NULL,
  296. NULL, NULL,
  297. NULL, NULL,
  298. NULL, NULL,
  299. NULL, NULL,
  300. NULL, NULL,
  301. NULL, NULL,
  302. NULL, NULL,
  303. NULL, NULL,
  304. NULL, NULL,
  305. NULL, NULL,
  306. NULL, NULL,
  307. NULL, NULL,
  308. NULL, NULL,
  309. NULL, NULL,
  310. #endif
  311. #ifndef NO_LAPACK
  312. cneg_tcopyTS, claswp_ncopyTS,
  313. #else
  314. NULL, NULL,
  315. #endif
  316. 0, 0, 0,
  317. ZGEMM_DEFAULT_UNROLL_M, ZGEMM_DEFAULT_UNROLL_N,
  318. #ifdef ZGEMM_DEFAULT_UNROLL_MN
  319. ZGEMM_DEFAULT_UNROLL_MN,
  320. #else
  321. MAX(ZGEMM_DEFAULT_UNROLL_M, ZGEMM_DEFAULT_UNROLL_N),
  322. #endif
  323. zamax_kTS, zamin_kTS, izamax_kTS, izamin_kTS,
  324. znrm2_kTS, zasum_kTS, zcopy_kTS,
  325. zdotu_kTS, zdotc_kTS, zdrot_kTS,
  326. zaxpy_kTS, zaxpyc_kTS, zscal_kTS, zswap_kTS,
  327. zgemv_nTS, zgemv_tTS, zgemv_rTS, zgemv_cTS,
  328. zgemv_oTS, zgemv_uTS, zgemv_sTS, zgemv_dTS,
  329. zgeru_kTS, zgerc_kTS, zgerv_kTS, zgerd_kTS,
  330. zsymv_LTS, zsymv_UTS,
  331. zhemv_LTS, zhemv_UTS, zhemv_MTS, zhemv_VTS,
  332. zgemm_kernel_nTS, zgemm_kernel_lTS, zgemm_kernel_rTS, zgemm_kernel_bTS,
  333. zgemm_betaTS,
  334. #if ZGEMM_DEFAULT_UNROLL_M != ZGEMM_DEFAULT_UNROLL_N
  335. zgemm_incopyTS, zgemm_itcopyTS,
  336. #else
  337. zgemm_oncopyTS, zgemm_otcopyTS,
  338. #endif
  339. zgemm_oncopyTS, zgemm_otcopyTS,
  340. ztrsm_kernel_LNTS, ztrsm_kernel_LTTS, ztrsm_kernel_LRTS, ztrsm_kernel_LCTS,
  341. ztrsm_kernel_RNTS, ztrsm_kernel_RTTS, ztrsm_kernel_RRTS, ztrsm_kernel_RCTS,
  342. #if ZGEMM_DEFAULT_UNROLL_M != ZGEMM_DEFAULT_UNROLL_N
  343. ztrsm_iunucopyTS, ztrsm_iunncopyTS, ztrsm_iutucopyTS, ztrsm_iutncopyTS,
  344. ztrsm_ilnucopyTS, ztrsm_ilnncopyTS, ztrsm_iltucopyTS, ztrsm_iltncopyTS,
  345. #else
  346. ztrsm_ounucopyTS, ztrsm_ounncopyTS, ztrsm_outucopyTS, ztrsm_outncopyTS,
  347. ztrsm_olnucopyTS, ztrsm_olnncopyTS, ztrsm_oltucopyTS, ztrsm_oltncopyTS,
  348. #endif
  349. ztrsm_ounucopyTS, ztrsm_ounncopyTS, ztrsm_outucopyTS, ztrsm_outncopyTS,
  350. ztrsm_olnucopyTS, ztrsm_olnncopyTS, ztrsm_oltucopyTS, ztrsm_oltncopyTS,
  351. ztrmm_kernel_RNTS, ztrmm_kernel_RTTS, ztrmm_kernel_RRTS, ztrmm_kernel_RCTS,
  352. ztrmm_kernel_LNTS, ztrmm_kernel_LTTS, ztrmm_kernel_LRTS, ztrmm_kernel_LCTS,
  353. #if ZGEMM_DEFAULT_UNROLL_M != ZGEMM_DEFAULT_UNROLL_N
  354. ztrmm_iunucopyTS, ztrmm_iunncopyTS, ztrmm_iutucopyTS, ztrmm_iutncopyTS,
  355. ztrmm_ilnucopyTS, ztrmm_ilnncopyTS, ztrmm_iltucopyTS, ztrmm_iltncopyTS,
  356. #else
  357. ztrmm_ounucopyTS, ztrmm_ounncopyTS, ztrmm_outucopyTS, ztrmm_outncopyTS,
  358. ztrmm_olnucopyTS, ztrmm_olnncopyTS, ztrmm_oltucopyTS, ztrmm_oltncopyTS,
  359. #endif
  360. ztrmm_ounucopyTS, ztrmm_ounncopyTS, ztrmm_outucopyTS, ztrmm_outncopyTS,
  361. ztrmm_olnucopyTS, ztrmm_olnncopyTS, ztrmm_oltucopyTS, ztrmm_oltncopyTS,
  362. #if ZGEMM_DEFAULT_UNROLL_M != ZGEMM_DEFAULT_UNROLL_N
  363. zsymm_iutcopyTS, zsymm_iltcopyTS,
  364. #else
  365. zsymm_outcopyTS, zsymm_oltcopyTS,
  366. #endif
  367. zsymm_outcopyTS, zsymm_oltcopyTS,
  368. #if ZGEMM_DEFAULT_UNROLL_M != ZGEMM_DEFAULT_UNROLL_N
  369. zhemm_iutcopyTS, zhemm_iltcopyTS,
  370. #else
  371. zhemm_outcopyTS, zhemm_oltcopyTS,
  372. #endif
  373. zhemm_outcopyTS, zhemm_oltcopyTS,
  374. 0, 0, 0,
  375. #if defined(USE_GEMM3M)
  376. #ifdef ZGEMM3M_DEFAULT_UNROLL_M
  377. ZGEMM3M_DEFAULT_UNROLL_M, ZGEMM3M_DEFAULT_UNROLL_N, MAX(ZGEMM3M_DEFAULT_UNROLL_M, ZGEMM3M_DEFAULT_UNROLL_N),
  378. #else
  379. DGEMM_DEFAULT_UNROLL_M, DGEMM_DEFAULT_UNROLL_N, MAX(DGEMM_DEFAULT_UNROLL_M, DGEMM_DEFAULT_UNROLL_N),
  380. #endif
  381. zgemm3m_kernelTS,
  382. zgemm3m_incopybTS, zgemm3m_incopyrTS,
  383. zgemm3m_incopyiTS, zgemm3m_itcopybTS,
  384. zgemm3m_itcopyrTS, zgemm3m_itcopyiTS,
  385. zgemm3m_oncopybTS, zgemm3m_oncopyrTS,
  386. zgemm3m_oncopyiTS, zgemm3m_otcopybTS,
  387. zgemm3m_otcopyrTS, zgemm3m_otcopyiTS,
  388. zsymm3m_iucopybTS, zsymm3m_ilcopybTS,
  389. zsymm3m_iucopyrTS, zsymm3m_ilcopyrTS,
  390. zsymm3m_iucopyiTS, zsymm3m_ilcopyiTS,
  391. zsymm3m_oucopybTS, zsymm3m_olcopybTS,
  392. zsymm3m_oucopyrTS, zsymm3m_olcopyrTS,
  393. zsymm3m_oucopyiTS, zsymm3m_olcopyiTS,
  394. zhemm3m_iucopybTS, zhemm3m_ilcopybTS,
  395. zhemm3m_iucopyrTS, zhemm3m_ilcopyrTS,
  396. zhemm3m_iucopyiTS, zhemm3m_ilcopyiTS,
  397. zhemm3m_oucopybTS, zhemm3m_olcopybTS,
  398. zhemm3m_oucopyrTS, zhemm3m_olcopyrTS,
  399. zhemm3m_oucopyiTS, zhemm3m_olcopyiTS,
  400. #else
  401. 0, 0, 0,
  402. NULL,
  403. NULL, NULL,
  404. NULL, NULL,
  405. NULL, NULL,
  406. NULL, NULL,
  407. NULL, NULL,
  408. NULL, NULL,
  409. NULL, NULL,
  410. NULL, NULL,
  411. NULL, NULL,
  412. NULL, NULL,
  413. NULL, NULL,
  414. NULL, NULL,
  415. NULL, NULL,
  416. NULL, NULL,
  417. NULL, NULL,
  418. NULL, NULL,
  419. NULL, NULL,
  420. NULL, NULL,
  421. #endif
  422. #ifndef NO_LAPACK
  423. zneg_tcopyTS, zlaswp_ncopyTS,
  424. #else
  425. NULL, NULL,
  426. #endif
  427. #ifdef EXPRECISION
  428. 0, 0, 0,
  429. XGEMM_DEFAULT_UNROLL_M, XGEMM_DEFAULT_UNROLL_N, MAX(XGEMM_DEFAULT_UNROLL_M, XGEMM_DEFAULT_UNROLL_N),
  430. xamax_kTS, xamin_kTS, ixamax_kTS, ixamin_kTS,
  431. xnrm2_kTS, xasum_kTS, xcopy_kTS,
  432. xdotu_kTS, xdotc_kTS, xqrot_kTS,
  433. xaxpy_kTS, xaxpyc_kTS, xscal_kTS, xswap_kTS,
  434. xgemv_nTS, xgemv_tTS, xgemv_rTS, xgemv_cTS,
  435. xgemv_oTS, xgemv_uTS, xgemv_sTS, xgemv_dTS,
  436. xgeru_kTS, xgerc_kTS, xgerv_kTS, xgerd_kTS,
  437. xsymv_LTS, xsymv_UTS,
  438. xhemv_LTS, xhemv_UTS, xhemv_MTS, xhemv_VTS,
  439. xgemm_kernel_nTS, xgemm_kernel_lTS, xgemm_kernel_rTS, xgemm_kernel_bTS,
  440. xgemm_betaTS,
  441. #if XGEMM_DEFAULT_UNROLL_M != XGEMM_DEFAULT_UNROLL_N
  442. xgemm_incopyTS, xgemm_itcopyTS,
  443. #else
  444. xgemm_oncopyTS, xgemm_otcopyTS,
  445. #endif
  446. xgemm_oncopyTS, xgemm_otcopyTS,
  447. xtrsm_kernel_LNTS, xtrsm_kernel_LTTS, xtrsm_kernel_LRTS, xtrsm_kernel_LCTS,
  448. xtrsm_kernel_RNTS, xtrsm_kernel_RTTS, xtrsm_kernel_RRTS, xtrsm_kernel_RCTS,
  449. #if XGEMM_DEFAULT_UNROLL_M != XGEMM_DEFAULT_UNROLL_N
  450. xtrsm_iunucopyTS, xtrsm_iunncopyTS, xtrsm_iutucopyTS, xtrsm_iutncopyTS,
  451. xtrsm_ilnucopyTS, xtrsm_ilnncopyTS, xtrsm_iltucopyTS, xtrsm_iltncopyTS,
  452. #else
  453. xtrsm_ounucopyTS, xtrsm_ounncopyTS, xtrsm_outucopyTS, xtrsm_outncopyTS,
  454. xtrsm_olnucopyTS, xtrsm_olnncopyTS, xtrsm_oltucopyTS, xtrsm_oltncopyTS,
  455. #endif
  456. xtrsm_ounucopyTS, xtrsm_ounncopyTS, xtrsm_outucopyTS, xtrsm_outncopyTS,
  457. xtrsm_olnucopyTS, xtrsm_olnncopyTS, xtrsm_oltucopyTS, xtrsm_oltncopyTS,
  458. xtrmm_kernel_RNTS, xtrmm_kernel_RTTS, xtrmm_kernel_RRTS, xtrmm_kernel_RCTS,
  459. xtrmm_kernel_LNTS, xtrmm_kernel_LTTS, xtrmm_kernel_LRTS, xtrmm_kernel_LCTS,
  460. #if XGEMM_DEFAULT_UNROLL_M != XGEMM_DEFAULT_UNROLL_N
  461. xtrmm_iunucopyTS, xtrmm_iunncopyTS, xtrmm_iutucopyTS, xtrmm_iutncopyTS,
  462. xtrmm_ilnucopyTS, xtrmm_ilnncopyTS, xtrmm_iltucopyTS, xtrmm_iltncopyTS,
  463. #else
  464. xtrmm_ounucopyTS, xtrmm_ounncopyTS, xtrmm_outucopyTS, xtrmm_outncopyTS,
  465. xtrmm_olnucopyTS, xtrmm_olnncopyTS, xtrmm_oltucopyTS, xtrmm_oltncopyTS,
  466. #endif
  467. xtrmm_ounucopyTS, xtrmm_ounncopyTS, xtrmm_outucopyTS, xtrmm_outncopyTS,
  468. xtrmm_olnucopyTS, xtrmm_olnncopyTS, xtrmm_oltucopyTS, xtrmm_oltncopyTS,
  469. #if XGEMM_DEFAULT_UNROLL_M != XGEMM_DEFAULT_UNROLL_N
  470. xsymm_iutcopyTS, xsymm_iltcopyTS,
  471. #else
  472. xsymm_outcopyTS, xsymm_oltcopyTS,
  473. #endif
  474. xsymm_outcopyTS, xsymm_oltcopyTS,
  475. #if XGEMM_DEFAULT_UNROLL_M != XGEMM_DEFAULT_UNROLL_N
  476. xhemm_iutcopyTS, xhemm_iltcopyTS,
  477. #else
  478. xhemm_outcopyTS, xhemm_oltcopyTS,
  479. #endif
  480. xhemm_outcopyTS, xhemm_oltcopyTS,
  481. 0, 0, 0,
  482. #if defined(USE_GEMM3M)
  483. QGEMM_DEFAULT_UNROLL_M, QGEMM_DEFAULT_UNROLL_N, MAX(QGEMM_DEFAULT_UNROLL_M, QGEMM_DEFAULT_UNROLL_N),
  484. xgemm3m_kernelTS,
  485. xgemm3m_incopybTS, xgemm3m_incopyrTS,
  486. xgemm3m_incopyiTS, xgemm3m_itcopybTS,
  487. xgemm3m_itcopyrTS, xgemm3m_itcopyiTS,
  488. xgemm3m_oncopybTS, xgemm3m_oncopyrTS,
  489. xgemm3m_oncopyiTS, xgemm3m_otcopybTS,
  490. xgemm3m_otcopyrTS, xgemm3m_otcopyiTS,
  491. xsymm3m_iucopybTS, xsymm3m_ilcopybTS,
  492. xsymm3m_iucopyrTS, xsymm3m_ilcopyrTS,
  493. xsymm3m_iucopyiTS, xsymm3m_ilcopyiTS,
  494. xsymm3m_oucopybTS, xsymm3m_olcopybTS,
  495. xsymm3m_oucopyrTS, xsymm3m_olcopyrTS,
  496. xsymm3m_oucopyiTS, xsymm3m_olcopyiTS,
  497. xhemm3m_iucopybTS, xhemm3m_ilcopybTS,
  498. xhemm3m_iucopyrTS, xhemm3m_ilcopyrTS,
  499. xhemm3m_iucopyiTS, xhemm3m_ilcopyiTS,
  500. xhemm3m_oucopybTS, xhemm3m_olcopybTS,
  501. xhemm3m_oucopyrTS, xhemm3m_olcopyrTS,
  502. xhemm3m_oucopyiTS, xhemm3m_olcopyiTS,
  503. #else
  504. 0, 0, 0,
  505. NULL,
  506. NULL, NULL,
  507. NULL, NULL,
  508. NULL, NULL,
  509. NULL, NULL,
  510. NULL, NULL,
  511. NULL, NULL,
  512. NULL, NULL,
  513. NULL, NULL,
  514. NULL, NULL,
  515. NULL, NULL,
  516. NULL, NULL,
  517. NULL, NULL,
  518. NULL, NULL,
  519. NULL, NULL,
  520. NULL, NULL,
  521. NULL, NULL,
  522. NULL, NULL,
  523. NULL, NULL,
  524. #endif
  525. #ifndef NO_LAPACK
  526. xneg_tcopyTS, xlaswp_ncopyTS,
  527. #else
  528. NULL, NULL,
  529. #endif
  530. #endif
  531. init_parameter,
  532. SNUMOPT, DNUMOPT, QNUMOPT,
  533. saxpby_kTS, daxpby_kTS, caxpby_kTS, zaxpby_kTS,
  534. somatcopy_k_cnTS, somatcopy_k_ctTS, somatcopy_k_rnTS, somatcopy_k_rtTS,
  535. domatcopy_k_cnTS, domatcopy_k_ctTS, domatcopy_k_rnTS, domatcopy_k_rtTS,
  536. comatcopy_k_cnTS, comatcopy_k_ctTS, comatcopy_k_rnTS, comatcopy_k_rtTS,
  537. comatcopy_k_cncTS, comatcopy_k_ctcTS, comatcopy_k_rncTS, comatcopy_k_rtcTS,
  538. zomatcopy_k_cnTS, zomatcopy_k_ctTS, zomatcopy_k_rnTS, zomatcopy_k_rtTS,
  539. zomatcopy_k_cncTS, zomatcopy_k_ctcTS, zomatcopy_k_rncTS, zomatcopy_k_rtcTS,
  540. simatcopy_k_cnTS, simatcopy_k_ctTS, simatcopy_k_rnTS, simatcopy_k_rtTS,
  541. dimatcopy_k_cnTS, dimatcopy_k_ctTS, dimatcopy_k_rnTS, dimatcopy_k_rtTS,
  542. cimatcopy_k_cnTS, cimatcopy_k_ctTS, cimatcopy_k_rnTS, cimatcopy_k_rtTS,
  543. cimatcopy_k_cncTS, cimatcopy_k_ctcTS, cimatcopy_k_rncTS, cimatcopy_k_rtcTS,
  544. zimatcopy_k_cnTS, zimatcopy_k_ctTS, zimatcopy_k_rnTS, zimatcopy_k_rtTS,
  545. zimatcopy_k_cncTS, zimatcopy_k_ctcTS, zimatcopy_k_rncTS, zimatcopy_k_rtcTS,
  546. sgeadd_kTS, dgeadd_kTS, cgeadd_kTS, zgeadd_kTS
  547. };
  548. #ifdef ARCH_X86
  549. static int get_l2_size_old(void){
  550. int i, eax, ebx, ecx, edx, cpuid_level;
  551. int info[15];
  552. cpuid(2, &eax, &ebx, &ecx, &edx);
  553. info[ 0] = BITMASK(eax, 8, 0xff);
  554. info[ 1] = BITMASK(eax, 16, 0xff);
  555. info[ 2] = BITMASK(eax, 24, 0xff);
  556. info[ 3] = BITMASK(ebx, 0, 0xff);
  557. info[ 4] = BITMASK(ebx, 8, 0xff);
  558. info[ 5] = BITMASK(ebx, 16, 0xff);
  559. info[ 6] = BITMASK(ebx, 24, 0xff);
  560. info[ 7] = BITMASK(ecx, 0, 0xff);
  561. info[ 8] = BITMASK(ecx, 8, 0xff);
  562. info[ 9] = BITMASK(ecx, 16, 0xff);
  563. info[10] = BITMASK(ecx, 24, 0xff);
  564. info[11] = BITMASK(edx, 0, 0xff);
  565. info[12] = BITMASK(edx, 8, 0xff);
  566. info[13] = BITMASK(edx, 16, 0xff);
  567. info[14] = BITMASK(edx, 24, 0xff);
  568. for (i = 0; i < 15; i++){
  569. switch (info[i]){
  570. /* This table is from http://www.sandpile.org/ia32/cpuid.htm */
  571. case 0x1a :
  572. return 96;
  573. case 0x39 :
  574. case 0x3b :
  575. case 0x41 :
  576. case 0x79 :
  577. case 0x81 :
  578. return 128;
  579. case 0x3a :
  580. return 192;
  581. case 0x21 :
  582. case 0x3c :
  583. case 0x42 :
  584. case 0x7a :
  585. case 0x7e :
  586. case 0x82 :
  587. return 256;
  588. case 0x3d :
  589. return 384;
  590. case 0x3e :
  591. case 0x43 :
  592. case 0x7b :
  593. case 0x7f :
  594. case 0x83 :
  595. case 0x86 :
  596. return 512;
  597. case 0x44 :
  598. case 0x78 :
  599. case 0x7c :
  600. case 0x84 :
  601. case 0x87 :
  602. return 1024;
  603. case 0x45 :
  604. case 0x7d :
  605. case 0x85 :
  606. return 2048;
  607. case 0x48 :
  608. return 3184;
  609. case 0x49 :
  610. return 4096;
  611. case 0x4e :
  612. return 6144;
  613. }
  614. }
  615. // return 0;
  616. fprintf (stderr,"OpenBLAS WARNING - could not determine the L2 cache size on this system, assuming 256k\n");
  617. return 256;
  618. }
  619. #endif
  620. static __inline__ int get_l2_size(void){
  621. int eax, ebx, ecx, edx, l2;
  622. cpuid(0x80000006, &eax, &ebx, &ecx, &edx);
  623. l2 = BITMASK(ecx, 16, 0xffff);
  624. #ifndef ARCH_X86
  625. if (l2 <= 0) {
  626. fprintf (stderr,"OpenBLAS WARNING - could not determine the L2 cache size on this system, assuming 256k\n");
  627. return 256;
  628. }
  629. return l2;
  630. #else
  631. if (l2 > 0) return l2;
  632. return get_l2_size_old();
  633. #endif
  634. }
  635. static __inline__ int get_l3_size(void){
  636. int eax, ebx, ecx, edx;
  637. cpuid(0x80000006, &eax, &ebx, &ecx, &edx);
  638. return BITMASK(edx, 18, 0x3fff) * 512;
  639. }
  640. static void init_parameter(void) {
  641. int l2 = get_l2_size();
  642. (void) l2; /* dirty trick to suppress unused variable warning for targets */
  643. /* where the GEMM unrolling parameters do not depend on l2 */
  644. TABLE_NAME.sgemm_q = SGEMM_DEFAULT_Q;
  645. TABLE_NAME.dgemm_q = DGEMM_DEFAULT_Q;
  646. TABLE_NAME.cgemm_q = CGEMM_DEFAULT_Q;
  647. TABLE_NAME.zgemm_q = ZGEMM_DEFAULT_Q;
  648. #ifdef CGEMM3M_DEFAULT_Q
  649. TABLE_NAME.cgemm3m_q = CGEMM3M_DEFAULT_Q;
  650. #else
  651. TABLE_NAME.cgemm3m_q = SGEMM_DEFAULT_Q;
  652. #endif
  653. #ifdef ZGEMM3M_DEFAULT_Q
  654. TABLE_NAME.zgemm3m_q = ZGEMM3M_DEFAULT_Q;
  655. #else
  656. TABLE_NAME.zgemm3m_q = DGEMM_DEFAULT_Q;
  657. #endif
  658. #ifdef EXPRECISION
  659. TABLE_NAME.qgemm_q = QGEMM_DEFAULT_Q;
  660. TABLE_NAME.xgemm_q = XGEMM_DEFAULT_Q;
  661. TABLE_NAME.xgemm3m_q = QGEMM_DEFAULT_Q;
  662. #endif
  663. #if defined(CORE_KATMAI) || defined(CORE_COPPERMINE) || defined(CORE_BANIAS) || defined(CORE_YONAH) || defined(CORE_ATHLON)
  664. #ifdef DEBUG
  665. fprintf(stderr, "Katmai, Coppermine, Banias, Athlon\n");
  666. #endif
  667. TABLE_NAME.sgemm_p = 64 * (l2 >> 7);
  668. TABLE_NAME.dgemm_p = 32 * (l2 >> 7);
  669. TABLE_NAME.cgemm_p = 32 * (l2 >> 7);
  670. TABLE_NAME.zgemm_p = 16 * (l2 >> 7);
  671. #ifdef EXPRECISION
  672. TABLE_NAME.qgemm_p = 16 * (l2 >> 7);
  673. TABLE_NAME.xgemm_p = 8 * (l2 >> 7);
  674. #endif
  675. #endif
  676. #ifdef CORE_NORTHWOOD
  677. #ifdef DEBUG
  678. fprintf(stderr, "Northwood\n");
  679. #endif
  680. TABLE_NAME.sgemm_p = 96 * (l2 >> 7);
  681. TABLE_NAME.dgemm_p = 48 * (l2 >> 7);
  682. TABLE_NAME.cgemm_p = 48 * (l2 >> 7);
  683. TABLE_NAME.zgemm_p = 24 * (l2 >> 7);
  684. #ifdef EXPRECISION
  685. TABLE_NAME.qgemm_p = 24 * (l2 >> 7);
  686. TABLE_NAME.xgemm_p = 12 * (l2 >> 7);
  687. #endif
  688. #endif
  689. #ifdef ATOM
  690. #ifdef DEBUG
  691. fprintf(stderr, "Atom\n");
  692. #endif
  693. TABLE_NAME.sgemm_p = 256;
  694. TABLE_NAME.dgemm_p = 128;
  695. TABLE_NAME.cgemm_p = 128;
  696. TABLE_NAME.zgemm_p = 64;
  697. #ifdef EXPRECISION
  698. TABLE_NAME.qgemm_p = 64;
  699. TABLE_NAME.xgemm_p = 32;
  700. #endif
  701. #endif
  702. #ifdef CORE_PRESCOTT
  703. #ifdef DEBUG
  704. fprintf(stderr, "Prescott\n");
  705. #endif
  706. TABLE_NAME.sgemm_p = 56 * (l2 >> 7);
  707. TABLE_NAME.dgemm_p = 28 * (l2 >> 7);
  708. TABLE_NAME.cgemm_p = 28 * (l2 >> 7);
  709. TABLE_NAME.zgemm_p = 14 * (l2 >> 7);
  710. #ifdef EXPRECISION
  711. TABLE_NAME.qgemm_p = 14 * (l2 >> 7);
  712. TABLE_NAME.xgemm_p = 7 * (l2 >> 7);
  713. #endif
  714. #endif
  715. #ifdef CORE2
  716. #ifdef DEBUG
  717. fprintf(stderr, "Core2\n");
  718. #endif
  719. TABLE_NAME.sgemm_p = 92 * (l2 >> 9) + 8;
  720. TABLE_NAME.dgemm_p = 46 * (l2 >> 9) + 8;
  721. TABLE_NAME.cgemm_p = 46 * (l2 >> 9) + 4;
  722. TABLE_NAME.zgemm_p = 23 * (l2 >> 9) + 4;
  723. #ifdef EXPRECISION
  724. TABLE_NAME.qgemm_p = 92 * (l2 >> 9) + 8;
  725. TABLE_NAME.xgemm_p = 46 * (l2 >> 9) + 4;
  726. #endif
  727. #endif
  728. #ifdef PENRYN
  729. #ifdef DEBUG
  730. fprintf(stderr, "Penryn\n");
  731. #endif
  732. TABLE_NAME.sgemm_p = 42 * (l2 >> 9) + 8;
  733. TABLE_NAME.dgemm_p = 42 * (l2 >> 9) + 8;
  734. TABLE_NAME.cgemm_p = 21 * (l2 >> 9) + 4;
  735. TABLE_NAME.zgemm_p = 21 * (l2 >> 9) + 4;
  736. #ifdef EXPRECISION
  737. TABLE_NAME.qgemm_p = 42 * (l2 >> 9) + 8;
  738. TABLE_NAME.xgemm_p = 21 * (l2 >> 9) + 4;
  739. #endif
  740. #endif
  741. #ifdef DUNNINGTON
  742. #ifdef DEBUG
  743. fprintf(stderr, "Dunnington\n");
  744. #endif
  745. TABLE_NAME.sgemm_p = 42 * (l2 >> 9) + 8;
  746. TABLE_NAME.dgemm_p = 42 * (l2 >> 9) + 8;
  747. TABLE_NAME.cgemm_p = 21 * (l2 >> 9) + 4;
  748. TABLE_NAME.zgemm_p = 21 * (l2 >> 9) + 4;
  749. #ifdef EXPRECISION
  750. TABLE_NAME.qgemm_p = 42 * (l2 >> 9) + 8;
  751. TABLE_NAME.xgemm_p = 21 * (l2 >> 9) + 4;
  752. #endif
  753. #endif
  754. #ifdef NEHALEM
  755. #ifdef DEBUG
  756. fprintf(stderr, "Nehalem\n");
  757. #endif
  758. TABLE_NAME.sgemm_p = SGEMM_DEFAULT_P;
  759. TABLE_NAME.dgemm_p = DGEMM_DEFAULT_P;
  760. TABLE_NAME.cgemm_p = CGEMM_DEFAULT_P;
  761. TABLE_NAME.zgemm_p = ZGEMM_DEFAULT_P;
  762. #ifdef EXPRECISION
  763. TABLE_NAME.qgemm_p = QGEMM_DEFAULT_P;
  764. TABLE_NAME.xgemm_p = XGEMM_DEFAULT_P;
  765. #endif
  766. #endif
  767. #ifdef SANDYBRIDGE
  768. #ifdef DEBUG
  769. fprintf(stderr, "Sandybridge\n");
  770. #endif
  771. TABLE_NAME.sgemm_p = SGEMM_DEFAULT_P;
  772. TABLE_NAME.dgemm_p = DGEMM_DEFAULT_P;
  773. TABLE_NAME.cgemm_p = CGEMM_DEFAULT_P;
  774. TABLE_NAME.zgemm_p = ZGEMM_DEFAULT_P;
  775. #ifdef EXPRECISION
  776. TABLE_NAME.qgemm_p = QGEMM_DEFAULT_P;
  777. TABLE_NAME.xgemm_p = XGEMM_DEFAULT_P;
  778. #endif
  779. #endif
  780. #ifdef HASWELL
  781. #ifdef DEBUG
  782. fprintf(stderr, "Haswell\n");
  783. #endif
  784. TABLE_NAME.sgemm_p = SGEMM_DEFAULT_P;
  785. TABLE_NAME.dgemm_p = DGEMM_DEFAULT_P;
  786. TABLE_NAME.cgemm_p = CGEMM_DEFAULT_P;
  787. TABLE_NAME.zgemm_p = ZGEMM_DEFAULT_P;
  788. #ifdef EXPRECISION
  789. TABLE_NAME.qgemm_p = QGEMM_DEFAULT_P;
  790. TABLE_NAME.xgemm_p = XGEMM_DEFAULT_P;
  791. #endif
  792. #endif
  793. #ifdef SKYLAKEX
  794. #ifdef DEBUG
  795. fprintf(stderr, "SkylakeX\n");
  796. #endif
  797. TABLE_NAME.sgemm_p = SGEMM_DEFAULT_P;
  798. TABLE_NAME.dgemm_p = DGEMM_DEFAULT_P;
  799. TABLE_NAME.cgemm_p = CGEMM_DEFAULT_P;
  800. TABLE_NAME.zgemm_p = ZGEMM_DEFAULT_P;
  801. #ifdef EXPRECISION
  802. TABLE_NAME.qgemm_p = QGEMM_DEFAULT_P;
  803. TABLE_NAME.xgemm_p = XGEMM_DEFAULT_P;
  804. #endif
  805. #endif
  806. #ifdef OPTERON
  807. #ifdef DEBUG
  808. fprintf(stderr, "Opteron\n");
  809. #endif
  810. TABLE_NAME.sgemm_p = 224 + 56 * (l2 >> 7);
  811. TABLE_NAME.dgemm_p = 112 + 28 * (l2 >> 7);
  812. TABLE_NAME.cgemm_p = 112 + 28 * (l2 >> 7);
  813. TABLE_NAME.zgemm_p = 56 + 14 * (l2 >> 7);
  814. #ifdef EXPRECISION
  815. TABLE_NAME.qgemm_p = 56 + 14 * (l2 >> 7);
  816. TABLE_NAME.xgemm_p = 28 + 7 * (l2 >> 7);
  817. #endif
  818. #endif
  819. #ifdef BARCELONA
  820. #ifdef DEBUG
  821. fprintf(stderr, "Barcelona\n");
  822. #endif
  823. TABLE_NAME.sgemm_p = SGEMM_DEFAULT_P;
  824. TABLE_NAME.dgemm_p = DGEMM_DEFAULT_P;
  825. TABLE_NAME.cgemm_p = CGEMM_DEFAULT_P;
  826. TABLE_NAME.zgemm_p = ZGEMM_DEFAULT_P;
  827. #ifdef EXPRECISION
  828. TABLE_NAME.qgemm_p = QGEMM_DEFAULT_P;
  829. TABLE_NAME.xgemm_p = XGEMM_DEFAULT_P;
  830. #endif
  831. #endif
  832. #ifdef BOBCAT
  833. #ifdef DEBUG
  834. fprintf(stderr, "Bobcate\n");
  835. #endif
  836. TABLE_NAME.sgemm_p = SGEMM_DEFAULT_P;
  837. TABLE_NAME.dgemm_p = DGEMM_DEFAULT_P;
  838. TABLE_NAME.cgemm_p = CGEMM_DEFAULT_P;
  839. TABLE_NAME.zgemm_p = ZGEMM_DEFAULT_P;
  840. #ifdef EXPRECISION
  841. TABLE_NAME.qgemm_p = QGEMM_DEFAULT_P;
  842. TABLE_NAME.xgemm_p = XGEMM_DEFAULT_P;
  843. #endif
  844. #endif
  845. #ifdef BULLDOZER
  846. #ifdef DEBUG
  847. fprintf(stderr, "Bulldozer\n");
  848. #endif
  849. TABLE_NAME.sgemm_p = SGEMM_DEFAULT_P;
  850. TABLE_NAME.dgemm_p = DGEMM_DEFAULT_P;
  851. TABLE_NAME.cgemm_p = CGEMM_DEFAULT_P;
  852. TABLE_NAME.zgemm_p = ZGEMM_DEFAULT_P;
  853. #ifdef EXPRECISION
  854. TABLE_NAME.qgemm_p = QGEMM_DEFAULT_P;
  855. TABLE_NAME.xgemm_p = XGEMM_DEFAULT_P;
  856. #endif
  857. #endif
  858. #ifdef EXCAVATOR
  859. #ifdef DEBUG
  860. fprintf(stderr, "Excavator\n");
  861. #endif
  862. TABLE_NAME.sgemm_p = SGEMM_DEFAULT_P;
  863. TABLE_NAME.dgemm_p = DGEMM_DEFAULT_P;
  864. TABLE_NAME.cgemm_p = CGEMM_DEFAULT_P;
  865. TABLE_NAME.zgemm_p = ZGEMM_DEFAULT_P;
  866. #ifdef EXPRECISION
  867. TABLE_NAME.qgemm_p = QGEMM_DEFAULT_P;
  868. TABLE_NAME.xgemm_p = XGEMM_DEFAULT_P;
  869. #endif
  870. #endif
  871. #ifdef PILEDRIVER
  872. #ifdef DEBUG
  873. fprintf(stderr, "Piledriver\n");
  874. #endif
  875. TABLE_NAME.sgemm_p = SGEMM_DEFAULT_P;
  876. TABLE_NAME.dgemm_p = DGEMM_DEFAULT_P;
  877. TABLE_NAME.cgemm_p = CGEMM_DEFAULT_P;
  878. TABLE_NAME.zgemm_p = ZGEMM_DEFAULT_P;
  879. #ifdef EXPRECISION
  880. TABLE_NAME.qgemm_p = QGEMM_DEFAULT_P;
  881. TABLE_NAME.xgemm_p = XGEMM_DEFAULT_P;
  882. #endif
  883. #endif
  884. #ifdef STEAMROLLER
  885. #ifdef DEBUG
  886. fprintf(stderr, "Steamroller\n");
  887. #endif
  888. TABLE_NAME.sgemm_p = SGEMM_DEFAULT_P;
  889. TABLE_NAME.dgemm_p = DGEMM_DEFAULT_P;
  890. TABLE_NAME.cgemm_p = CGEMM_DEFAULT_P;
  891. TABLE_NAME.zgemm_p = ZGEMM_DEFAULT_P;
  892. #ifdef EXPRECISION
  893. TABLE_NAME.qgemm_p = QGEMM_DEFAULT_P;
  894. TABLE_NAME.xgemm_p = XGEMM_DEFAULT_P;
  895. #endif
  896. #endif
  897. #ifdef ZEN
  898. #ifdef DEBUG
  899. fprintf(stderr, "Zen\n");
  900. #endif
  901. TABLE_NAME.sgemm_p = SGEMM_DEFAULT_P;
  902. TABLE_NAME.dgemm_p = DGEMM_DEFAULT_P;
  903. TABLE_NAME.cgemm_p = CGEMM_DEFAULT_P;
  904. TABLE_NAME.zgemm_p = ZGEMM_DEFAULT_P;
  905. #ifdef EXPRECISION
  906. TABLE_NAME.qgemm_p = QGEMM_DEFAULT_P;
  907. TABLE_NAME.xgemm_p = XGEMM_DEFAULT_P;
  908. #endif
  909. #endif
  910. #ifdef NANO
  911. #ifdef DEBUG
  912. fprintf(stderr, "NANO\n");
  913. #endif
  914. TABLE_NAME.sgemm_p = SGEMM_DEFAULT_P;
  915. TABLE_NAME.dgemm_p = DGEMM_DEFAULT_P;
  916. TABLE_NAME.cgemm_p = CGEMM_DEFAULT_P;
  917. TABLE_NAME.zgemm_p = ZGEMM_DEFAULT_P;
  918. #ifdef EXPRECISION
  919. TABLE_NAME.qgemm_p = QGEMM_DEFAULT_P;
  920. TABLE_NAME.xgemm_p = XGEMM_DEFAULT_P;
  921. #endif
  922. #endif
  923. #ifdef CGEMM3M_DEFAULT_P
  924. TABLE_NAME.cgemm3m_p = CGEMM3M_DEFAULT_P;
  925. #else
  926. TABLE_NAME.cgemm3m_p = TABLE_NAME.sgemm_p;
  927. #endif
  928. #ifdef ZGEMM3M_DEFAULT_P
  929. TABLE_NAME.zgemm3m_p = ZGEMM3M_DEFAULT_P;
  930. #else
  931. TABLE_NAME.zgemm3m_p = TABLE_NAME.dgemm_p;
  932. #endif
  933. #ifdef EXPRECISION
  934. TABLE_NAME.xgemm3m_p = TABLE_NAME.qgemm_p;
  935. #endif
  936. TABLE_NAME.sgemm_p = ((TABLE_NAME.sgemm_p + SGEMM_DEFAULT_UNROLL_M - 1)/SGEMM_DEFAULT_UNROLL_M) * SGEMM_DEFAULT_UNROLL_M;
  937. TABLE_NAME.dgemm_p = ((TABLE_NAME.dgemm_p + DGEMM_DEFAULT_UNROLL_M - 1)/DGEMM_DEFAULT_UNROLL_M) * DGEMM_DEFAULT_UNROLL_M;
  938. TABLE_NAME.cgemm_p = ((TABLE_NAME.cgemm_p + CGEMM_DEFAULT_UNROLL_M - 1)/CGEMM_DEFAULT_UNROLL_M) * CGEMM_DEFAULT_UNROLL_M;
  939. TABLE_NAME.zgemm_p = ((TABLE_NAME.zgemm_p + ZGEMM_DEFAULT_UNROLL_M - 1)/ZGEMM_DEFAULT_UNROLL_M) * ZGEMM_DEFAULT_UNROLL_M;
  940. #ifdef CGEMM3M_DEFAULT_UNROLL_M
  941. TABLE_NAME.cgemm3m_p = ((TABLE_NAME.cgemm3m_p + CGEMM3M_DEFAULT_UNROLL_M - 1)/CGEMM3M_DEFAULT_UNROLL_M) * CGEMM3M_DEFAULT_UNROLL_M;
  942. #else
  943. TABLE_NAME.cgemm3m_p = ((TABLE_NAME.cgemm3m_p + SGEMM_DEFAULT_UNROLL_M - 1)/SGEMM_DEFAULT_UNROLL_M) * SGEMM_DEFAULT_UNROLL_M;
  944. #endif
  945. #ifdef ZGEMM3M_DEFAULT_UNROLL_M
  946. TABLE_NAME.zgemm3m_p = ((TABLE_NAME.zgemm3m_p + ZGEMM3M_DEFAULT_UNROLL_M - 1)/ZGEMM3M_DEFAULT_UNROLL_M) * ZGEMM3M_DEFAULT_UNROLL_M;
  947. #else
  948. TABLE_NAME.zgemm3m_p = ((TABLE_NAME.zgemm3m_p + DGEMM_DEFAULT_UNROLL_M - 1)/DGEMM_DEFAULT_UNROLL_M) * DGEMM_DEFAULT_UNROLL_M;
  949. #endif
  950. #ifdef QUAD_PRECISION
  951. TABLE_NAME.qgemm_p = ((TABLE_NAME.qgemm_p + QGEMM_DEFAULT_UNROLL_M - 1)/QGEMM_DEFAULT_UNROLL_M) * QGEMM_DEFAULT_UNROLL_M;
  952. TABLE_NAME.xgemm_p = ((TABLE_NAME.xgemm_p + XGEMM_DEFAULT_UNROLL_M - 1)/XGEMM_DEFAULT_UNROLL_M) * XGEMM_DEFAULT_UNROLL_M;
  953. TABLE_NAME.xgemm3m_p = ((TABLE_NAME.xgemm3m_p + QGEMM_DEFAULT_UNROLL_M - 1)/QGEMM_DEFAULT_UNROLL_M) * QGEMM_DEFAULT_UNROLL_M;
  954. #endif
  955. #ifdef DEBUG
  956. fprintf(stderr, "L2 = %8d DGEMM_P .. %d\n", l2, TABLE_NAME.dgemm_p);
  957. #endif
  958. TABLE_NAME.sgemm_r = (((BUFFER_SIZE -
  959. ((TABLE_NAME.sgemm_p * TABLE_NAME.sgemm_q * 4 + TABLE_NAME.offsetA
  960. + TABLE_NAME.align) & ~TABLE_NAME.align)
  961. ) / (TABLE_NAME.sgemm_q * 4) - 15) & ~15);
  962. TABLE_NAME.dgemm_r = (((BUFFER_SIZE -
  963. ((TABLE_NAME.dgemm_p * TABLE_NAME.dgemm_q * 8 + TABLE_NAME.offsetA
  964. + TABLE_NAME.align) & ~TABLE_NAME.align)
  965. ) / (TABLE_NAME.dgemm_q * 8) - 15) & ~15);
  966. #ifdef EXPRECISION
  967. TABLE_NAME.qgemm_r = (((BUFFER_SIZE -
  968. ((TABLE_NAME.qgemm_p * TABLE_NAME.qgemm_q * 16 + TABLE_NAME.offsetA
  969. + TABLE_NAME.align) & ~TABLE_NAME.align)
  970. ) / (TABLE_NAME.qgemm_q * 16) - 15) & ~15);
  971. #endif
  972. TABLE_NAME.cgemm_r = (((BUFFER_SIZE -
  973. ((TABLE_NAME.cgemm_p * TABLE_NAME.cgemm_q * 8 + TABLE_NAME.offsetA
  974. + TABLE_NAME.align) & ~TABLE_NAME.align)
  975. ) / (TABLE_NAME.cgemm_q * 8) - 15) & ~15);
  976. TABLE_NAME.zgemm_r = (((BUFFER_SIZE -
  977. ((TABLE_NAME.zgemm_p * TABLE_NAME.zgemm_q * 16 + TABLE_NAME.offsetA
  978. + TABLE_NAME.align) & ~TABLE_NAME.align)
  979. ) / (TABLE_NAME.zgemm_q * 16) - 15) & ~15);
  980. TABLE_NAME.cgemm3m_r = (((BUFFER_SIZE -
  981. ((TABLE_NAME.cgemm3m_p * TABLE_NAME.cgemm3m_q * 8 + TABLE_NAME.offsetA
  982. + TABLE_NAME.align) & ~TABLE_NAME.align)
  983. ) / (TABLE_NAME.cgemm3m_q * 8) - 15) & ~15);
  984. TABLE_NAME.zgemm3m_r = (((BUFFER_SIZE -
  985. ((TABLE_NAME.zgemm3m_p * TABLE_NAME.zgemm3m_q * 16 + TABLE_NAME.offsetA
  986. + TABLE_NAME.align) & ~TABLE_NAME.align)
  987. ) / (TABLE_NAME.zgemm3m_q * 16) - 15) & ~15);
  988. #ifdef EXPRECISION
  989. TABLE_NAME.xgemm_r = (((BUFFER_SIZE -
  990. ((TABLE_NAME.xgemm_p * TABLE_NAME.xgemm_q * 32 + TABLE_NAME.offsetA
  991. + TABLE_NAME.align) & ~TABLE_NAME.align)
  992. ) / (TABLE_NAME.xgemm_q * 32) - 15) & ~15);
  993. TABLE_NAME.xgemm3m_r = (((BUFFER_SIZE -
  994. ((TABLE_NAME.xgemm3m_p * TABLE_NAME.xgemm3m_q * 32 + TABLE_NAME.offsetA
  995. + TABLE_NAME.align) & ~TABLE_NAME.align)
  996. ) / (TABLE_NAME.xgemm3m_q * 32) - 15) & ~15);
  997. #endif
  998. }