1 Star 0 Fork 1

flice/OpenBLAS

加入 Gitee
与超过 1200万 开发者一起发现、参与优秀开源项目,私有仓库也完全免费 :)
免费加入
文件
克隆/下载
symcopy.h 29.46 KB
一键复制 编辑 原始数据 按行查看 历史
1234567891011121314151617181920212223242526272829303132333435363738394041424344454647484950515253545556575859606162636465666768697071727374757677787980818283848586878889909192939495969798991001011021031041051061071081091101111121131141151161171181191201211221231241251261271281291301311321331341351361371381391401411421431441451461471481491501511521531541551561571581591601611621631641651661671681691701711721731741751761771781791801811821831841851861871881891901911921931941951961971981992002012022032042052062072082092102112122132142152162172182192202212222232242252262272282292302312322332342352362372382392402412422432442452462472482492502512522532542552562572582592602612622632642652662672682692702712722732742752762772782792802812822832842852862872882892902912922932942952962972982993003013023033043053063073083093103113123133143153163173183193203213223233243253263273283293303313323333343353363373383393403413423433443453463473483493503513523533543553563573583593603613623633643653663673683693703713723733743753763773783793803813823833843853863873883893903913923933943953963973983994004014024034044054064074084094104114124134144154164174184194204214224234244254264274284294304314324334344354364374384394404414424434444454464474484494504514524534544554564574584594604614624634644654664674684694704714724734744754764774784794804814824834844854864874884894904914924934944954964974984995005015025035045055065075085095105115125135145155165175185195205215225235245255265275285295305315325335345355365375385395405415425435445455465475485495505515525535545555565575585595605615625635645655665675685695705715725735745755765775785795805815825835845855865875885895905915925935945955965975985996006016026036046056066076086096106116126136146156166176186196206216226236246256266276286296306316326336346356366376386396406416426436446456466476486496506516526536546556566576586596606616626636646656666676686696706716726736746756766776786796806816826836846856866876886896906916926936946956966976986997007017027037047057067077087097107117127137147157167177187197207217227237247257267277287297307317327337347357367377387397407417427437447457467477487497507517527537547557567577587597607617627637647657667677687697707717727737747757767777787797807817827837847857867877887897907917927937947957967977987998008018028038048058068078088098108118128138148158168178188198208218228238248258268278288298308318328338348358368378388398408418428438448458468478488498508518528538548558568578588598608618628638648658668678688698708718728738748758768778788798808818828838848858868878888898908918928938948958968978988999009019029039049059069079089099109119129139149159169179189199209219229239249259269279289299309319329339349359369379389399409419429439449459469479489499509519529539549559569579589599609619629639649659669679689699709719729739749759769779789799809819829839849859869879889899909919929939949959969979989991000100110021003100410051006100710081009101010111012101310141015101610171018101910201021102210231024102510261027102810291030103110321033103410351036103710381039104010411042104310441045104610471048104910501051105210531054105510561057105810591060106110621063106410651066106710681069107010711072107310741075107610771078107910801081108210831084108510861087108810891090109110921093109410951096109710981099110011011102110311041105110611071108110911101111111211131114111511161117111811191120112111221123112411251126112711281129113011311132113311341135113611371138113911401141114211431144114511461147114811491150115111521153115411551156115711581159116011611162116311641165116611671168116911701171117211731174117511761177117811791180118111821183118411851186118711881189119011911192119311941195119611971198119912001201120212031204120512061207120812091210121112121213121412151216121712181219122012211222122312241225122612271228122912301231123212331234123512361237123812391240124112421243124412451246124712481249125012511252125312541255125612571258125912601261126212631264126512661267126812691270127112721273127412751276127712781279128012811282128312841285128612871288128912901291129212931294129512961297129812991300130113021303130413051306130713081309131013111312131313141315131613171318131913201321132213231324132513261327132813291330133113321333133413351336133713381339134013411342134313441345134613471348134913501351135213531354135513561357135813591360136113621363136413651366136713681369137013711372137313741375137613771378137913801381138213831384138513861387138813891390139113921393139413951396139713981399140014011402140314041405140614071408140914101411141214131414141514161417141814191420142114221423142414251426142714281429143014311432143314341435143614371438143914401441144214431444144514461447144814491450145114521453145414551456145714581459146014611462146314641465146614671468146914701471147214731474147514761477147814791480148114821483148414851486148714881489149014911492149314941495149614971498149915001501150215031504150515061507150815091510151115121513151415151516151715181519152015211522152315241525152615271528152915301531153215331534153515361537153815391540154115421543154415451546154715481549155015511552155315541555155615571558155915601561156215631564156515661567156815691570157115721573157415751576157715781579158015811582158315841585158615871588158915901591159215931594159515961597159815991600160116021603160416051606160716081609161016111612161316141615161616171618161916201621162216231624162516261627162816291630163116321633163416351636163716381639164016411642164316441645164616471648164916501651165216531654165516561657165816591660166116621663166416651666166716681669167016711672167316741675167616771678167916801681168216831684168516861687168816891690169116921693169416951696169716981699170017011702170317041705170617071708170917101711171217131714171517161717171817191720172117221723172417251726172717281729173017311732173317341735173617371738173917401741174217431744174517461747174817491750175117521753175417551756175717581759176017611762176317641765176617671768176917701771177217731774177517761777177817791780178117821783178417851786178717881789179017911792179317941795179617971798179918001801180218031804180518061807180818091810181118121813181418151816181718181819182018211822182318241825182618271828182918301831183218331834183518361837183818391840184118421843184418451846184718481849185018511852185318541855185618571858185918601861186218631864186518661867186818691870187118721873
/*********************************************************************/
/* Copyright 2009, 2010 The University of Texas at Austin. */
/* All rights reserved. */
/* */
/* Redistribution and use in source and binary forms, with or */
/* without modification, are permitted provided that the following */
/* conditions are met: */
/* */
/* 1. Redistributions of source code must retain the above */
/* copyright notice, this list of conditions and the following */
/* disclaimer. */
/* */
/* 2. Redistributions in binary form must reproduce the above */
/* copyright notice, this list of conditions and the following */
/* disclaimer in the documentation and/or other materials */
/* provided with the distribution. */
/* */
/* THIS SOFTWARE IS PROVIDED BY THE UNIVERSITY OF TEXAS AT */
/* AUSTIN ``AS IS'' AND ANY EXPRESS OR IMPLIED WARRANTIES, */
/* INCLUDING, BUT NOT LIMITED TO, THE IMPLIED WARRANTIES OF */
/* MERCHANTABILITY AND FITNESS FOR A PARTICULAR PURPOSE ARE */
/* DISCLAIMED. IN NO EVENT SHALL THE UNIVERSITY OF TEXAS AT */
/* AUSTIN OR CONTRIBUTORS BE LIABLE FOR ANY DIRECT, INDIRECT, */
/* INCIDENTAL, SPECIAL, EXEMPLARY, OR CONSEQUENTIAL DAMAGES */
/* (INCLUDING, BUT NOT LIMITED TO, PROCUREMENT OF SUBSTITUTE */
/* GOODS OR SERVICES; LOSS OF USE, DATA, OR PROFITS; OR */
/* BUSINESS INTERRUPTION) HOWEVER CAUSED AND ON ANY THEORY OF */
/* LIABILITY, WHETHER IN CONTRACT, STRICT LIABILITY, OR TORT */
/* (INCLUDING NEGLIGENCE OR OTHERWISE) ARISING IN ANY WAY OUT */
/* OF THE USE OF THIS SOFTWARE, EVEN IF ADVISED OF THE */
/* POSSIBILITY OF SUCH DAMAGE. */
/* */
/* The views and conclusions contained in the software and */
/* documentation are those of the authors and should not be */
/* interpreted as representing official policies, either expressed */
/* or implied, of The University of Texas at Austin. */
/*********************************************************************/
/* This implementation is completely wrong. I'll rewrite this */
#ifndef SYMCOPY_H
#define SYMCOPY_H
#if !defined(XDOUBLE) || !defined(QUAD_PRECISION)
static __inline void SYMCOPY_L(BLASLONG m, FLOAT *a, BLASLONG lda, FLOAT *b){
BLASLONG is, js;
FLOAT *aa1, *aa2;
FLOAT *b1, *b2;
FLOAT *bb1, *bb2;
FLOAT *cc1, *cc2;
FLOAT a11, a12;
FLOAT a21, a22;
b1 = b;
b2 = b;
for (js = 0; js < m; js += 2){
aa1 = a + 0 * lda;
aa2 = a + 1 * lda;
a += 2 * lda + 2;
bb1 = b1 + 0 * m;
bb2 = b1 + 1 * m;
b1 += 2 * m + 2;
cc1 = b2 + 0 * m;
cc2 = b2 + 1 * m;
b2 += 2 * m + 2;
if (m - js >= 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a22 = *(aa2 + 1);
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb2 + 0) = a21;
*(bb2 + 1) = a22;
aa1 += 2;
aa2 += 2;
bb1 += 2;
bb2 += 2;
cc1 += 2 * m;
cc2 += 2 * m;
is = ((m - js - 2) >> 1);
while (is > 0){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
aa1 += 2;
aa2 += 2;
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(cc1 + 0) = a11;
*(cc1 + 1) = a12;
*(cc2 + 0) = a21;
*(cc2 + 1) = a22;
bb1 += 2;
bb2 += 2;
cc1 += 2 * m;
cc2 += 2 * m;
is --;
}
is = ((m - js - 2) & 1);
if (is == 1){
a11 = *(aa1 + 0);
a12 = *(aa2 + 0);
*(bb1 + 0) = a11;
*(bb2 + 0) = a12;
*(cc1 + 0) = a11;
*(cc1 + 1) = a12;
}
}
if (m - js == 1){
a11 = *(aa1 + 0);
*(bb1 + 0) = a11;
}
}
}
static __inline void SYMCOPY_U(BLASLONG m, FLOAT *a, BLASLONG lda, FLOAT *b){
BLASLONG is, js;
FLOAT *aa1, *aa2;
FLOAT *b1, *b2;
FLOAT *bb1, *bb2;
FLOAT *cc1, *cc2;
FLOAT a11, a12;
FLOAT a21, a22;
b1 = b;
b2 = b;
for (js = 0; js < m; js += 2){
aa1 = a + 0 * lda;
aa2 = a + 1 * lda;
a += 2 * lda;
bb1 = b1 + 0 * m;
bb2 = b1 + 1 * m;
b1 += 2 * m;
cc1 = b2 + 0 * m;
cc2 = b2 + 1 * m;
b2 += 2;
if (m - js >= 2){
for (is = 0; is < js; is += 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
aa1 += 2;
aa2 += 2;
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(cc1 + 0) = a11;
*(cc1 + 1) = a12;
*(cc2 + 0) = a21;
*(cc2 + 1) = a22;
bb1 += 2;
bb2 += 2;
cc1 += 2 * m;
cc2 += 2 * m;
}
a11 = *(aa1 + 0);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
*(bb1 + 0) = a11;
*(bb1 + 1) = a12;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
}
if (m - js == 1){
for (is = 0; is < js; is += 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
aa1 += 2;
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(cc1 + 0) = a11;
*(cc2 + 0) = a21;
bb1 += 2;
cc1 += 2 * m;
cc2 += 2 * m;
}
a11 = *(aa1 + 0);
*(bb1 + 0) = a11;
}
}
}
static __inline void ZSYMCOPY_L(BLASLONG m, FLOAT *a, BLASLONG lda, FLOAT *b){
BLASLONG is, js;
FLOAT *aa1, *aa2;
FLOAT *b1, *b2;
FLOAT *bb1, *bb2;
FLOAT *cc1, *cc2;
FLOAT a11, a21, a31, a41;
FLOAT a12, a22, a32, a42;
b1 = b;
b2 = b;
lda *= 2;
for (js = 0; js < m; js += 2){
aa1 = a + 0 * lda;
aa2 = a + 1 * lda;
a += 2 * lda + 4;
bb1 = b1 + 0 * m;
bb2 = b1 + 2 * m;
b1 += 4 * m + 4;
cc1 = b2 + 0 * m;
cc2 = b2 + 2 * m;
b2 += 4 * m + 4;
if (m - js >= 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a31 = *(aa1 + 2);
a41 = *(aa1 + 3);
a12 = *(aa2 + 2);
a22 = *(aa2 + 3);
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb1 + 2) = a31;
*(bb1 + 3) = a41;
*(bb2 + 0) = a31;
*(bb2 + 1) = a41;
*(bb2 + 2) = a12;
*(bb2 + 3) = a22;
aa1 += 4;
aa2 += 4;
bb1 += 4;
bb2 += 4;
cc1 += 4 * m;
cc2 += 4 * m;
is = ((m - js - 2) >> 1);
while (is > 0){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a31 = *(aa1 + 2);
a41 = *(aa1 + 3);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
a32 = *(aa2 + 2);
a42 = *(aa2 + 3);
aa1 += 4;
aa2 += 4;
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb1 + 2) = a31;
*(bb1 + 3) = a41;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(bb2 + 2) = a32;
*(bb2 + 3) = a42;
*(cc1 + 0) = a11;
*(cc1 + 1) = a21;
*(cc1 + 2) = a12;
*(cc1 + 3) = a22;
*(cc2 + 0) = a31;
*(cc2 + 1) = a41;
*(cc2 + 2) = a32;
*(cc2 + 3) = a42;
bb1 += 4;
bb2 += 4;
cc1 += 4 * m;
cc2 += 4 * m;
is --;
}
if (m & 1){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(cc1 + 0) = a11;
*(cc1 + 1) = a21;
*(cc1 + 2) = a12;
*(cc1 + 3) = a22;
}
}
if (m - js == 1){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
}
}
}
static __inline void ZSYMCOPY_U(BLASLONG m, FLOAT *a, BLASLONG lda, FLOAT *b){
BLASLONG is, js;
FLOAT *aa1, *aa2;
FLOAT *b1, *b2;
FLOAT *bb1, *bb2;
FLOAT *cc1, *cc2;
FLOAT a11, a21, a31, a41;
FLOAT a12, a22, a32, a42;
b1 = b;
b2 = b;
lda *= 2;
for (js = 0; js < m; js += 2){
aa1 = a + 0 * lda;
aa2 = a + 1 * lda;
a += 2 * lda;
bb1 = b1 + 0 * m;
bb2 = b1 + 2 * m;
b1 += 4 * m;
cc1 = b2 + 0 * m;
cc2 = b2 + 2 * m;
b2 += 4;
if (m - js >= 2){
for (is = 0; is < js; is += 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a31 = *(aa1 + 2);
a41 = *(aa1 + 3);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
a32 = *(aa2 + 2);
a42 = *(aa2 + 3);
aa1 += 4;
aa2 += 4;
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb1 + 2) = a31;
*(bb1 + 3) = a41;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(bb2 + 2) = a32;
*(bb2 + 3) = a42;
*(cc1 + 0) = a11;
*(cc1 + 1) = a21;
*(cc1 + 2) = a12;
*(cc1 + 3) = a22;
*(cc2 + 0) = a31;
*(cc2 + 1) = a41;
*(cc2 + 2) = a32;
*(cc2 + 3) = a42;
bb1 += 4;
bb2 += 4;
cc1 += 4 * m;
cc2 += 4 * m;
}
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
a32 = *(aa2 + 2);
a42 = *(aa2 + 3);
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb1 + 2) = a12;
*(bb1 + 3) = a22;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(bb2 + 2) = a32;
*(bb2 + 3) = a42;
}
if (m - js == 1){
for (is = 0; is < js; is += 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a31 = *(aa1 + 2);
a41 = *(aa1 + 3);
aa1 += 4;
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb1 + 2) = a31;
*(bb1 + 3) = a41;
*(cc1 + 0) = a11;
*(cc1 + 1) = a21;
*(cc2 + 0) = a31;
*(cc2 + 1) = a41;
bb1 += 4;
cc1 += 4 * m;
cc2 += 4 * m;
}
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
}
}
}
static __inline void ZHEMCOPY_L(BLASLONG m, FLOAT *a, BLASLONG lda, FLOAT *b){
BLASLONG is, js;
FLOAT *aa1, *aa2;
FLOAT *b1, *b2;
FLOAT *bb1, *bb2;
FLOAT *cc1, *cc2;
FLOAT a11, a21, a31, a41;
FLOAT a12, a22, a32, a42;
b1 = b;
b2 = b;
lda *= 2;
for (js = 0; js < m; js += 2){
aa1 = a + 0 * lda;
aa2 = a + 1 * lda;
a += 2 * lda + 4;
bb1 = b1 + 0 * m;
bb2 = b1 + 2 * m;
b1 += 4 * m + 4;
cc1 = b2 + 0 * m;
cc2 = b2 + 2 * m;
b2 += 4 * m + 4;
if (m - js >= 2){
a11 = *(aa1 + 0);
a31 = *(aa1 + 2);
a41 = *(aa1 + 3);
a12 = *(aa2 + 2);
*(bb1 + 0) = a11;
*(bb1 + 1) = 0.;
*(bb1 + 2) = a31;
*(bb1 + 3) = a41;
*(bb2 + 0) = a31;
*(bb2 + 1) = -a41;
*(bb2 + 2) = a12;
*(bb2 + 3) = 0.;
aa1 += 4;
aa2 += 4;
bb1 += 4;
bb2 += 4;
cc1 += 4 * m;
cc2 += 4 * m;
is = ((m - js - 2) >> 1);
while (is > 0){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a31 = *(aa1 + 2);
a41 = *(aa1 + 3);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
a32 = *(aa2 + 2);
a42 = *(aa2 + 3);
aa1 += 4;
aa2 += 4;
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb1 + 2) = a31;
*(bb1 + 3) = a41;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(bb2 + 2) = a32;
*(bb2 + 3) = a42;
*(cc1 + 0) = a11;
*(cc1 + 1) = -a21;
*(cc1 + 2) = a12;
*(cc1 + 3) = -a22;
*(cc2 + 0) = a31;
*(cc2 + 1) = -a41;
*(cc2 + 2) = a32;
*(cc2 + 3) = -a42;
bb1 += 4;
bb2 += 4;
cc1 += 4 * m;
cc2 += 4 * m;
is --;
}
if (m & 1){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(cc1 + 0) = a11;
*(cc1 + 1) = -a21;
*(cc1 + 2) = a12;
*(cc1 + 3) = -a22;
}
}
if (m - js == 1){
a11 = *(aa1 + 0);
*(bb1 + 0) = a11;
*(bb1 + 1) = 0.;
}
}
}
static __inline void ZHEMCOPY_U(BLASLONG m, FLOAT *a, BLASLONG lda, FLOAT *b){
BLASLONG is, js;
FLOAT *aa1, *aa2;
FLOAT *b1, *b2;
FLOAT *bb1, *bb2;
FLOAT *cc1, *cc2;
FLOAT a11, a21, a31, a41;
FLOAT a12, a22, a32, a42;
b1 = b;
b2 = b;
lda *= 2;
for (js = 0; js < m; js += 2){
aa1 = a + 0 * lda;
aa2 = a + 1 * lda;
a += 2 * lda;
bb1 = b1 + 0 * m;
bb2 = b1 + 2 * m;
b1 += 4 * m;
cc1 = b2 + 0 * m;
cc2 = b2 + 2 * m;
b2 += 4;
if (m - js >= 2){
for (is = 0; is < js; is += 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a31 = *(aa1 + 2);
a41 = *(aa1 + 3);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
a32 = *(aa2 + 2);
a42 = *(aa2 + 3);
aa1 += 4;
aa2 += 4;
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb1 + 2) = a31;
*(bb1 + 3) = a41;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(bb2 + 2) = a32;
*(bb2 + 3) = a42;
*(cc1 + 0) = a11;
*(cc1 + 1) = -a21;
*(cc1 + 2) = a12;
*(cc1 + 3) = -a22;
*(cc2 + 0) = a31;
*(cc2 + 1) = -a41;
*(cc2 + 2) = a32;
*(cc2 + 3) = -a42;
bb1 += 4;
bb2 += 4;
cc1 += 4 * m;
cc2 += 4 * m;
}
a11 = *(aa1 + 0);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
a32 = *(aa2 + 2);
*(bb1 + 0) = a11;
*(bb1 + 1) = 0.;
*(bb1 + 2) = a12;
*(bb1 + 3) = -a22;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(bb2 + 2) = a32;
*(bb2 + 3) = 0.;
}
if (m - js == 1){
for (is = 0; is < js; is += 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a31 = *(aa1 + 2);
a41 = *(aa1 + 3);
aa1 += 4;
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb1 + 2) = a31;
*(bb1 + 3) = a41;
*(cc1 + 0) = a11;
*(cc1 + 1) = -a21;
*(cc2 + 0) = a31;
*(cc2 + 1) = -a41;
bb1 += 4;
cc1 += 4 * m;
cc2 += 4 * m;
}
a11 = *(aa1 + 0);
*(bb1 + 0) = a11;
*(bb1 + 1) = 0.;
}
}
}
static __inline void ZHEMCOPY_M(BLASLONG m, FLOAT *a, BLASLONG lda, FLOAT *b){
BLASLONG is, js;
FLOAT *aa1, *aa2;
FLOAT *b1, *b2;
FLOAT *bb1, *bb2;
FLOAT *cc1, *cc2;
FLOAT a11, a21, a31, a41;
FLOAT a12, a22, a32, a42;
b1 = b;
b2 = b;
lda *= 2;
for (js = 0; js < m; js += 2){
aa1 = a + 0 * lda;
aa2 = a + 1 * lda;
a += 2 * lda + 4;
bb1 = b1 + 0 * m;
bb2 = b1 + 2 * m;
b1 += 4 * m + 4;
cc1 = b2 + 0 * m;
cc2 = b2 + 2 * m;
b2 += 4 * m + 4;
if (m - js >= 2){
a11 = *(aa1 + 0);
a31 = *(aa1 + 2);
a41 = *(aa1 + 3);
a12 = *(aa2 + 2);
*(bb1 + 0) = a11;
*(bb1 + 1) = 0.;
*(bb1 + 2) = a31;
*(bb1 + 3) = -a41;
*(bb2 + 0) = a31;
*(bb2 + 1) = a41;
*(bb2 + 2) = a12;
*(bb2 + 3) = 0.;
aa1 += 4;
aa2 += 4;
bb1 += 4;
bb2 += 4;
cc1 += 4 * m;
cc2 += 4 * m;
is = ((m - js - 2) >> 1);
while (is > 0){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a31 = *(aa1 + 2);
a41 = *(aa1 + 3);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
a32 = *(aa2 + 2);
a42 = *(aa2 + 3);
aa1 += 4;
aa2 += 4;
*(bb1 + 0) = a11;
*(bb1 + 1) = -a21;
*(bb1 + 2) = a31;
*(bb1 + 3) = -a41;
*(bb2 + 0) = a12;
*(bb2 + 1) = -a22;
*(bb2 + 2) = a32;
*(bb2 + 3) = -a42;
*(cc1 + 0) = a11;
*(cc1 + 1) = a21;
*(cc1 + 2) = a12;
*(cc1 + 3) = a22;
*(cc2 + 0) = a31;
*(cc2 + 1) = a41;
*(cc2 + 2) = a32;
*(cc2 + 3) = a42;
bb1 += 4;
bb2 += 4;
cc1 += 4 * m;
cc2 += 4 * m;
is --;
}
if (m & 1){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
*(bb1 + 0) = a11;
*(bb1 + 1) = -a21;
*(bb2 + 0) = a12;
*(bb2 + 1) = -a22;
*(cc1 + 0) = a11;
*(cc1 + 1) = a21;
*(cc1 + 2) = a12;
*(cc1 + 3) = a22;
}
}
if (m - js == 1){
a11 = *(aa1 + 0);
*(bb1 + 0) = a11;
*(bb1 + 1) = 0.;
}
}
}
static __inline void ZHEMCOPY_V(BLASLONG m, FLOAT *a, BLASLONG lda, FLOAT *b){
BLASLONG is, js;
FLOAT *aa1, *aa2;
FLOAT *b1, *b2;
FLOAT *bb1, *bb2;
FLOAT *cc1, *cc2;
FLOAT a11, a21, a31, a41;
FLOAT a12, a22, a32, a42;
b1 = b;
b2 = b;
lda *= 2;
for (js = 0; js < m; js += 2){
aa1 = a + 0 * lda;
aa2 = a + 1 * lda;
a += 2 * lda;
bb1 = b1 + 0 * m;
bb2 = b1 + 2 * m;
b1 += 4 * m;
cc1 = b2 + 0 * m;
cc2 = b2 + 2 * m;
b2 += 4;
if (m - js >= 2){
for (is = 0; is < js; is += 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a31 = *(aa1 + 2);
a41 = *(aa1 + 3);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
a32 = *(aa2 + 2);
a42 = *(aa2 + 3);
aa1 += 4;
aa2 += 4;
*(bb1 + 0) = a11;
*(bb1 + 1) = -a21;
*(bb1 + 2) = a31;
*(bb1 + 3) = -a41;
*(bb2 + 0) = a12;
*(bb2 + 1) = -a22;
*(bb2 + 2) = a32;
*(bb2 + 3) = -a42;
*(cc1 + 0) = a11;
*(cc1 + 1) = a21;
*(cc1 + 2) = a12;
*(cc1 + 3) = a22;
*(cc2 + 0) = a31;
*(cc2 + 1) = a41;
*(cc2 + 2) = a32;
*(cc2 + 3) = a42;
bb1 += 4;
bb2 += 4;
cc1 += 4 * m;
cc2 += 4 * m;
}
a11 = *(aa1 + 0);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
a32 = *(aa2 + 2);
*(bb1 + 0) = a11;
*(bb1 + 1) = 0.;
*(bb1 + 2) = a12;
*(bb1 + 3) = a22;
*(bb2 + 0) = a12;
*(bb2 + 1) = -a22;
*(bb2 + 2) = a32;
*(bb2 + 3) = 0.;
}
if (m - js == 1){
for (is = 0; is < js; is += 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a31 = *(aa1 + 2);
a41 = *(aa1 + 3);
aa1 += 4;
*(bb1 + 0) = a11;
*(bb1 + 1) = -a21;
*(bb1 + 2) = a31;
*(bb1 + 3) = -a41;
*(cc1 + 0) = a11;
*(cc1 + 1) = a21;
*(cc2 + 0) = a31;
*(cc2 + 1) = a41;
bb1 += 4;
cc1 += 4 * m;
cc2 += 4 * m;
}
a11 = *(aa1 + 0);
*(bb1 + 0) = a11;
*(bb1 + 1) = 0.;
}
}
}
static __inline void TRMCOPY_NL(BLASLONG m, FLOAT *a, BLASLONG lda, FLOAT *b){
BLASLONG is, js;
FLOAT *aa1, *aa2;
FLOAT *b1, *b2;
FLOAT *bb1, *bb2;
FLOAT *cc1, *cc2;
FLOAT a11, a12;
FLOAT a21, a22;
b1 = b;
b2 = b;
for (js = 0; js < m; js += 2){
aa1 = a + 0 * lda;
aa2 = a + 1 * lda;
a += 2 * lda + 2;
bb1 = b1 + 0 * m;
bb2 = b1 + 1 * m;
b1 += 2 * m + 2;
cc1 = b2 + 0 * m;
cc2 = b2 + 1 * m;
b2 += 2 * m + 2;
if (m - js >= 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a22 = *(aa2 + 1);
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb2 + 0) = a21;
*(bb2 + 1) = a22;
aa1 += 2;
aa2 += 2;
bb1 += 2;
bb2 += 2;
cc1 += 2 * m;
cc2 += 2 * m;
is = ((m - js - 2) >> 1);
while (is > 0){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
aa1 += 2;
aa2 += 2;
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(cc1 + 0) = a11;
*(cc1 + 1) = a12;
*(cc2 + 0) = a21;
*(cc2 + 1) = a22;
bb1 += 2;
bb2 += 2;
cc1 += 2 * m;
cc2 += 2 * m;
is --;
}
is = ((m - js - 2) & 1);
if (is == 1){
a11 = *(aa1 + 0);
a12 = *(aa2 + 0);
*(bb1 + 0) = a11;
*(bb2 + 0) = a12;
*(cc1 + 0) = a11;
*(cc1 + 1) = a12;
}
}
if (m - js == 1){
a11 = *(aa1 + 0);
*(bb1 + 0) = a11;
}
}
}
static __inline void TRMCOPY_TL(BLASLONG m, FLOAT *a, BLASLONG lda, FLOAT *b){
BLASLONG is, js;
FLOAT *aa1, *aa2;
FLOAT *b1, *b2;
FLOAT *bb1, *bb2;
FLOAT *cc1, *cc2;
FLOAT a11, a12;
FLOAT a21, a22;
b1 = b;
b2 = b;
for (js = 0; js < m; js += 2){
aa1 = a + 0 * lda;
aa2 = a + 1 * lda;
a += 2 * lda + 2;
bb1 = b1 + 0 * m;
bb2 = b1 + 1 * m;
b1 += 2 * m + 2;
cc1 = b2 + 0 * m;
cc2 = b2 + 1 * m;
b2 += 2 * m + 2;
if (m - js >= 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a22 = *(aa2 + 1);
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb2 + 0) = a21;
*(bb2 + 1) = a22;
aa1 += 2;
aa2 += 2;
bb1 += 2;
bb2 += 2;
cc1 += 2 * m;
cc2 += 2 * m;
is = ((m - js - 2) >> 1);
while (is > 0){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
aa1 += 2;
aa2 += 2;
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(cc1 + 0) = a11;
*(cc1 + 1) = a12;
*(cc2 + 0) = a21;
*(cc2 + 1) = a22;
bb1 += 2;
bb2 += 2;
cc1 += 2 * m;
cc2 += 2 * m;
is --;
}
is = ((m - js - 2) & 1);
if (is == 1){
a11 = *(aa1 + 0);
a12 = *(aa2 + 0);
*(bb1 + 0) = a11;
*(bb2 + 0) = a12;
*(cc1 + 0) = a11;
*(cc1 + 1) = a12;
}
}
if (m - js == 1){
a11 = *(aa1 + 0);
*(bb1 + 0) = a11;
}
}
}
static __inline void TRMCOPY_NU(BLASLONG m, FLOAT *a, BLASLONG lda, FLOAT *b){
BLASLONG is, js;
FLOAT *aa1, *aa2;
FLOAT *b1, *b2;
FLOAT *bb1, *bb2;
FLOAT *cc1, *cc2;
FLOAT a11, a12;
FLOAT a21, a22;
b1 = b;
b2 = b;
for (js = 0; js < m; js += 2){
aa1 = a + 0 * lda;
aa2 = a + 1 * lda;
a += 2 * lda;
bb1 = b1 + 0 * m;
bb2 = b1 + 1 * m;
b1 += 2 * m;
cc1 = b2 + 0 * m;
cc2 = b2 + 1 * m;
b2 += 2;
if (m - js >= 2){
for (is = 0; is < js; is += 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
aa1 += 2;
aa2 += 2;
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(cc1 + 0) = a11;
*(cc1 + 1) = a12;
*(cc2 + 0) = a21;
*(cc2 + 1) = a22;
bb1 += 2;
bb2 += 2;
cc1 += 2 * m;
cc2 += 2 * m;
}
a11 = *(aa1 + 0);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
*(bb1 + 0) = a11;
*(bb1 + 1) = a12;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
}
if (m - js == 1){
for (is = 0; is < js; is += 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
aa1 += 2;
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(cc1 + 0) = a11;
*(cc2 + 0) = a21;
bb1 += 2;
cc1 += 2 * m;
cc2 += 2 * m;
}
a11 = *(aa1 + 0);
*(bb1 + 0) = a11;
}
}
}
static __inline void TRMCOPY_TU(BLASLONG m, FLOAT *a, BLASLONG lda, FLOAT *b){
BLASLONG is, js;
FLOAT *aa1, *aa2;
FLOAT *b1, *b2;
FLOAT *bb1, *bb2;
FLOAT *cc1, *cc2;
FLOAT a11, a12;
FLOAT a21, a22;
b1 = b;
b2 = b;
for (js = 0; js < m; js += 2){
aa1 = a + 0 * lda;
aa2 = a + 1 * lda;
a += 2 * lda;
bb1 = b1 + 0 * m;
bb2 = b1 + 1 * m;
b1 += 2 * m;
cc1 = b2 + 0 * m;
cc2 = b2 + 1 * m;
b2 += 2;
if (m - js >= 2){
for (is = 0; is < js; is += 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
aa1 += 2;
aa2 += 2;
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(cc1 + 0) = a11;
*(cc1 + 1) = a12;
*(cc2 + 0) = a21;
*(cc2 + 1) = a22;
bb1 += 2;
bb2 += 2;
cc1 += 2 * m;
cc2 += 2 * m;
}
a11 = *(aa1 + 0);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
*(bb1 + 0) = a11;
*(bb1 + 1) = a12;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
}
if (m - js == 1){
for (is = 0; is < js; is += 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
aa1 += 2;
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(cc1 + 0) = a11;
*(cc2 + 0) = a21;
bb1 += 2;
cc1 += 2 * m;
cc2 += 2 * m;
}
a11 = *(aa1 + 0);
*(bb1 + 0) = a11;
}
}
}
static __inline void ZTRMCOPY_NL(BLASLONG m, FLOAT *a, BLASLONG lda, FLOAT *b){
BLASLONG is, js;
FLOAT *aa1, *aa2;
FLOAT *b1, *b2;
FLOAT *bb1, *bb2;
FLOAT *cc1, *cc2;
FLOAT a11, a21, a31, a41;
FLOAT a12, a22, a32, a42;
b1 = b;
b2 = b;
lda *= 2;
for (js = 0; js < m; js += 2){
aa1 = a + 0 * lda;
aa2 = a + 1 * lda;
a += 2 * lda + 4;
bb1 = b1 + 0 * m;
bb2 = b1 + 2 * m;
b1 += 4 * m + 4;
cc1 = b2 + 0 * m;
cc2 = b2 + 2 * m;
b2 += 4 * m + 4;
if (m - js >= 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a31 = *(aa1 + 2);
a41 = *(aa1 + 3);
a12 = *(aa2 + 2);
a22 = *(aa2 + 3);
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb1 + 2) = a31;
*(bb1 + 3) = a41;
*(bb2 + 0) = a31;
*(bb2 + 1) = a41;
*(bb2 + 2) = a12;
*(bb2 + 3) = a22;
aa1 += 4;
aa2 += 4;
bb1 += 4;
bb2 += 4;
cc1 += 4 * m;
cc2 += 4 * m;
is = ((m - js - 2) >> 1);
while (is > 0){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a31 = *(aa1 + 2);
a41 = *(aa1 + 3);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
a32 = *(aa2 + 2);
a42 = *(aa2 + 3);
aa1 += 4;
aa2 += 4;
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb1 + 2) = a31;
*(bb1 + 3) = a41;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(bb2 + 2) = a32;
*(bb2 + 3) = a42;
*(cc1 + 0) = a11;
*(cc1 + 1) = a21;
*(cc1 + 2) = a12;
*(cc1 + 3) = a22;
*(cc2 + 0) = a31;
*(cc2 + 1) = a41;
*(cc2 + 2) = a32;
*(cc2 + 3) = a42;
bb1 += 4;
bb2 += 4;
cc1 += 4 * m;
cc2 += 4 * m;
is --;
}
if (m & 1){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(cc1 + 0) = a11;
*(cc1 + 1) = a21;
*(cc1 + 2) = a12;
*(cc1 + 3) = a22;
}
}
if (m - js == 1){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
}
}
}
static __inline void ZTRMCOPY_TL(BLASLONG m, FLOAT *a, BLASLONG lda, FLOAT *b){
BLASLONG is, js;
FLOAT *aa1, *aa2;
FLOAT *b1, *b2;
FLOAT *bb1, *bb2;
FLOAT *cc1, *cc2;
FLOAT a11, a21, a31, a41;
FLOAT a12, a22, a32, a42;
b1 = b;
b2 = b;
lda *= 2;
for (js = 0; js < m; js += 2){
aa1 = a + 0 * lda;
aa2 = a + 1 * lda;
a += 2 * lda + 4;
bb1 = b1 + 0 * m;
bb2 = b1 + 2 * m;
b1 += 4 * m + 4;
cc1 = b2 + 0 * m;
cc2 = b2 + 2 * m;
b2 += 4 * m + 4;
if (m - js >= 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a31 = *(aa1 + 2);
a41 = *(aa1 + 3);
a12 = *(aa2 + 2);
a22 = *(aa2 + 3);
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb1 + 2) = a31;
*(bb1 + 3) = a41;
*(bb2 + 0) = a31;
*(bb2 + 1) = a41;
*(bb2 + 2) = a12;
*(bb2 + 3) = a22;
aa1 += 4;
aa2 += 4;
bb1 += 4;
bb2 += 4;
cc1 += 4 * m;
cc2 += 4 * m;
is = ((m - js - 2) >> 1);
while (is > 0){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a31 = *(aa1 + 2);
a41 = *(aa1 + 3);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
a32 = *(aa2 + 2);
a42 = *(aa2 + 3);
aa1 += 4;
aa2 += 4;
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb1 + 2) = a31;
*(bb1 + 3) = a41;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(bb2 + 2) = a32;
*(bb2 + 3) = a42;
*(cc1 + 0) = a11;
*(cc1 + 1) = a21;
*(cc1 + 2) = a12;
*(cc1 + 3) = a22;
*(cc2 + 0) = a31;
*(cc2 + 1) = a41;
*(cc2 + 2) = a32;
*(cc2 + 3) = a42;
bb1 += 4;
bb2 += 4;
cc1 += 4 * m;
cc2 += 4 * m;
is --;
}
if (m & 1){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(cc1 + 0) = a11;
*(cc1 + 1) = a21;
*(cc1 + 2) = a12;
*(cc1 + 3) = a22;
}
}
if (m - js == 1){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
}
}
}
static __inline void ZTRMCOPY_NU(BLASLONG m, FLOAT *a, BLASLONG lda, FLOAT *b){
BLASLONG is, js;
FLOAT *aa1, *aa2;
FLOAT *b1, *b2;
FLOAT *bb1, *bb2;
FLOAT *cc1, *cc2;
FLOAT a11, a21, a31, a41;
FLOAT a12, a22, a32, a42;
b1 = b;
b2 = b;
lda *= 2;
for (js = 0; js < m; js += 2){
aa1 = a + 0 * lda;
aa2 = a + 1 * lda;
a += 2 * lda;
bb1 = b1 + 0 * m;
bb2 = b1 + 2 * m;
b1 += 4 * m;
cc1 = b2 + 0 * m;
cc2 = b2 + 2 * m;
b2 += 4;
if (m - js >= 2){
for (is = 0; is < js; is += 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a31 = *(aa1 + 2);
a41 = *(aa1 + 3);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
a32 = *(aa2 + 2);
a42 = *(aa2 + 3);
aa1 += 4;
aa2 += 4;
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb1 + 2) = a31;
*(bb1 + 3) = a41;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(bb2 + 2) = a32;
*(bb2 + 3) = a42;
*(cc1 + 0) = a11;
*(cc1 + 1) = a21;
*(cc1 + 2) = a12;
*(cc1 + 3) = a22;
*(cc2 + 0) = a31;
*(cc2 + 1) = a41;
*(cc2 + 2) = a32;
*(cc2 + 3) = a42;
bb1 += 4;
bb2 += 4;
cc1 += 4 * m;
cc2 += 4 * m;
}
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
a32 = *(aa2 + 2);
a42 = *(aa2 + 3);
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb1 + 2) = a12;
*(bb1 + 3) = a22;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(bb2 + 2) = a32;
*(bb2 + 3) = a42;
}
if (m - js == 1){
for (is = 0; is < js; is += 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a31 = *(aa1 + 2);
a41 = *(aa1 + 3);
aa1 += 4;
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb1 + 2) = a31;
*(bb1 + 3) = a41;
*(cc1 + 0) = a11;
*(cc1 + 1) = a21;
*(cc2 + 0) = a31;
*(cc2 + 1) = a41;
bb1 += 4;
cc1 += 4 * m;
cc2 += 4 * m;
}
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
}
}
}
static __inline void ZTRMCOPY_TU(BLASLONG m, FLOAT *a, BLASLONG lda, FLOAT *b){
BLASLONG is, js;
FLOAT *aa1, *aa2;
FLOAT *b1, *b2;
FLOAT *bb1, *bb2;
FLOAT *cc1, *cc2;
FLOAT a11, a21, a31, a41;
FLOAT a12, a22, a32, a42;
b1 = b;
b2 = b;
lda *= 2;
for (js = 0; js < m; js += 2){
aa1 = a + 0 * lda;
aa2 = a + 1 * lda;
a += 2 * lda;
bb1 = b1 + 0 * m;
bb2 = b1 + 2 * m;
b1 += 4 * m;
cc1 = b2 + 0 * m;
cc2 = b2 + 2 * m;
b2 += 4;
if (m - js >= 2){
for (is = 0; is < js; is += 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a31 = *(aa1 + 2);
a41 = *(aa1 + 3);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
a32 = *(aa2 + 2);
a42 = *(aa2 + 3);
aa1 += 4;
aa2 += 4;
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb1 + 2) = a31;
*(bb1 + 3) = a41;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(bb2 + 2) = a32;
*(bb2 + 3) = a42;
*(cc1 + 0) = a11;
*(cc1 + 1) = a21;
*(cc1 + 2) = a12;
*(cc1 + 3) = a22;
*(cc2 + 0) = a31;
*(cc2 + 1) = a41;
*(cc2 + 2) = a32;
*(cc2 + 3) = a42;
bb1 += 4;
bb2 += 4;
cc1 += 4 * m;
cc2 += 4 * m;
}
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a12 = *(aa2 + 0);
a22 = *(aa2 + 1);
a32 = *(aa2 + 2);
a42 = *(aa2 + 3);
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb1 + 2) = a12;
*(bb1 + 3) = a22;
*(bb2 + 0) = a12;
*(bb2 + 1) = a22;
*(bb2 + 2) = a32;
*(bb2 + 3) = a42;
}
if (m - js == 1){
for (is = 0; is < js; is += 2){
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
a31 = *(aa1 + 2);
a41 = *(aa1 + 3);
aa1 += 4;
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
*(bb1 + 2) = a31;
*(bb1 + 3) = a41;
*(cc1 + 0) = a11;
*(cc1 + 1) = a21;
*(cc2 + 0) = a31;
*(cc2 + 1) = a41;
bb1 += 4;
cc1 += 4 * m;
cc2 += 4 * m;
}
a11 = *(aa1 + 0);
a21 = *(aa1 + 1);
*(bb1 + 0) = a11;
*(bb1 + 1) = a21;
}
}
}
#endif
#endif
马建仓 AI 助手
尝试更多
代码解读
代码找茬
代码优化
1
https://gitee.com/flice/OpenBLAS.git
git@gitee.com:flice/OpenBLAS.git
flice
OpenBLAS
OpenBLAS
master

搜索帮助

0d507c66 1850385 C8b1a773 1850385