Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

perf: reduce the read and write of shared memory in the FusedAddRMSNormKernel #592

Merged
merged 4 commits into from
Nov 9, 2024

Conversation

Abatom
Copy link
Contributor

@Abatom Abatom commented Nov 8, 2024

Use vec_t<float, VEC_SIZE> x_vec to reduce the number of read and write operations to shared memory.

@@ -133,6 +133,8 @@ __global__ void FusedAddRMSNormKernel(T* __restrict__ input, T* __restrict__ res
input_vec.fill(0.f);
vec_t<T, VEC_SIZE> residual_vec;
residual_vec.fill(0.f);
vec_t<float, VEC_SIZE> x_vec;
Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I wrote this kernel in August #419, and you can actually use https://pytorch.org/docs/stable/benchmark_utils.html to add a benchmark. This way, you can know whether there is a performance improvement before and after the changes.

Copy link
Member

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

I wrote this kernel in August #419, and you can actually use https://pytorch.org/docs/stable/benchmark_utils.html to add a benchmark. This way, you can know whether there is a performance improvement before and after the changes.

Okay, I'll look into this, but I've analyzed this PR using Nsign Compute and found that the performance is about the same as the code before the precision improvement(#587).

Copy link
Collaborator

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Copy link
Contributor Author

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

Okay,I'll try to write a benchmark test like this.

@Abatom
Copy link
Contributor Author

Abatom commented Nov 8, 2024

On A800

Original

#python3 benchmarks/bench_fused_add_rmsnorm.py
batch_size:   1, hidden_size:   111, dtype: float16 , latency:  9us, throughput:   0.127GB/s
batch_size:   1, hidden_size:   500, dtype: float16 , latency:  9us, throughput:   0.556GB/s
batch_size:   1, hidden_size:  1024, dtype: float16 , latency:  8us, throughput:   1.212GB/s
batch_size:   1, hidden_size:  3072, dtype: float16 , latency:  9us, throughput:   3.593GB/s
batch_size:   1, hidden_size:  4096, dtype: float16 , latency:  8us, throughput:   5.004GB/s
batch_size:   1, hidden_size:  8192, dtype: float16 , latency:  9us, throughput:   8.731GB/s
---
batch_size:  19, hidden_size:   111, dtype: float16 , latency:  8us, throughput:   2.138GB/s
batch_size:  19, hidden_size:   500, dtype: float16 , latency:  8us, throughput:  10.063GB/s
batch_size:  19, hidden_size:  1024, dtype: float16 , latency:  8us, throughput:  20.380GB/s
batch_size:  19, hidden_size:  3072, dtype: float16 , latency:  9us, throughput:  54.253GB/s
batch_size:  19, hidden_size:  4096, dtype: float16 , latency:  9us, throughput:  71.412GB/s
batch_size:  19, hidden_size:  8192, dtype: float16 , latency:  9us, throughput: 136.994GB/s
---
batch_size:  99, hidden_size:   111, dtype: float16 , latency:  8us, throughput:  11.466GB/s
batch_size:  99, hidden_size:   500, dtype: float16 , latency:  9us, throughput:  45.655GB/s
batch_size:  99, hidden_size:  1024, dtype: float16 , latency:  8us, throughput:  98.419GB/s
batch_size:  99, hidden_size:  3072, dtype: float16 , latency:  9us, throughput: 261.573GB/s
batch_size:  99, hidden_size:  4096, dtype: float16 , latency: 10us, throughput: 337.045GB/s
batch_size:  99, hidden_size:  8192, dtype: float16 , latency: 12us, throughput: 533.967GB/s
---
batch_size: 989, hidden_size:   111, dtype: float16 , latency:  9us, throughput:  97.215GB/s
batch_size: 989, hidden_size:   500, dtype: float16 , latency: 10us, throughput: 392.762GB/s
batch_size: 989, hidden_size:  1024, dtype: float16 , latency: 12us, throughput: 670.266GB/s
batch_size: 989, hidden_size:  3072, dtype: float16 , latency: 24us, throughput: 1025.302GB/s
batch_size: 989, hidden_size:  4096, dtype: float16 , latency: 28us, throughput: 1155.931GB/s
batch_size: 989, hidden_size:  8192, dtype: float16 , latency: 48us, throughput: 1351.178GB/s
---

v1

#python3 benchmarks/bench_fused_add_rmsnorm.py
batch_size:   1, hidden_size:   111, dtype: float16 , latency:  9us, throughput:   0.124GB/s
batch_size:   1, hidden_size:   500, dtype: float16 , latency:  8us, throughput:   0.592GB/s
batch_size:   1, hidden_size:  1024, dtype: float16 , latency:  7us, throughput:   1.369GB/s
batch_size:   1, hidden_size:  3072, dtype: float16 , latency:  8us, throughput:   3.775GB/s
batch_size:   1, hidden_size:  4096, dtype: float16 , latency:  8us, throughput:   4.885GB/s
batch_size:   1, hidden_size:  8192, dtype: float16 , latency:  9us, throughput:   8.759GB/s
---
batch_size:  19, hidden_size:   111, dtype: float16 , latency:  8us, throughput:   2.244GB/s
batch_size:  19, hidden_size:   500, dtype: float16 , latency:  8us, throughput:   9.345GB/s
batch_size:  19, hidden_size:  1024, dtype: float16 , latency:  8us, throughput:  20.047GB/s
batch_size:  19, hidden_size:  3072, dtype: float16 , latency:  8us, throughput:  56.730GB/s
batch_size:  19, hidden_size:  4096, dtype: float16 , latency:  9us, throughput:  73.456GB/s
batch_size:  19, hidden_size:  8192, dtype: float16 , latency: 10us, throughput: 129.034GB/s
---
batch_size:  99, hidden_size:   111, dtype: float16 , latency:  8us, throughput:  11.182GB/s
batch_size:  99, hidden_size:   500, dtype: float16 , latency:  9us, throughput:  45.244GB/s
batch_size:  99, hidden_size:  1024, dtype: float16 , latency:  8us, throughput:  97.719GB/s
batch_size:  99, hidden_size:  3072, dtype: float16 , latency: 10us, throughput: 255.074GB/s
batch_size:  99, hidden_size:  4096, dtype: float16 , latency: 10us, throughput: 309.801GB/s
batch_size:  99, hidden_size:  8192, dtype: float16 , latency: 13us, throughput: 498.308GB/s
---
batch_size: 989, hidden_size:   111, dtype: float16 , latency:  9us, throughput: 100.372GB/s
batch_size: 989, hidden_size:   500, dtype: float16 , latency: 11us, throughput: 346.132GB/s
batch_size: 989, hidden_size:  1024, dtype: float16 , latency: 14us, throughput: 586.412GB/s
batch_size: 989, hidden_size:  3072, dtype: float16 , latency: 25us, throughput: 971.120GB/s
batch_size: 989, hidden_size:  4096, dtype: float16 , latency: 30us, throughput: 1094.190GB/s
batch_size: 989, hidden_size:  8192, dtype: float16 , latency: 50us, throughput: 1295.277GB/s
---

v2: Improve the precision of the FusedAddRMSNormKernel function #587

#python3 benchmarks/bench_fused_add_rmsnorm.py
batch_size:   1, hidden_size:   111, dtype: float16 , latency:  9us, throughput:   0.127GB/s
batch_size:   1, hidden_size:   500, dtype: float16 , latency:  9us, throughput:   0.560GB/s
batch_size:   1, hidden_size:  1024, dtype: float16 , latency:  9us, throughput:   1.121GB/s
batch_size:   1, hidden_size:  3072, dtype: float16 , latency:  9us, throughput:   3.494GB/s
batch_size:   1, hidden_size:  4096, dtype: float16 , latency:  9us, throughput:   4.588GB/s
batch_size:   1, hidden_size:  8192, dtype: float16 , latency: 11us, throughput:   7.656GB/s
---
batch_size:  19, hidden_size:   111, dtype: float16 , latency:  8us, throughput:   2.120GB/s
batch_size:  19, hidden_size:   500, dtype: float16 , latency:  8us, throughput:   9.461GB/s
batch_size:  19, hidden_size:  1024, dtype: float16 , latency:  8us, throughput:  18.881GB/s
batch_size:  19, hidden_size:  3072, dtype: float16 , latency:  9us, throughput:  54.199GB/s
batch_size:  19, hidden_size:  4096, dtype: float16 , latency: 10us, throughput:  64.051GB/s
batch_size:  19, hidden_size:  8192, dtype: float16 , latency: 12us, throughput: 106.248GB/s
---
batch_size:  99, hidden_size:   111, dtype: float16 , latency:  8us, throughput:  11.392GB/s
batch_size:  99, hidden_size:   500, dtype: float16 , latency:  8us, throughput:  49.162GB/s
batch_size:  99, hidden_size:  1024, dtype: float16 , latency:  9us, throughput:  92.030GB/s
batch_size:  99, hidden_size:  3072, dtype: float16 , latency: 10us, throughput: 235.945GB/s
batch_size:  99, hidden_size:  4096, dtype: float16 , latency: 11us, throughput: 303.909GB/s
batch_size:  99, hidden_size:  8192, dtype: float16 , latency: 13us, throughput: 483.708GB/s
---
batch_size: 989, hidden_size:   111, dtype: float16 , latency:  8us, throughput: 103.845GB/s
batch_size: 989, hidden_size:   500, dtype: float16 , latency: 10us, throughput: 382.028GB/s
batch_size: 989, hidden_size:  1024, dtype: float16 , latency: 14us, throughput: 563.498GB/s
batch_size: 989, hidden_size:  3072, dtype: float16 , latency: 26us, throughput: 927.479GB/s
batch_size: 989, hidden_size:  4096, dtype: float16 , latency: 32us, throughput: 1009.666GB/s
batch_size: 989, hidden_size:  8192, dtype: float16 , latency: 54us, throughput: 1207.263GB/s
---

v3: This PR

#python3 benchmarks/bench_fused_add_rmsnorm.py
batch_size:   1, hidden_size:   111, dtype: float16 , latency:  9us, throughput:   0.126GB/s
batch_size:   1, hidden_size:   500, dtype: float16 , latency:  9us, throughput:   0.584GB/s
batch_size:   1, hidden_size:  1024, dtype: float16 , latency:  9us, throughput:   1.123GB/s
batch_size:   1, hidden_size:  3072, dtype: float16 , latency:  9us, throughput:   3.578GB/s
batch_size:   1, hidden_size:  4096, dtype: float16 , latency:  9us, throughput:   4.644GB/s
batch_size:   1, hidden_size:  8192, dtype: float16 , latency: 10us, throughput:   8.458GB/s
---
batch_size:  19, hidden_size:   111, dtype: float16 , latency:  8us, throughput:   2.136GB/s
batch_size:  19, hidden_size:   500, dtype: float16 , latency:  8us, throughput:   9.398GB/s
batch_size:  19, hidden_size:  1024, dtype: float16 , latency:  8us, throughput:  19.087GB/s
batch_size:  19, hidden_size:  3072, dtype: float16 , latency:  9us, throughput:  54.074GB/s
batch_size:  19, hidden_size:  4096, dtype: float16 , latency:  9us, throughput:  70.508GB/s
batch_size:  19, hidden_size:  8192, dtype: float16 , latency: 10us, throughput: 124.474GB/s
---
batch_size:  99, hidden_size:   111, dtype: float16 , latency:  8us, throughput:  10.703GB/s
batch_size:  99, hidden_size:   500, dtype: float16 , latency:  8us, throughput:  48.536GB/s
batch_size:  99, hidden_size:  1024, dtype: float16 , latency:  9us, throughput:  93.118GB/s
batch_size:  99, hidden_size:  3072, dtype: float16 , latency: 10us, throughput: 247.778GB/s
batch_size:  99, hidden_size:  4096, dtype: float16 , latency: 10us, throughput: 316.723GB/s
batch_size:  99, hidden_size:  8192, dtype: float16 , latency: 12us, throughput: 543.516GB/s
---
batch_size: 989, hidden_size:   111, dtype: float16 , latency:  8us, throughput: 104.044GB/s
batch_size: 989, hidden_size:   500, dtype: float16 , latency: 11us, throughput: 369.634GB/s
batch_size: 989, hidden_size:  1024, dtype: float16 , latency: 12us, throughput: 658.708GB/s
batch_size: 989, hidden_size:  3072, dtype: float16 , latency: 24us, throughput: 1003.270GB/s
batch_size: 989, hidden_size:  4096, dtype: float16 , latency: 29us, throughput: 1127.658GB/s
batch_size: 989, hidden_size:  8192, dtype: float16 , latency: 49us, throughput: 1317.664GB/s
---

This PR has essentially the same performance as the original, slightly lower than the original, but for the sake of precision, it can be improved by an order of magnitude (from 1e-2 to 1e-3), which I believe is acceptable.
cc @yzh119 @zhyncs

@Abatom Abatom requested review from yzh119 and zhyncs November 8, 2024 08:42
@jeejeelee
Copy link
Contributor

This PR has essentially the same performance as the original, slightly lower than the original, but for the sake of precision, it can be improved by an order of magnitude (from 1e-2 to 1e-3), which I believe is acceptable.

I tested #587 , and if dtype=torch.bfloat16, 1e-3 is too strict

@zhyncs
Copy link
Member

zhyncs commented Nov 8, 2024

slightly lower than the original

@Abatom Could you clarify the main purpose of this PR, since it doesn't appear to focus on performance improvements? I'm also curious about any other potential benefits, as the accuracy improvements don't seem significant.

@zhyncs zhyncs changed the title perf: reduce the read and write of shared memory in the FusedAddRMSNormKernel minor: reduce the read and write of shared memory in the FusedAddRMSNormKernel Nov 8, 2024
@Abatom
Copy link
Contributor Author

Abatom commented Nov 8, 2024

The main purpose of this PR is to fix the performance degradation issue caused by #587. Now, compared to the version before #587, the performance is almost unchanged, but the precision has been improved (from 1e-2 to 1e-3), and benchmark tests for fused_add_rmsnorm have been added.

@yzh119
Copy link
Collaborator

yzh119 commented Nov 8, 2024

@zhyncs this PR uses vectorized load/store to shared memory to accelerate #587 .

Copy link
Collaborator

@yzh119 yzh119 left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

The achieved bandwidth might have been underestimated because we didn't count memory access to residual.

benchmarks/bench_fused_add_rmsnorm.py Outdated Show resolved Hide resolved
@Abatom Abatom requested a review from yzh119 November 9, 2024 02:23
Copy link
Collaborator

@yzh119 yzh119 left a comment

Choose a reason for hiding this comment

The reason will be displayed to describe this comment to others. Learn more.

LGTM, updated results on H100:

---
batch_size:   1, hidden_size:   111, dtype: float16 , latency:  6us, throughput:   0.180GB/s
batch_size:   1, hidden_size:   500, dtype: float16 , latency:  6us, throughput:   0.850GB/s
batch_size:   1, hidden_size:  1024, dtype: float16 , latency:  6us, throughput:   1.706GB/s
batch_size:   1, hidden_size:  3072, dtype: float16 , latency:  6us, throughput:   4.962GB/s
batch_size:   1, hidden_size:  4096, dtype: float16 , latency:  6us, throughput:   6.540GB/s
batch_size:   1, hidden_size:  8192, dtype: float16 , latency:  7us, throughput:  11.960GB/s
---
batch_size:  19, hidden_size:   111, dtype: float16 , latency:  6us, throughput:   2.968GB/s
batch_size:  19, hidden_size:   500, dtype: float16 , latency:  6us, throughput:  13.038GB/s
batch_size:  19, hidden_size:  1024, dtype: float16 , latency:  6us, throughput:  26.839GB/s
batch_size:  19, hidden_size:  3072, dtype: float16 , latency:  6us, throughput:  76.400GB/s
batch_size:  19, hidden_size:  4096, dtype: float16 , latency:  6us, throughput:  98.126GB/s
batch_size:  19, hidden_size:  8192, dtype: float16 , latency:  7us, throughput: 182.439GB/s
---
batch_size:  99, hidden_size:   111, dtype: float16 , latency:  6us, throughput:  15.074GB/s
batch_size:  99, hidden_size:   500, dtype: float16 , latency:  6us, throughput:  65.834GB/s
batch_size:  99, hidden_size:  1024, dtype: float16 , latency:  6us, throughput: 131.731GB/s
batch_size:  99, hidden_size:  3072, dtype: float16 , latency:  7us, throughput: 364.698GB/s
batch_size:  99, hidden_size:  4096, dtype: float16 , latency:  7us, throughput: 463.593GB/s
batch_size:  99, hidden_size:  8192, dtype: float16 , latency:  8us, throughput: 769.438GB/s
---
batch_size: 989, hidden_size:   111, dtype: float16 , latency:  6us, throughput: 137.252GB/s
batch_size: 989, hidden_size:   500, dtype: float16 , latency:  7us, throughput: 547.651GB/s
batch_size: 989, hidden_size:  1024, dtype: float16 , latency:  8us, throughput: 964.887GB/s
batch_size: 989, hidden_size:  3072, dtype: float16 , latency: 14us, throughput: 1701.064GB/s
batch_size: 989, hidden_size:  4096, dtype: float16 , latency: 17us, throughput: 1914.022GB/s
batch_size: 989, hidden_size:  8192, dtype: float16 , latency: 29us, throughput: 2242.014GB/s
---

@yzh119 yzh119 merged commit 2043ca2 into flashinfer-ai:main Nov 9, 2024
@zhyncs
Copy link
Member

zhyncs commented Nov 9, 2024

@zhyncs this PR uses vectorized load/store to shared memory to accelerate #587 .

cool!

@Abatom Abatom changed the title minor: reduce the read and write of shared memory in the FusedAddRMSNormKernel perf: reduce the read and write of shared memory in the FusedAddRMSNormKernel Nov 11, 2024
yzh119 added a commit that referenced this pull request Nov 24, 2024
gemma-style rmsnorm kernels (introduced in #477 ) are similar to
original rmsnorm kernel, and we should use the same kernel for them.
This PR cleans up duplicate code and unifies the kernels for gemma-style
and original rmsnorm kernels.

The precision improvements
(#587,
#592) are kept in this
PR.
yzh119 added a commit that referenced this pull request Nov 25, 2024
…apes (#636)

This PR fixes the issue #634, which is brought by #592 .
If we want to use 16-bytes vectorized read/write, we need to confirm the
address is aligned to 16 bytes.
When `num_warps` is not a multiple of 4 (4*sizeof(float) = 16), the
address of `smem + num_warps` might not align to 16 bytes.

We can fix this by shifting the start offset of vectorized read/write to
`smem + ceil_div(num_warps, 4) * 4` to force the alignment.

cc @ovowei @Abatom
yzh119 added a commit that referenced this pull request Dec 17, 2024
🤖 I have created a release *beep* *boop*
---


##
[0.2.0](v0.1.6...v0.2.0)
(2024-12-17)

[Release
Blog](https://flashinfer.ai/2024/12/16/flashinfer-v02-release.html).

### Features

* add `rotary_dim` argument to rope APIs for partial apply rope
([#599](#599))
([eb9bc71](eb9bc71))
* add a `use_softmax` field in variant class
([#533](#533))
([d81af97](d81af97))
* add an option `non_blocking` to plan function
([#622](#622))
([560af6f](560af6f))
* add gemma_rmsnorm and gemma_fused_add_rmsnorm
([#477](#477))
([1a6b17e](1a6b17e))
* add group size 3 to GQA decode dispatch
([#558](#558))
([6227562](6227562))
* add JIT compilation support for FA3 templates
([#672](#672))
([d4e8d79](d4e8d79))
* allow the cascade kernels to be executed using varying sequence
lenghts ([#627](#627))
([92ac440](92ac440))
* CUDAGraph compatibility of multi-level cascade inference APIs
([#586](#586))
([2332e8a](2332e8a))
* fix the maximal grid dimension in prefill planning with CUDA graphs
([#639](#639))
([86ca89a](86ca89a))
* improve the precision of the FusedAddRMSNormKernel function
([#587](#587))
([c7dc921](c7dc921))
* JIT compilation
([#507](#507))
([3613a5b](3613a5b))
* modify group-gemm stage number
([#497](#497))
([52dab1d](52dab1d))
* non-contiguous query with paged kv cache
([#553](#553))
([89f2c4a](89f2c4a))
* pass a dynamic token count to the cascade kernels
([#635](#635))
([5fe9f7d](5fe9f7d))
* simplify prefill JIT compilation
([#605](#605))
([fe4f898](fe4f898))
* specify gemm backend
([#648](#648))
([0cc1a51](0cc1a51))
* support cached cos/sin in rope APIs
([#585](#585))
([83e541d](83e541d))
* support huggingface transformer style rope interface
([#568](#568))
([4f40420](4f40420))
* support sm90 cutlass group gemm
([#509](#509))
([794bdda](794bdda))
* torch custom_op fix for rope
([#569](#569))
([3e104bc](3e104bc))
* torch custom_op support: norm
([#552](#552))
([f6e0010](f6e0010))
* torch.compile and custom_op support
([#554](#554))
([9bf916f](9bf916f))
* warmup for jit kernel tests
([#629](#629))
([8f5f349](8f5f349))


### Bug Fixes

* AOT compiler flags on non-sm90
([#522](#522))
([0aa4726](0aa4726))
* batch decode kernel redundant store output to gmem
([#505](#505))
([90e42a7](90e42a7))
* compatible with torch 2.2
([#478](#478))
([ac41d1b](ac41d1b))
* #452
([b53a46f](b53a46f))
* remove redundant load
([#495](#495))
([2de16b0](2de16b0))
* update bmm fp8 test
([#487](#487))
([45eac04](45eac04))


### Performance Improvements

* accelerate JIT compilation speed
([#618](#618))
([eaf73fd](eaf73fd))
* Dense and sparse customizable flashattention-3 template
([#667](#667))
([51236c9](51236c9))
* fix prefill kernel performance degradation (step 1)
([#602](#602))
([595cf60](595cf60))
* fix the performance issue of `append_paged_kv_cache`
([#588](#588))
([e15f7c9](e15f7c9))
* improve parallelism in RoPE with pos_ids
([#609](#609))
([ff05155](ff05155))
* improve plan performance by using non-blocking memcpy
([#547](#547))
([41ebe6d](41ebe6d))
* reduce the read and write of shared memory in the
FusedAddRMSNormKernel
([#592](#592))
([2043ca2](2043ca2))
* reduce total_num_tiles_q by one
([#644](#644))
([553ace5](553ace5))
* remove unnecessary contiguous operation in block sparse attention
([#561](#561))
([7a7ad46](7a7ad46))
* speedup jit compilation of prefill attention kernels
([#632](#632))
([a059586](a059586))
* use cuda-core implemention for io-bound block-sparse attention
([#560](#560))
([3fbf028](3fbf028))

---
This PR was generated with [Release
Please](https://github.com/googleapis/release-please). See
[documentation](https://github.com/googleapis/release-please#release-please).

---------

Co-authored-by: github-actions[bot] <41898282+github-actions[bot]@users.noreply.github.com>
Co-authored-by: Zihao Ye <expye@outlook.com>
Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

Successfully merging this pull request may close these issues.

4 participants