From patchwork Thu Apr 18 09:42:22 2024 Content-Type: text/plain; charset="utf-8" MIME-Version: 1.0 Content-Transfer-Encoding: 7bit X-Patchwork-Submitter: pengxu X-Patchwork-Id: 48149 Delivered-To: ffmpegpatchwork2@gmail.com Received: by 2002:a05:6a20:ce4e:b0:1a9:af23:56c1 with SMTP id id14csp1570834pzb; Thu, 18 Apr 2024 02:42:50 -0700 (PDT) X-Forwarded-Encrypted: i=2; AJvYcCVp8Aejoknpu0uO8F8U0vcsIuMHiksh/JVVTAixSNT/A3Vjl+jfJlv7EgImp+089ILksZqN4vcOc+OtFVW//QBz/OzsqXcJjSeFWQ== X-Google-Smtp-Source: AGHT+IHYM1XtusRhxsnC4UtZOVtp+QEqJEa38kXHbcSxYnDkyErYXJPYOOdt9RpC+zmOf5nQmogS X-Received: by 2002:a17:906:1389:b0:a4e:410e:9525 with SMTP id f9-20020a170906138900b00a4e410e9525mr1298656ejc.30.1713433370423; Thu, 18 Apr 2024 02:42:50 -0700 (PDT) ARC-Seal: i=1; a=rsa-sha256; t=1713433370; cv=none; d=google.com; s=arc-20160816; b=fMfNeTTMDikZUC8XIvy6ciEbp+5wpdGk4B39X63qusjwHoMe0b9Dg5DWR8RdBy/1TH Ij8uhN50D9Mq7x9UPFSxMi5J2SXqYnOJ4oclT7or51653qk9bx2MyhwHVBkLZgGBSB7T QlN/Z4XYSbrP9jPuwAQRpzRNoOslLz0tn01g8GAdwCErKhWM0Z1gXoQKdBYuUGkW9AL8 nn6L12D24hvy8iQCc79JulePXcCamZfnhWGcyvsLB2ZmDH7MGICUlFh/bt1O9ZMbzoft JdbOK9Iw5TlfWgUGeSCPngOcrOv4WuAGpFD+4DlOJx6MwesLvwJQLNEfUw5xT4mJyQ1L NzzA== ARC-Message-Signature: i=1; a=rsa-sha256; c=relaxed/relaxed; d=google.com; s=arc-20160816; h=sender:errors-to:content-transfer-encoding:reply-to:list-subscribe :list-help:list-post:list-archive:list-unsubscribe:list-id :precedence:subject:mime-version:references:in-reply-to:message-id :date:to:from:delivered-to; bh=TONJP+MU6uukXZrtcSNlQIwsDOj/SOKhkSUosV1vuKo=; fh=YOA8vD9MJZuwZ71F/05pj6KdCjf6jQRmzLS+CATXUQk=; b=xYhejmOfXF9a35ybTThAMzWl+t3/DRqbC5fFbv7geGPNrkHEE9pnKXYtBWxiMe5TLP RuHn4/xl4j4VdxfTmOPd7ORgJvq9tY3Kn5MauslEQ47ox+386Lo/nHodTuJksFnDFmCI akQRn2i47SYTuxfoOGnUwP9HESLehMG9T4liGt9MRI6o44+k1WpVnTBZ8WHERWeEK8ig itxP+mBKZXJaBjLLfe1RSo08yIBG2SD686TC9d+9V92bP5Uf+WtWWMMGT3fi9htKrmw6 n8rY6P16Y+LEYmzjDP3c7lmcAqoIYLUiL/tWK4uBsLU4gAqGAWWE5yw3aO9NmpeZVVxL m+jw==; dara=google.com ARC-Authentication-Results: i=1; mx.google.com; spf=pass (google.com: domain of ffmpeg-devel-bounces@ffmpeg.org designates 79.124.17.100 as permitted sender) smtp.mailfrom=ffmpeg-devel-bounces@ffmpeg.org Return-Path: Received: from ffbox0-bg.mplayerhq.hu (ffbox0-bg.ffmpeg.org. [79.124.17.100]) by mx.google.com with ESMTP id k2-20020a170906578200b00a51a11372a2si643741ejq.301.2024.04.18.02.42.50; Thu, 18 Apr 2024 02:42:50 -0700 (PDT) Received-SPF: pass (google.com: domain of ffmpeg-devel-bounces@ffmpeg.org designates 79.124.17.100 as permitted sender) client-ip=79.124.17.100; Authentication-Results: mx.google.com; spf=pass (google.com: domain of ffmpeg-devel-bounces@ffmpeg.org designates 79.124.17.100 as permitted sender) smtp.mailfrom=ffmpeg-devel-bounces@ffmpeg.org Received: from [127.0.1.1] (localhost [127.0.0.1]) by ffbox0-bg.mplayerhq.hu (Postfix) with ESMTP id 7A9AC68D441; Thu, 18 Apr 2024 12:42:39 +0300 (EEST) X-Original-To: ffmpeg-devel@ffmpeg.org Delivered-To: ffmpeg-devel@ffmpeg.org Received: from mail.loongson.cn (mail.loongson.cn [114.242.206.163]) by ffbox0-bg.mplayerhq.hu (Postfix) with ESMTP id 7290668D42F for ; Thu, 18 Apr 2024 12:42:32 +0300 (EEST) Received: from loongson.cn (unknown [36.33.26.33]) by gateway (Coremail) with SMTP id _____8BxXbsF6yBmZzUpAA--.12322S3; Thu, 18 Apr 2024 17:42:29 +0800 (CST) Received: from localhost (unknown [36.33.26.33]) by localhost.localdomain (Coremail) with SMTP id AQAAf8AxTRMC6yBm8Ph+AA--.45196S3; Thu, 18 Apr 2024 17:42:27 +0800 (CST) From: pengxu To: ffmpeg-devel@ffmpeg.org Date: Thu, 18 Apr 2024 17:42:22 +0800 Message-Id: <20240418094223.21492-2-pengxu@loongson.cn> X-Mailer: git-send-email 2.20.1 In-Reply-To: <20240418094223.21492-1-pengxu@loongson.cn> References: <20240418094223.21492-1-pengxu@loongson.cn> MIME-Version: 1.0 X-CM-TRANSID: AQAAf8AxTRMC6yBm8Ph+AA--.45196S3 X-CM-SenderInfo: pshqw53x6o00pqjv00gofq/ X-Coremail-Antispam: 1Uk129KBj93XoW3uFWxZr4xGFy7uryrJrWfXrc_yoWDXF15pF W3uw1jqw4kKFyS9F4kZ3s2vr1rXr97CF1SgF98WFW8ZFW8Gr4UXrZrtF9xCFyxX3yUAF4Y 9ayrKa4SyFyrAwcCm3ZEXasCq-sJn29KB7ZKAUJUUUUU529EdanIXcx71UUUUU7KY7ZEXa sCq-sGcSsGvfJ3Ic02F40EFcxC0VAKzVAqx4xG6I80ebIjqfuFe4nvWSU5nxnvy29KBjDU 0xBIdaVrnRJUUUkFb4IE77IF4wAFF20E14v26r1j6r4UM7CY07I20VC2zVCF04k26cxKx2 IYs7xG6rWj6s0DM7CIcVAFz4kK6r1j6r18M28lY4IEw2IIxxk0rwA2F7IY1VAKz4vEj48v e4kI8wA2z4x0Y4vE2Ix0cI8IcVAFwI0_JFI_Gr1l84ACjcxK6xIIjxv20xvEc7CjxVAFwI 0_Jr0_Gr1l84ACjcxK6I8E87Iv67AKxVW8Jr0_Cr1UM28EF7xvwVC2z280aVCY1x0267AK xVW8Jr0_Cr1UM2AIxVAIcxkEcVAq07x20xvEncxIr21l57IF6xkI12xvs2x26I8E6xACxx 1l5I8CrVACY4xI64kE6c02F40Ex7xfMcIj6xIIjxv20xvE14v26r1Y6r17McIj6I8E87Iv 67AKxVW8JVWxJwAm72CE4IkC6x0Yz7v_Jr0_Gr1lF7xvr2IYc2Ij64vIr41l42xK82IYc2 Ij64vIr41l4I8I3I0E4IkC6x0Yz7v_Jr0_Gr1lx2IqxVAqx4xG67AKxVWUJVWUGwC20s02 6x8GjcxK67AKxVWUGVWUWwC2zVAF1VAY17CE14v26r1j6r15MIIYrxkI7VAKI48JMIIF0x vE2Ix0cI8IcVAFwI0_Jr0_JF4lIxAIcVC0I7IYx2IY6xkF7I0E14v26r1j6r4UMIIF0xvE 42xK8VAvwI8IcIk0rVWUJVWUCwCI42IY6I8E87Iv67AKxVWUJVW8JwCI42IY6I8E87Iv6x kF7I0E14v26r1j6r4UYxBIdaVFxhVjvjDU0xZFpf9x07jb_-PUUUUU= Subject: [FFmpeg-devel] [PATCH v3 1/2] avutil/loongarch: add LSX optimization for aac audio decode X-BeenThere: ffmpeg-devel@ffmpeg.org X-Mailman-Version: 2.1.29 Precedence: list List-Id: FFmpeg development discussions and patches List-Unsubscribe: , List-Archive: List-Post: List-Help: List-Subscribe: , Reply-To: FFmpeg development discussions and patches Errors-To: ffmpeg-devel-bounces@ffmpeg.org Sender: "ffmpeg-devel" X-TUID: cAPb7nVadVXU Add functions: vector_fmul_window_lsx butterflies_float_lsx vector_fmul_scalar_lsx ./ffmpeg -i ../../1.aac -f null - before:482x after:523x --- libavutil/float_dsp.c | 2 + libavutil/float_dsp.h | 1 + libavutil/loongarch/Makefile | 5 +- libavutil/loongarch/float_dsp.S | 287 ++++++++++++++++++ libavutil/loongarch/float_dsp.h | 32 ++ .../loongarch/float_dsp_init_loongarch.c | 35 +++ 6 files changed, 361 insertions(+), 1 deletion(-) create mode 100644 libavutil/loongarch/float_dsp.S create mode 100644 libavutil/loongarch/float_dsp.h create mode 100644 libavutil/loongarch/float_dsp_init_loongarch.c diff --git a/libavutil/float_dsp.c b/libavutil/float_dsp.c index e9fb023466..7128ff3f96 100644 --- a/libavutil/float_dsp.c +++ b/libavutil/float_dsp.c @@ -162,6 +162,8 @@ av_cold AVFloatDSPContext *avpriv_float_dsp_alloc(int bit_exact) ff_float_dsp_init_x86(fdsp); #elif ARCH_MIPS ff_float_dsp_init_mips(fdsp); +#elif ARCH_LOONGARCH64 + ff_float_dsp_init_loongarch(fdsp); #endif return fdsp; } diff --git a/libavutil/float_dsp.h b/libavutil/float_dsp.h index 342a8715c5..679a930eab 100644 --- a/libavutil/float_dsp.h +++ b/libavutil/float_dsp.h @@ -206,6 +206,7 @@ void ff_float_dsp_init_ppc(AVFloatDSPContext *fdsp, int strict); void ff_float_dsp_init_riscv(AVFloatDSPContext *fdsp); void ff_float_dsp_init_x86(AVFloatDSPContext *fdsp); void ff_float_dsp_init_mips(AVFloatDSPContext *fdsp); +void ff_float_dsp_init_loongarch(AVFloatDSPContext *fdsp); /** * Allocate a float DSP context. diff --git a/libavutil/loongarch/Makefile b/libavutil/loongarch/Makefile index 2addd9351c..ae710f0515 100644 --- a/libavutil/loongarch/Makefile +++ b/libavutil/loongarch/Makefile @@ -1 +1,4 @@ -OBJS += loongarch/cpu.o +OBJS += loongarch/cpu.o \ + loongarch/float_dsp_init_loongarch.o + +LSX-OBJS += loongarch/float_dsp.o diff --git a/libavutil/loongarch/float_dsp.S b/libavutil/loongarch/float_dsp.S new file mode 100644 index 0000000000..5073c8424f --- /dev/null +++ b/libavutil/loongarch/float_dsp.S @@ -0,0 +1,287 @@ +/* + * Loongarch LASX/LSX optimizeds dsp + * + * Copyright (c) 2024 Loongson Technology Corporation Limited + * Contributed by PengXu + * + * This file is part of FFmpeg. + * + * FFmpeg is free software; you can redistribute it and/or + * modify it under the terms of the GNU Lesser General Public + * License as published by the Free Software Foundation; either + * version 2.1 of the License, or (at your option) any later version. + * + * FFmpeg is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + * Lesser General Public License for more details. + * + * You should have received a copy of the GNU Lesser General Public + * License along with FFmpeg; if not, write to the Free Software + * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA + */ + +#include "libavcodec/loongarch/loongson_asm.S" + + +/* void vector_fmul_window(float *dst, const float *src0, + const float *src1, const float *win, int len) */ +function vector_fmul_window_lsx + addi.d sp, sp, -8 + st.d $r23, sp, 0 + + move t4, a0 + move t5, a1 + move t6, a2 + move t7, a3 + move t8, a4 + slli.d t8, t8, 2 + + add.d t4, t4, t8 + add.d t7, t7, t8 + add.d t5, t5, t8 + + add.d a6, $r0, t8 + addi.d a7, t8, -4 + + move a5, $r0 + srai.d t0, a4, 2 + beq a5, t0, .VFW02 + +.VFW01: + sub.d t1, t5, a6 + addi.d t2, a7, -12 + vld vr1, t1, 0x00 //s0 + vldx vr2, a2, t2 //s1 + + sub.d t1, t7, a6 + vld vr3, t1, 0x00 //wi + vldx vr4, t7, t2 //wj + + vpermi.w vr2, vr2, 0x1b + vpermi.w vr4, vr4, 0x1b + + vfmul.s vr5, vr2, vr3 + vfmsub.s vr5, vr1, vr4, vr5 //dsti + + vfmul.s vr6, vr2, vr4 + vfmadd.s vr6, vr1, vr3, vr6 //dstj + + vpermi.w vr6, vr6, 0x1b + + sub.d t1, t4, a6 + vst vr5, t1, 0x00 + vstx vr6, t4, t2 + + addi.d a6, a6, -16 + addi.d a7, a7, -16 + + addi.d a5, a5, 1 + blt a5, t0, .VFW01 + +.VFW02: + andi t0, a4, 2 + beq $r0, t0, .VFW03 + + sub.d t0, t5, a6 + addi.d t1, a7, -4 + add.d t1, t5, t1 + + sub.d t2, t7, a6 + addi.d t3, a7, -4 + add.d t3, t7, t3 + + fld.s f0, t0, 0x00 //s0 + fld.s f1, t0, 0x04 + + fld.s f2, t1, 0x04 //s1 + fld.s f3, t1, 0x00 + + fld.s f4, t2, 0x00 //wi + fld.s f5, t2, 0x04 + + fld.s f6, t3, 0x04 //wj + fld.s f7, t3, 0x00 + + fmul.s f8, f2, f4 + fmsub.s f8, f0, f6, f8 //dsti + fmul.s f9, f3, f5 + fmsub.s f9, f1, f7, f9 + + fmul.s f10, f2, f6 + fmadd.s f10, f0, f4, f10 //dstj + fmul.s f11, f3, f7 + fmadd.s f11, f1, f5, f11 + + sub.d t2, t4, a6 + add.d t3, t4, a7 + addi.d t3, t3, -4 + + fst.s f8, t2, 0x00 + fst.s f9, t2, 0x04 + fst.s f10, t3, 0x04 + fst.s f11, t3, 0x00 + + addi.d a6, a6, -2 + addi.d a7, a7, -2 + +.VFW03: + andi t0, a4, 1 + beq $r0, t0, .VFW04 + + sub.d t0, t5, a6 + + fldx.s f0, t5, t0 //s0 + fldx.s f2, t6, a7 //s1 + fldx.s f4, t7, t0 //wi + fldx.s f6, t7, a7 //wj + + fmul.s f8, f2, f4 + fmsub.s f8, f0, f6, f8 //dsti + + fmul.s f10, f2, f6 + fmadd.s f10, f0, f4, f10 //dstj + + sub.d t0, t4, a6 + + fst.s f8, t0, 0x00 + fstx.s f10, t4, a7 + + addi.d a6, a6, -1 + addi.d a7, a7, -1 + +.VFW04: + ld.d $r23, sp, 0 + addi.d sp, sp, 8 + +endfunc + + +/* void butterflies_float(float *restrict v1, float *restrict v2, + int len) */ +function butterflies_float_lsx + move a6, $r0 + move a7, $r0 + + move t4, a0 + move t5, a1 + move t6, a2 + + srai.d t0, t6, 2 + beq a6, t0, .BFL02 + +.BFL01: + vldx vr0, t4, a7 + vldx vr1, t5, a7 + + vfsub.s vr3, vr0, vr1 + vfadd.s vr4, vr0, vr1 + + vstx vr4, t4, a7 + vstx vr3, t5, a7 + + addi.d a7, a7, 16 + addi.d a6, a6, 1 + blt a6, t0, .BFL01 + +.BFL02: + andi t0, t6, 2 + beq $r0, t0, .BFL03 + + add.d t1, t4, a7 + add.d t2, t5, a7 + + fld.s f0, t1, 0x00 + fld.s f1, t1, 0x04 + fld.s f2, t2, 0x00 + fld.s f3, t2, 0x04 + + fsub.s f4, f0, f2 + fsub.s f5, f1, f3 + fadd.s f6, f0, f2 + fadd.s f7, f1, f3 + + fst.s f6, t1, 0x00 + fst.s f7, t1, 0x04 + fst.s f4, t2, 0x00 + fst.s f5, t2, 0x04 + + addi.d a7, a7, 8 + +.BFL03: + andi t0, t6, 1 + beq $r0, t0, .BFL04 + + fldx.s f0, t4, a7 + fldx.s f2, t5, a7 + + fsub.s f4, f0, f2 + fadd.s f6, f0, f2 + + fstx.s f6, t4, a7 + fstx.s f4, t5, a7 + + addi.d a7, a7, 4 + +.BFL04: +endfunc + + +/* void vector_fmul_scalar_lsx(float *dst, const float *src, float mul, + int len) */ +function vector_fmul_scalar_lsx + move a6, $r0 + move a7, $r0 + + move t4, a0 + move t5, a1 + move t6, a2 + + vpermi.w vr0, vr0, 0x00 + + srai.d t0, t6, 2 + beq a6, t0, .BFS02 + +.BFS01: + vldx vr1, t5, a7 + + vfmul.s vr2, vr1, vr0 + + vstx vr2, t4, a7 + + addi.d a7, a7, 16 + addi.d a6, a6, 1 + blt a6, t0, .BFS01 + +.BFS02: + andi t0, t6, 2 + beq $r0, t0, .BFS03 + + add.d t1, t5, a7 + add.d t2, t4, a7 + + fld.s f1, t1, 0x00 + fld.s f2, t1, 0x04 + + fmul.s f3, f1, f0 + fmul.s f4, f2, f0 + + fst.s f3, t2, 0x00 + fst.s f4, t2, 0x04 + + addi.d a7, a7, 8 + +.BFS03: + andi t0, t6, 1 + beq $r0, t0, .BFS04 + + fldx.s f1, t5, a7 + + fmul.s f3, f1, f0 + + fstx.s f3, t4, a7 + + addi.d a7, a7, 4 + +.BFS04: +endfunc \ No newline at end of file diff --git a/libavutil/loongarch/float_dsp.h b/libavutil/loongarch/float_dsp.h new file mode 100644 index 0000000000..644c1f3713 --- /dev/null +++ b/libavutil/loongarch/float_dsp.h @@ -0,0 +1,32 @@ +/* + * Copyright (c) 2024 Loongson Technology Corporation Limited + * This file is part of FFmpeg. + * + * FFmpeg is free software; you can redistribute it and/or + * modify it under the terms of the GNU Lesser General Public + * License as published by the Free Software Foundation; either + * version 2.1 of the License, or (at your option) any later version. + * + * FFmpeg is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + * Lesser General Public License for more details. + * + * You should have received a copy of the GNU Lesser General Public + * License along with FFmpeg; if not, write to the Free Software + * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA + */ + +#ifndef AVUTIL_LOONGARCH_FLOAT_DSP_H +#define AVUTIL_LOONGARCH_FLOAT_DSP_H + +#include "libavutil/float_dsp.h" + +void vector_fmul_window_lsx(float *dst, const float *src0, + const float *src1, const float *win, int len); + +void butterflies_float_lsx(float *restrict v1, float *restrict v2, int len); + +void vector_fmul_scalar_lsx(float *dst, const float *src, float mul, int len); + +#endif /* AVUTIL_LOONGARCH_FLOAT_DSP_H */ \ No newline at end of file diff --git a/libavutil/loongarch/float_dsp_init_loongarch.c b/libavutil/loongarch/float_dsp_init_loongarch.c new file mode 100644 index 0000000000..592ba78058 --- /dev/null +++ b/libavutil/loongarch/float_dsp_init_loongarch.c @@ -0,0 +1,35 @@ +/* + * This file is part of FFmpeg. + * + * FFmpeg is free software; you can redistribute it and/or + * modify it under the terms of the GNU Lesser General Public + * License as published by the Free Software Foundation; either + * version 2.1 of the License, or (at your option) any later version. + * + * FFmpeg is distributed in the hope that it will be useful, + * but WITHOUT ANY WARRANTY; without even the implied warranty of + * MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the GNU + * Lesser General Public License for more details. + * + * You should have received a copy of the GNU Lesser General Public + * License along with FFmpeg; if not, write to the Free Software + * Foundation, Inc., 51 Franklin Street, Fifth Floor, Boston, MA 02110-1301 USA + */ + +#include "config.h" + +#include "libavutil/attributes.h" +#include "libavutil/cpu.h" +#include "float_dsp.h" +#include "libavutil/loongarch/cpu.h" + +av_cold void ff_float_dsp_init_loongarch(AVFloatDSPContext *fdsp) +{ + int cpu_flags = av_get_cpu_flags(); + + if (have_lsx(cpu_flags)) { + fdsp->vector_fmul_window = vector_fmul_window_lsx; + fdsp->butterflies_float = butterflies_float_lsx; + fdsp->vector_fmul_scalar = vector_fmul_scalar_lsx; + } +}