OSchip
/
OpenBLAS
Not watched
1
0
Fork 0
Code
Releases 66 Wiki evaluate Activity
Issues 0 Pull Requests 0 Datasets Model Cloudbrain HPC
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.
Branch: develop
Branches Tags
${ item.name }
Create branch ${ searchTerm }
from 'develop'
${ noResults }
OpenBLAS/utest/test_potrs.c

574 lines
28 kB
Raw Permalink Blame History 

  1. /*****************************************************************************

  2. Copyright (c) 2011-2016, The OpenBLAS Project

  3. All rights reserved.

  4. 

  5. Redistribution and use in source and binary forms, with or without

  6. modification, are permitted provided that the following conditions are

  7. met:

  8. 

  9.    1. Redistributions of source code must retain the above copyright

  10.       notice, this list of conditions and the following disclaimer.

  11. 

  12.    2. Redistributions in binary form must reproduce the above copyright

  13.       notice, this list of conditions and the following disclaimer in

  14.       the documentation and/or other materials provided with the

  15.       distribution.

  16.    3. Neither the name of the OpenBLAS project nor the names of 

  17.       its contributors may be used to endorse or promote products 

  18.       derived from this software without specific prior written 

  19.       permission.

  20. 

  21. THIS SOFTWARE IS PROVIDED BY THE COPYRIGHT HOLDERS AND CONTRIBUTORS "AS IS"

  22. AND ANY EXPRESS OR IMPLIED WARRANTIES, INCLUDING, BUT NOT LIMITED TO, THE

  23. IMPLIED WARRANTIES OF MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE

  24. ARE DISCLAIMED. IN NO EVENT SHALL THE COPYRIGHT OWNER OR CONTRIBUTORS BE

  25. LIABLE FOR ANY DIRECT, INDIRECT, INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL

  26. DAMAGES (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE GOODS OR

  27. SERVICES; LOSS OF USE, DATA, OR PROFITS; OR BUSINESS INTERRUPTION) HOWEVER

  28. CAUSED AND ON ANY THEORY OF LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY,

  29. OR TORT (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT OF THE

  30. USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE POSSIBILITY OF SUCH DAMAGE.

  31. 

  32. **********************************************************************************/

  33. 

  34. #include "openblas_utest.h"

  35. #if defined(ARCH_LOONGARCH64)

  36. #pragma GCC optimize("no-gcse")

  37. #endif

  38. /*

  39. void BLASFUNC(cpotrf)(char*, BLASINT*, complex float*, BLASINT*, BLASINT*);

  40. void BLASFUNC(zpotrs_(char*, BLASINT*, BLASINT*, complex double*,

  41. 	     BLASINT*, complex double*, BLASINT*, BLASINT*);

  42. */

  43. 

  44. //https://github.com/xianyi/OpenBLAS/issues/695

  45. CTEST(potrf, bug_695){

  46. 

  47.   openblas_complex_float A1[100] = 

  48.   {

  49.     openblas_make_complex_float(5.8525753, +0.0),

  50.     openblas_make_complex_float(-0.79540455, -0.7066077),

  51.     openblas_make_complex_float(0.98274714, -1.3824869),

  52.     openblas_make_complex_float(2.619998, -1.8532984),

  53.     openblas_make_complex_float(-1.8306153, +1.2336911),

  54.     openblas_make_complex_float(0.32275113, -0.015575029),

  55.     openblas_make_complex_float(2.1968813, -1.0640624),

  56.     openblas_make_complex_float(0.27894387, -0.97911835),

  57.     openblas_make_complex_float(3.0476584, -0.18548489),

  58.     openblas_make_complex_float(0.3842994, -0.7050991),

  59.     openblas_make_complex_float(-0.79540455, +0.7066077),

  60.     openblas_make_complex_float(8.313246, +0.0),

  61.     openblas_make_complex_float(-1.8076122, +0.8882447),

  62.     openblas_make_complex_float(0.47806996, -0.48494184),

  63.     openblas_make_complex_float(0.5096429, +0.5395974),

  64.     openblas_make_complex_float(-0.7285097, +0.10360408),

  65.     openblas_make_complex_float(-1.1760061, +2.7146957),

  66.     openblas_make_complex_float(-0.4271084, -0.042899966),

  67.     openblas_make_complex_float(-1.7228563, -2.8335886),

  68.     openblas_make_complex_float(1.8942566, -0.6389735),

  69.     openblas_make_complex_float(0.98274714, +1.3824869),

  70.     openblas_make_complex_float(-1.8076122, -0.8882447),

  71.     openblas_make_complex_float(9.367975, +0.0),

  72.     openblas_make_complex_float(-0.1838578, -0.6468568),

  73.     openblas_make_complex_float(-1.8338387, -0.7064959),

  74.     openblas_make_complex_float(0.041852742, +0.6556877),

  75.     openblas_make_complex_float(2.5673025, -1.9732997),

  76.     openblas_make_complex_float(-1.1148382, +0.15693812),

  77.     openblas_make_complex_float(2.4704504, +1.0389464),

  78.     openblas_make_complex_float(1.0858271, +1.298006),

  79.     openblas_make_complex_float(2.619998, +1.8532984),

  80.     openblas_make_complex_float(0.47806996, +0.48494184),

  81.     openblas_make_complex_float(-0.1838578, +0.6468568),

  82.     openblas_make_complex_float(3.1117508, +0.0),

  83.     openblas_make_complex_float(-1.956626, -0.22825956),

  84.     openblas_make_complex_float(0.07081801, +0.31801307),

  85.     openblas_make_complex_float(0.3698375, +0.5400855),

  86.     openblas_make_complex_float(0.80686307, -1.5315914),

  87.     openblas_make_complex_float(1.5649154, +1.6229297),

  88.     openblas_make_complex_float(-0.112077385, -1.2014246),

  89.     openblas_make_complex_float(-1.8306153, -1.2336911),

  90.     openblas_make_complex_float(0.5096429, -0.5395974),

  91.     openblas_make_complex_float(-1.8338387, +0.7064959),

  92.     openblas_make_complex_float(-1.956626, +0.22825956),

  93.     openblas_make_complex_float(3.6439795, +0.0),

  94.     openblas_make_complex_float(-0.2594722, -0.48786148),

  95.     openblas_make_complex_float(-0.47636223, +0.27821827),

  96.     openblas_make_complex_float(-0.61608654, +2.01858),

  97.     openblas_make_complex_float(-2.7767487, -1.7693765),

  98.     openblas_make_complex_float(0.048102796, +0.9741874),

  99.     openblas_make_complex_float(0.32275113, +0.015575029),

  100.     openblas_make_complex_float(-0.7285097, -0.10360408),

  101.     openblas_make_complex_float(0.041852742, -0.6556877),

  102.     openblas_make_complex_float(0.07081801, -0.31801307),

  103.     openblas_make_complex_float(-0.2594722, +0.48786148),

  104.     openblas_make_complex_float(3.624376, +0.0),

  105.     openblas_make_complex_float(-1.6697118, -0.4017511),

  106.     openblas_make_complex_float(-1.4397877, +0.7550918),

  107.     openblas_make_complex_float(-0.31456697, +1.0403451),

  108.     openblas_make_complex_float(-0.31978557, -0.13701046),

  109.     openblas_make_complex_float(2.1968813, +1.0640624),

  110.     openblas_make_complex_float(-1.1760061, -2.7146957),

  111.     openblas_make_complex_float(2.5673025, +1.9732997),

  112.     openblas_make_complex_float(0.3698375, -0.5400855),

  113.     openblas_make_complex_float(-0.47636223, -0.27821827),

  114.     openblas_make_complex_float(-1.6697118, +0.4017511),

  115.     openblas_make_complex_float(6.8273163, +0.0),

  116.     openblas_make_complex_float(-0.10051322, -0.24303961),

  117.     openblas_make_complex_float(1.4415971, -0.29750675),

  118.     openblas_make_complex_float(1.221786, +0.85654986),

  119.     openblas_make_complex_float(0.27894387, +0.97911835),

  120.     openblas_make_complex_float(-0.4271084, +0.042899966),

  121.     openblas_make_complex_float(-1.1148382, -0.15693812),

  122.     openblas_make_complex_float(0.80686307, +1.5315914),

  123.     openblas_make_complex_float(-0.61608654, -2.01858),

  124.     openblas_make_complex_float(-1.4397877, -0.7550918),

  125.     openblas_make_complex_float(-0.10051322, +0.24303961),

  126.     openblas_make_complex_float(3.4057708, +0.0),

  127.     openblas_make_complex_float(-0.5856801, +1.0203559),

  128.     openblas_make_complex_float(0.7103452, -0.8422135),

  129.     openblas_make_complex_float(3.0476584, +0.18548489),

  130.     openblas_make_complex_float(-1.7228563, +2.8335886),

  131.     openblas_make_complex_float(2.4704504, -1.0389464),

  132.     openblas_make_complex_float(1.5649154, -1.6229297),

  133.     openblas_make_complex_float(-2.7767487, +1.7693765),

  134.     openblas_make_complex_float(-0.31456697, -1.0403451),

  135.     openblas_make_complex_float(1.4415971, +0.29750675),

  136.     openblas_make_complex_float(-0.5856801, -1.0203559),

  137.     openblas_make_complex_float(7.005772, +0.0),

  138.     openblas_make_complex_float(-0.9617417, +1.2486815),

  139.     openblas_make_complex_float(0.3842994, +0.7050991),

  140.     openblas_make_complex_float(1.8942566, +0.6389735),

  141.     openblas_make_complex_float(1.0858271, -1.298006),

  142.     openblas_make_complex_float(-0.112077385, +1.2014246),

  143.     openblas_make_complex_float(0.048102796, -0.9741874),

  144.     openblas_make_complex_float(-0.31978557, +0.13701046),

  145.     openblas_make_complex_float(1.221786, -0.85654986),

  146.     openblas_make_complex_float(0.7103452, +0.8422135),

  147.     openblas_make_complex_float(-0.9617417, -1.2486815),

  148.     openblas_make_complex_float(3.4629636, +0.0)

  149.   };

  150. 

  151.   char up = 'U';

  152. 

  153.   blasint n=10;

  154.   blasint info[1];

  155. #ifdef BUILD_COMPLEX

  156.   BLASFUNC(cpotrf)(&up, &n, (float*)(A1), &n, info);

  157.   //printf("%g+%g*I\n", creal(A1[91]), cimag(A1[91]));

  158. #endif

  159. 

  160.   openblas_complex_double A2[100] = 

  161.   {

  162.     openblas_make_complex_double(3.0607147216796875, +0.0),

  163.     openblas_make_complex_double(-0.5905849933624268, -0.29020825028419495),

  164.     openblas_make_complex_double(0.321084201335907, +0.45168760418891907),

  165.     openblas_make_complex_double(0.8387917876243591, -0.644718587398529),

  166.     openblas_make_complex_double(-0.3642411530017853, +0.051274992525577545),

  167.     openblas_make_complex_double(0.8071482181549072, +0.33944568037986755),

  168.     openblas_make_complex_double(0.013674172572791576, +0.21422699093818665),

  169.     openblas_make_complex_double(0.35476258397102356, +0.42408594489097595),

  170.     openblas_make_complex_double(-0.5991537570953369, -0.23082709312438965),

  171.     openblas_make_complex_double(-0.0600702166557312, -0.2113417387008667),

  172.     openblas_make_complex_double(-0.7954045534133911, +0.7066076993942261),

  173.     openblas_make_complex_double(2.807175397872925, +0.0),

  174.     openblas_make_complex_double(-0.1691000759601593, +0.313548743724823),

  175.     openblas_make_complex_double(-0.30911174416542053, +0.7447023987770081),

  176.     openblas_make_complex_double(-0.22347848117351532, +0.03316075727343559),

  177.     openblas_make_complex_double(-0.4088296890258789, -1.0214389562606812),

  178.     openblas_make_complex_double(-0.2344931811094284, +0.08056317269802094),

  179.     openblas_make_complex_double(0.793269693851471, -0.17507623136043549),

  180.     openblas_make_complex_double(0.03163455054163933, +0.20559945702552795),

  181.     openblas_make_complex_double(0.13581633567810059, -0.2110036462545395),

  182.     openblas_make_complex_double(0.9827471375465393, +1.3824869394302368),

  183.     openblas_make_complex_double(-1.8076121807098389, -0.8882446885108948),

  184.     openblas_make_complex_double(2.3277781009674072, +0.0),

  185.     openblas_make_complex_double(0.830405056476593, -0.19296252727508545),

  186.     openblas_make_complex_double(0.1394239068031311, -0.5260677933692932),

  187.     openblas_make_complex_double(1.239942193031311, -0.09915469586849213),

  188.     openblas_make_complex_double(0.06731037050485611, -0.059320636093616486),

  189.     openblas_make_complex_double(0.11507681757211685, -0.1984301060438156),

  190.     openblas_make_complex_double(-0.6843825578689575, +0.4647614359855652),

  191.     openblas_make_complex_double(1.213119387626648, -0.7757048010826111),

  192.     openblas_make_complex_double(2.619997978210449, +1.8532984256744385),

  193.     openblas_make_complex_double(0.4780699610710144, +0.48494184017181396),

  194.     openblas_make_complex_double(-0.18385779857635498, +0.6468567848205566),

  195.     openblas_make_complex_double(2.0811400413513184, +0.0),

  196.     openblas_make_complex_double(-0.035075582563877106, +0.09732913225889206),

  197.     openblas_make_complex_double(0.27337002754211426, -0.9032229781150818),

  198.     openblas_make_complex_double(-0.8374675512313843, +0.0479498989880085),

  199.     openblas_make_complex_double(0.6916252374649048, +0.45711082220077515),

  200.     openblas_make_complex_double(0.1883818507194519, +0.06482727080583572),

  201.     openblas_make_complex_double(-0.32384994626045227, +0.05857187137007713),

  202.     openblas_make_complex_double(-1.8306152820587158, -1.2336910963058472),

  203.     openblas_make_complex_double(0.5096428990364075, -0.5395973920822144),

  204.     openblas_make_complex_double(-1.833838701248169, +0.7064958810806274),

  205.     openblas_make_complex_double(-1.956626057624817, +0.22825956344604492),

  206.     openblas_make_complex_double(1.706615924835205, +0.0),

  207.     openblas_make_complex_double(-0.2895336151123047, +0.17579378187656403),

  208.     openblas_make_complex_double(-0.923172116279602, -0.4530014097690582),

  209.     openblas_make_complex_double(0.5040621757507324, -0.37026339769363403),

  210.     openblas_make_complex_double(-0.2824432849884033, -1.0374568700790405),

  211.     openblas_make_complex_double(0.1399831622838974, +0.4977008104324341),

  212.     openblas_make_complex_double(0.32275113463401794, +0.015575028955936432),

  213.     openblas_make_complex_double(-0.7285097241401672, -0.10360407829284668),

  214.     openblas_make_complex_double(0.041852742433547974, -0.655687689781189),

  215.     openblas_make_complex_double(0.07081800699234009, -0.318013072013855),

  216.     openblas_make_complex_double(-0.25947219133377075, +0.4878614842891693),

  217.     openblas_make_complex_double(1.5735365152359009, +0.0),

  218.     openblas_make_complex_double(-0.2647853195667267, -0.26654252409935),

  219.     openblas_make_complex_double(-0.6190430521965027, -0.24699924886226654),

  220.     openblas_make_complex_double(-0.6288471221923828, +0.48154571652412415),

  221.     openblas_make_complex_double(0.02446540631353855, -0.2611822783946991),

  222.     openblas_make_complex_double(2.1968812942504883, +1.0640623569488525),

  223.     openblas_make_complex_double(-1.1760060787200928, -2.714695692062378),

  224.     openblas_make_complex_double(2.5673024654388428, +1.9732997417449951),

  225.     openblas_make_complex_double(0.3698374927043915, -0.54008549451828),

  226.     openblas_make_complex_double(-0.4763622283935547, -0.27821826934814453),

  227.     openblas_make_complex_double(-1.6697118282318115, +0.4017511010169983),

  228.     openblas_make_complex_double(1.2674795389175415, +0.0),

  229.     openblas_make_complex_double(0.3079095482826233, -0.07258892804384232),

  230.     openblas_make_complex_double(-0.5929520130157471, -0.038360968232154846),

  231.     openblas_make_complex_double(0.04388086497783661, -0.025549031794071198),

  232.     openblas_make_complex_double(0.27894386649131775, +0.9791183471679688),

  233.     openblas_make_complex_double(-0.42710840702056885, +0.0428999662399292),

  234.     openblas_make_complex_double(-1.1148382425308228, -0.1569381207227707),

  235.     openblas_make_complex_double(0.8068630695343018, +1.5315914154052734),

  236.     openblas_make_complex_double(-0.6160865426063538, -2.0185799598693848),

  237.     openblas_make_complex_double(-1.439787745475769, -0.7550917863845825),

  238.     openblas_make_complex_double(-0.10051321983337402, +0.24303960800170898),

  239.     openblas_make_complex_double(0.9066106081008911, +0.0),

  240.     openblas_make_complex_double(0.05315789580345154, -0.06136537343263626),

  241.     openblas_make_complex_double(-0.21304509043693542, +0.6494344472885132),

  242.     openblas_make_complex_double(3.0476584434509277, +0.1854848861694336),

  243.     openblas_make_complex_double(-1.7228562831878662, +2.8335886001586914),

  244.     openblas_make_complex_double(2.4704504013061523, -1.0389463901519775),

  245.     openblas_make_complex_double(1.564915418624878, -1.6229296922683716),

  246.     openblas_make_complex_double(-2.7767486572265625, +1.769376516342163),

  247.     openblas_make_complex_double(-0.314566969871521, -1.0403450727462769),

  248.     openblas_make_complex_double(1.4415971040725708, +0.29750674962997437),

  249.     openblas_make_complex_double(-0.5856801271438599, -1.0203559398651123),

  250.     openblas_make_complex_double(0.5668219923973083, +0.0),

  251.     openblas_make_complex_double(0.033351436257362366, -0.07832501083612442),

  252.     openblas_make_complex_double(0.3842993974685669, +0.7050991058349609),

  253.     openblas_make_complex_double(1.894256591796875, +0.6389734745025635),

  254.     openblas_make_complex_double(1.085827112197876, -1.2980060577392578),

  255.     openblas_make_complex_double(-0.11207738518714905, +1.2014245986938477),

  256.     openblas_make_complex_double(0.04810279607772827, -0.9741873741149902),

  257.     openblas_make_complex_double(-0.31978556513786316, +0.13701045513153076),

  258.     openblas_make_complex_double(1.2217860221862793, -0.856549859046936),

  259.     openblas_make_complex_double(0.7103452086448669, +0.84221351146698),

  260.     openblas_make_complex_double(-0.9617416858673096, -1.2486815452575684),

  261.     openblas_make_complex_double(0.0756804421544075, +0.0)

  262.   };

  263.   openblas_complex_double B[20] = 

  264.   {

  265.     openblas_make_complex_double(-0.21782716937787788, -0.9222220085490986),

  266.     openblas_make_complex_double(-0.7620356655676837, +0.15533508334193666),

  267.     openblas_make_complex_double(-0.905011814118756, +0.2847570854574069),

  268.     openblas_make_complex_double(-0.3451346708401685, +1.076948486041297),

  269.     openblas_make_complex_double(0.25336108035924787, +0.975317836492159),

  270.     openblas_make_complex_double(0.11192755545114, -0.1603741874112385),

  271.     openblas_make_complex_double(-0.20604111555491242, +0.10570814584017311),

  272.     openblas_make_complex_double(-1.0568488936791578, -0.06025820467086475),

  273.     openblas_make_complex_double(-0.6650468984506477, -0.5000967284800251),

  274.     openblas_make_complex_double(-1.0509472322215125, +0.5022165705328413),

  275.     openblas_make_complex_double(-0.727775859267237, +0.50638268521728),

  276.     openblas_make_complex_double(0.39947219167701153, -0.4576746001199889),

  277.     openblas_make_complex_double(-0.7122162951294634, -0.630289556702497),

  278.     openblas_make_complex_double(0.9870834574024372, -0.2825689605519449),

  279.     openblas_make_complex_double(0.0628393808469436, -0.1253397353973715),

  280.     openblas_make_complex_double(0.8439562576196216, +1.0850814110398734),

  281.     openblas_make_complex_double(0.562377322638969, -0.2578030745663871),

  282.     openblas_make_complex_double(0.12696236014017806, -0.09853584666755086),

  283.     openblas_make_complex_double(-0.023682508769195098, +0.18093440285319276),

  284.     openblas_make_complex_double(-0.7264975746431271, +0.31670415674097235)

  285.   };

  286.   char lo = 'L';

  287.   blasint nrhs = 2;

  288. #ifdef BUILD_COMPLEX16

  289.   BLASFUNC(zpotrs)(&lo, &n, &nrhs, (double*)(A2), &n, (double*)(B), &n, info);

  290. #endif

  291.   // note that this is exactly equal to A1

  292.   openblas_complex_float A3[100] = 

  293.   {

  294.     openblas_make_complex_float(5.8525753, +0.0),

  295.     openblas_make_complex_float(-0.79540455, -0.7066077),

  296.     openblas_make_complex_float(0.98274714, -1.3824869),

  297.     openblas_make_complex_float(2.619998, -1.8532984),

  298.     openblas_make_complex_float(-1.8306153, +1.2336911),

  299.     openblas_make_complex_float(0.32275113, -0.015575029),

  300.     openblas_make_complex_float(2.1968813, -1.0640624),

  301.     openblas_make_complex_float(0.27894387, -0.97911835),

  302.     openblas_make_complex_float(3.0476584, -0.18548489),

  303.     openblas_make_complex_float(0.3842994, -0.7050991),

  304.     openblas_make_complex_float(-0.79540455, +0.7066077),

  305.     openblas_make_complex_float(8.313246, +0.0),

  306.     openblas_make_complex_float(-1.8076122, +0.8882447),

  307.     openblas_make_complex_float(0.47806996, -0.48494184),

  308.     openblas_make_complex_float(0.5096429, +0.5395974),

  309.     openblas_make_complex_float(-0.7285097, +0.10360408),

  310.     openblas_make_complex_float(-1.1760061, +2.7146957),

  311.     openblas_make_complex_float(-0.4271084, -0.042899966),

  312.     openblas_make_complex_float(-1.7228563, -2.8335886),

  313.     openblas_make_complex_float(1.8942566, -0.6389735),

  314.     openblas_make_complex_float(0.98274714, +1.3824869),

  315.     openblas_make_complex_float(-1.8076122, -0.8882447),

  316.     openblas_make_complex_float(9.367975, +0.0),

  317.     openblas_make_complex_float(-0.1838578, -0.6468568),

  318.     openblas_make_complex_float(-1.8338387, -0.7064959),

  319.     openblas_make_complex_float(0.041852742, +0.6556877),

  320.     openblas_make_complex_float(2.5673025, -1.9732997),

  321.     openblas_make_complex_float(-1.1148382, +0.15693812),

  322.     openblas_make_complex_float(2.4704504, +1.0389464),

  323.     openblas_make_complex_float(1.0858271, +1.298006),

  324.     openblas_make_complex_float(2.619998, +1.8532984),

  325.     openblas_make_complex_float(0.47806996, +0.48494184),

  326.     openblas_make_complex_float(-0.1838578, +0.6468568),

  327.     openblas_make_complex_float(3.1117508, +0.0),

  328.     openblas_make_complex_float(-1.956626, -0.22825956),

  329.     openblas_make_complex_float(0.07081801, +0.31801307),

  330.     openblas_make_complex_float(0.3698375, +0.5400855),

  331.     openblas_make_complex_float(0.80686307, -1.5315914),

  332.     openblas_make_complex_float(1.5649154, +1.6229297),

  333.     openblas_make_complex_float(-0.112077385, -1.2014246),

  334.     openblas_make_complex_float(-1.8306153, -1.2336911),

  335.     openblas_make_complex_float(0.5096429, -0.5395974),

  336.     openblas_make_complex_float(-1.8338387, +0.7064959),

  337.     openblas_make_complex_float(-1.956626, +0.22825956),

  338.     openblas_make_complex_float(3.6439795, +0.0),

  339.     openblas_make_complex_float(-0.2594722, -0.48786148),

  340.     openblas_make_complex_float(-0.47636223, +0.27821827),

  341.     openblas_make_complex_float(-0.61608654, +2.01858),

  342.     openblas_make_complex_float(-2.7767487, -1.7693765),

  343.     openblas_make_complex_float(0.048102796, +0.9741874),

  344.     openblas_make_complex_float(0.32275113, +0.015575029),

  345.     openblas_make_complex_float(-0.7285097, -0.10360408),

  346.     openblas_make_complex_float(0.041852742, -0.6556877),

  347.     openblas_make_complex_float(0.07081801, -0.31801307),

  348.     openblas_make_complex_float(-0.2594722, +0.48786148),

  349.     openblas_make_complex_float(3.624376, +0.0),

  350.     openblas_make_complex_float(-1.6697118, -0.4017511),

  351.     openblas_make_complex_float(-1.4397877, +0.7550918),

  352.     openblas_make_complex_float(-0.31456697, +1.0403451),

  353.     openblas_make_complex_float(-0.31978557, -0.13701046),

  354.     openblas_make_complex_float(2.1968813, +1.0640624),

  355.     openblas_make_complex_float(-1.1760061, -2.7146957),

  356.     openblas_make_complex_float(2.5673025, +1.9732997),

  357.     openblas_make_complex_float(0.3698375, -0.5400855),

  358.     openblas_make_complex_float(-0.47636223, -0.27821827),

  359.     openblas_make_complex_float(-1.6697118, +0.4017511),

  360.     openblas_make_complex_float(6.8273163, +0.0),

  361.     openblas_make_complex_float(-0.10051322, -0.24303961),

  362.     openblas_make_complex_float(1.4415971, -0.29750675),

  363.     openblas_make_complex_float(1.221786, +0.85654986),

  364.     openblas_make_complex_float(0.27894387, +0.97911835),

  365.     openblas_make_complex_float(-0.4271084, +0.042899966),

  366.     openblas_make_complex_float(-1.1148382, -0.15693812),

  367.     openblas_make_complex_float(0.80686307, +1.5315914),

  368.     openblas_make_complex_float(-0.61608654, -2.01858),

  369.     openblas_make_complex_float(-1.4397877, -0.7550918),

  370.     openblas_make_complex_float(-0.10051322, +0.24303961),

  371.     openblas_make_complex_float(3.4057708, +0.0),

  372.     openblas_make_complex_float(-0.5856801, +1.0203559),

  373.     openblas_make_complex_float(0.7103452, -0.8422135),

  374.     openblas_make_complex_float(3.0476584, +0.18548489),

  375.     openblas_make_complex_float(-1.7228563, +2.8335886),

  376.     openblas_make_complex_float(2.4704504, -1.0389464),

  377.     openblas_make_complex_float(1.5649154, -1.6229297),

  378.     openblas_make_complex_float(-2.7767487, +1.7693765),

  379.     openblas_make_complex_float(-0.31456697, -1.0403451),

  380.     openblas_make_complex_float(1.4415971, +0.29750675),

  381.     openblas_make_complex_float(-0.5856801, -1.0203559),

  382.     openblas_make_complex_float(7.005772, +0.0),

  383.     openblas_make_complex_float(-0.9617417, +1.2486815),

  384.     openblas_make_complex_float(0.3842994, +0.7050991),

  385.     openblas_make_complex_float(1.8942566, +0.6389735),

  386.     openblas_make_complex_float(1.0858271, -1.298006),

  387.     openblas_make_complex_float(-0.112077385, +1.2014246),

  388.     openblas_make_complex_float(0.048102796, -0.9741874),

  389.     openblas_make_complex_float(-0.31978557, +0.13701046),

  390.     openblas_make_complex_float(1.221786, -0.85654986),

  391.     openblas_make_complex_float(0.7103452, +0.8422135),

  392.     openblas_make_complex_float(-0.9617417, -1.2486815),

  393.     openblas_make_complex_float(3.4629636, +0.0)

  394.   };

  395. #ifdef BUILD_COMPLEX

  396.   BLASFUNC(cpotrf)(&up, &n, (float*)(A3), &n, info);

  397.   //  printf("%g+%g*I\n", creal(A3[91]), cimag(A3[91]));

  398.   if(isnan(CREAL(A3[91])) || isnan(CIMAG(A3[91]))) {

  399.     CTEST_ERR("%s:%d  got NaN", __FILE__, __LINE__);

  400.   }

  401. #endif

  402. }

  403. 

  404. // Check potrf factorizes a small problem correctly

  405. CTEST(potrf, smoketest_trivial){

  406.   float A1s[4] = {2, 0.3, 0.3, 3};

  407.   double A1d[4] = {2, 0.3, 0.3, 3};

  408.   openblas_complex_float A1c[4] = {

  409.     openblas_make_complex_float(2,0),

  410.     openblas_make_complex_float(0.3,0.1),

  411.     openblas_make_complex_float(0.3,-0.1),

  412.     openblas_make_complex_float(3,0)

  413.   };

  414.   openblas_complex_double A1z[4] = {

  415.     openblas_make_complex_double(2,0),

  416.     openblas_make_complex_double(0.3,0.1),

  417.     openblas_make_complex_double(0.3,-0.1),

  418.     openblas_make_complex_double(3,0)

  419.   };

  420.   float zeros = 0, ones = 1;

  421.   double zerod = 0, oned = 1;

  422.   openblas_complex_float zeroc = openblas_make_complex_float(0, 0),

  423.                          onec = openblas_make_complex_float(1, 0);

  424.   openblas_complex_double zeroz = openblas_make_complex_double(0, 0),

  425.                           onez = openblas_make_complex_float(1, 0);

  426. 

  427.   char uplo, trans1, trans2;

  428.   blasint nv = 4;

  429.   blasint n = 2;

  430.   blasint inc = 1;

  431.   blasint info = 0;

  432.   int i, j, cycle;

  433. 

  434.   float As[4], Bs[4];

  435.   double Ad[4], Bd[4];

  436.   openblas_complex_float Ac[4], Bc[4];

  437.   openblas_complex_double Az[4], Bz[4];

  438. 

  439.   for (cycle = 0; cycle < 2; ++cycle) {

  440.     if (cycle == 0) {

  441.       uplo = 'L';

  442.     }

  443.     else {

  444.       uplo = 'U';

  445.     }

  446. 

  447. #ifdef BUILD_SINGLE

  448.     BLASFUNC(scopy)(&nv, A1s, &inc, As, &inc);

  449. #endif

  450. #ifdef BUILD_DOUBLE

  451.     BLASFUNC(dcopy)(&nv, A1d, &inc, Ad, &inc);

  452. #endif

  453. #ifdef BUILD_COMPLEX

  454.     BLASFUNC(ccopy)(&nv, (float *)A1c, &inc, (float *)Ac, &inc);

  455. #endif

  456. #ifdef BUILD_COMPLEX16

  457.     BLASFUNC(zcopy)(&nv, (double *)A1z, &inc, (double *)Az, &inc);

  458. #endif

  459. 

  460. #ifdef BUILD_SINGLE

  461.     BLASFUNC(spotrf)(&uplo, &n, As, &n, &info);

  462.     if (info != 0) {

  463.       CTEST_ERR("%s:%d  info != 0", __FILE__, __LINE__);

  464.     }

  465. #endif

  466. #ifdef BUILD_DOUBLE

  467.     BLASFUNC(dpotrf)(&uplo, &n, Ad, &n, &info);

  468.     if (info != 0) {

  469.       CTEST_ERR("%s:%d  info != 0", __FILE__, __LINE__);

  470.     }

  471. #endif

  472. #ifdef BUILD_COMPLEX

  473.     BLASFUNC(cpotrf)(&uplo, &n, (float *)Ac, &n, &info);

  474.     if (info != 0) {

  475.       CTEST_ERR("%s:%d  info != 0", __FILE__, __LINE__);

  476.     }

  477. #endif

  478. #ifdef BUILD_COMPLEX16

  479.     BLASFUNC(zpotrf)(&uplo, &n, (double *)Az, &n, &info);

  480.     if (info != 0) {

  481.       CTEST_ERR("%s:%d  info != 0", __FILE__, __LINE__);

  482.     }

  483. #endif

  484.     /* Fill the other triangle */

  485.     if (uplo == 'L') {

  486.       for (i = 0; i < n; ++i) {

  487.         for (j = i+1; j < n; ++j) {

  488.           As[i+n*j] = 0;

  489.           Ad[i+n*j] = 0;

  490.           Ac[i+n*j] = zeroc;

  491.           Az[i+n*j] = zeroz;

  492.         }

  493.       }

  494.     }

  495.     else {

  496.       for (i = 0; i < n; ++i) {

  497.         for (j = 0; j < i; ++j) {

  498.           As[i+n*j] = 0;

  499.           Ad[i+n*j] = 0;

  500.           Ac[i+n*j] = zeroc;

  501.           Az[i+n*j] = zeroz;

  502.         }

  503.       }

  504.     }

  505. 

  506.     /* B = A A^H or A^H A */

  507.     if (uplo == 'L') {

  508.       trans1 = 'N';

  509.       trans2 = 'C';

  510.     }

  511.     else {

  512.       trans1 = 'C';

  513.       trans2 = 'N';

  514.     }

  515. #ifdef BUILD_SINGLE

  516.     BLASFUNC(sgemm)(&trans1, &trans2, &n, &n, &n, &ones, As, &n, As, &n, &zeros, Bs, &n);

  517. #endif

  518. #ifdef BUILD_DOUBLE

  519.     BLASFUNC(dgemm)(&trans1, &trans2, &n, &n, &n, &oned, Ad, &n, Ad, &n, &zerod, Bd, &n);

  520. #endif

  521. #ifdef BUILD_COMPLEX

  522.     BLASFUNC(cgemm)(&trans1, &trans2, &n, &n, &n, (float *)&onec,

  523.                     (float *)Ac, &n, (float *)Ac, &n, (float *)&zeroc, (float *)Bc, &n);

  524. #endif

  525. #ifdef BUILD_COMPLEX16

  526.     BLASFUNC(zgemm)(&trans1, &trans2, &n, &n, &n, (double *)&onez,

  527.                     (double *)Az, &n, (double *)Az, &n, (double *)&zeroz, (double *)Bz, &n);

  528. #endif

  529.     /* Check result is close to original */

  530.     for (i = 0; i < n; ++i) {

  531.       for (j = 0; j < n; ++j) {

  532.         double err;

  533. 

  534. #ifdef BUILD_SINGLE

  535.         err = fabs(A1s[i+n*j] - Bs[i+n*j]);

  536.         if (err > 1e-5) {

  537.           CTEST_ERR("%s:%d  %c s(%d,%d) difference: %g", __FILE__, __LINE__, uplo, i, j, err);

  538.         }

  539. #endif

  540. #ifdef BUILD_DOUBLE

  541.         err = fabs(A1d[i+n*j] - Bd[i+n*j]);

  542.         if (err > 1e-12) {

  543.           CTEST_ERR("%s:%d  %c d(%d,%d) difference: %g", __FILE__, __LINE__, uplo, i, j, err);

  544.         }

  545. #endif

  546. 

  547. #ifdef BUILD_COMPLEX

  548. #ifdef OPENBLAS_COMPLEX_C99

  549.         err = cabsf(A1c[i+n*j] - Bc[i+n*j]);

  550. #else

  551.         err = hypot(A1c[i+n*j].real - Bc[i+n*j].real,

  552.                     A1c[i+n*j].imag - Bc[i+n*j].imag);

  553. #endif

  554.         if (err > 1e-5) {

  555.           CTEST_ERR("%s:%d  %c c(%d,%d) difference: %g", __FILE__, __LINE__, uplo, i, j, err);

  556.         }

  557. #endif

  558. 

  559. #ifdef BUILD_COMPLEX16

  560. #ifdef OPENBLAS_COMPLEX_C99

  561.         err = cabs(A1z[i+n*j] - Bz[i+n*j]);

  562. #else

  563.         err = hypot(A1z[i+n*j].real - Bz[i+n*j].real,

  564.                     A1z[i+n*j].imag - Bz[i+n*j].imag);

  565. #endif

  566.         if (err > 1e-12) {

  567.           CTEST_ERR("%s:%d  %c z(%d,%d) difference: %g", __FILE__, __LINE__, uplo, i, j, err);

  568.         }

  569. #endif

  570.       }

  571.     }

  572.   }

  573. }