ArthurZucker/best_benchmark.py

## best_benchmark.py
FRANCE_ARTICLE = """<s>Marseille, France (CNN)The French prosecutor leading an investigation into the crash of Germanwings Flight 9525 insisted Wednesday that he was not aware of any video footage from on board the plane. Marseille prosecutor Brice Robin told CNN that "so far no videos were used in the crash investigation." He added, "A person who has such a video needs to immediately give it to the investigators." Robin\'s comments follow claims by two magazines, German daily Bild and French Paris Match, of a cell phone video showing the harrowing final seconds from on board Germanwings Flight 9525 as it crashed into the French Alps. All 150 on board were killed. Paris Match and Bild reported that the video was recovered from a phone at the wreckage site. The two publications described the supposed video, but did not post it on their websites. The publications said that they watched the video, which was found by a source close to the investigation. \"One can hear cries of 'My God' in several languages,\" Paris Match reported. "Metallic banging can also be heard more than three times, perhaps of the pilot trying to open the cockpit door with a heavy object.  Towards the end, after a heavy shake, stronger than the others, the screaming intensifies. Then nothing." "It is a very disturbing scene," said Julian Reichelt, editor-in-chief of Bild online. An official with France's accident investigation agency, the BEA, said the agency is not aware of any such video. Lt. Col. Jean-Marc Menichini, a French Gendarmerie spokesman in charge of communications on rescue efforts around the Germanwings crash site, told CNN that the reports were "completely wrong" and "unwarranted." Cell phones have been collected at the site, he said, but that they "hadn\'t been exploited yet." Menichini said he believed the cell phones would need to be sent to the Criminal Research Institute in Rosny sous-Bois, near Paris, in order to be analyzed by specialized technicians working hand-in-hand with investigators. But none of the cell phones found so far have been sent to the institute, Menichini said. Asked whether staff involved in the search could have leaked a memory card to the media, Menichini answered with a categorical "no." Reichelt told "Erin Burnett: Outfront" that he had watched the video and stood by the report, saying Bild and Paris Match are "very confident" that the clip is real. He noted that investigators only revealed they\'d recovered cell phones from the crash site after Bild and Paris Match published their reports. "That is something we did not know before. ... Overall we can say many things of the investigation weren't revealed by the investigation at the beginning," he said. What was mental state of Germanwings co-pilot? German airline Lufthansa confirmed Tuesday that co-pilot Andreas Lubitz had battled depression years before he took the controls of Germanwings Flight 9525, which he's accused of deliberately crashing last week in the French Alps. Lubitz told his Lufthansa flight training school in 2009 that he had a "previous episode of severe depression," the airline said Tuesday. Email correspondence between Lubitz and the school discovered in an internal investigation, Lufthansa said, included medical documents he submitted in connection with resuming his flight training. The announcement indicates that Lufthansa, the parent company of Germanwings, knew of Lubitz's battle with depression, allowed him to continue training and ultimately put him in the cockpit. Lufthansa, whose CEO Carsten Spohr previously said Lubitz was 100% fit to fly, described its statement Tuesday as a "swift and seamless clarification" and said it was sharing the information and documents -- including training and medical records -- with public prosecutors. Spohr traveled to the crash site Wednesday, where recovery teams have been working for the past week to recover human remains and plane debris scattered across a steep mountainside. He saw the crisis center set up in Seyne-les-Alpes, laid a wreath in the village of Le Vernet, closer to the crash site, where grieving families have left flowers at a simple stone memorial. Menichini told CNN late Tuesday that no visible human remains were left at the site but recovery teams would keep searching. French President Francois Hollande, speaking Tuesday, said that it should be possible to identify all the victims using DNA analysis by the end of the week, sooner than authorities had previously suggested. In the meantime, the recovery of the victims' personal belongings will start Wednesday, Menichini said. Among those personal belongings could be more cell phones belonging to the 144 passengers and six crew on board. Check out the latest from our correspondents . The details about Lubitz's correspondence with the flight school during his training were among several developments as investigators continued to delve into what caused the crash and Lubitz's possible motive for downing the jet. A Lufthansa spokesperson told CNN on Tuesday that Lubitz had a valid medical certificate, had passed all his examinations and "held all the licenses required." Earlier, a spokesman for the prosecutor\'s office in Dusseldorf, Christoph Kumpa, said medical records reveal Lubitz suffered from suicidal tendencies at some point before his aviation career and underwent psychotherapy before he got his pilot's license. Kumpa emphasized there's no evidence suggesting Lubitz was suicidal or acting aggressively before the crash. Investigators are looking into whether Lubitz feared his medical condition would cause him to lose his pilot's license, a European government official briefed on the investigation told CNN on Tuesday. While flying was "a big part of his life," the source said, it\'s only one theory being considered. Another source, a law enforcement official briefed on the investigation, also told CNN that authorities believe the primary motive for Lubitz to bring down the plane was that he feared he would not be allowed to fly because of his medical problems. Lubitz's girlfriend told investigators he had seen an eye doctor and a neuropsychologist, both of whom deemed him unfit to work recently and concluded he had psychological issues, the European government official said. But no matter what details emerge about his previous mental health struggles, there's more to the story, said Brian Russell, a forensic psychologist. "Psychology can explain why somebody would turn rage inward on themselves about the fact that maybe they weren't going to keep doing their job and they're upset about that and so they're suicidal," he said. "But there is no mental illness that explains why somebody then feels entitled to also take that rage and turn it outward on 149 other people who had nothing to do with the person's problems." Germanwings crash compensation: What we know . Who was the captain of Germanwings Flight 9525? CNN's Margot Haddad reported from Marseille and Pamela Brown from Dusseldorf, while Laura Smith-Spark wrote from London. CNN's Frederik Pleitgen, Pamela Boykoff, Antonia Mortensen, Sandrine Amiel and Anna-Maja Rappard contributed to this report.""" * 40
import torch, transformers, time, tokenizers, os
os.environ["TOKENIZERS_PARALLELISM"] = "0"
import torch._dynamo.config
import torch._inductor.config
from transformers import AutoModelForCausalLM, AutoTokenizer, DynamicCache, set_seed  ,StaticCache
from typing import Optional
# use https://github.com/huggingface/transformers/pull/30476 branch of transformers
print(tokenizers.__version__)
print(transformers.__version__)

def sample(logits, temperature: float = 0.6, top_k: Optional[int] = 5):
    logits = logits / max(temperature, 1e-5)
    if top_k is not None:
        v, _ = torch.topk(logits, min(top_k, logits.size(-1)))
        pivot = v.select(-1, -1).unsqueeze(-1)
        logits = torch.where(logits < pivot, -float("Inf"), logits)
    probs_sort = torch.nn.functional.softmax(logits, dim=-1)[:, -1]
    probs = torch.empty_like(probs_sort).exponential_(1)
    idx_next = torch.argmax(probs / probs_sort, dim=-1, keepdim=True).to(dtype=torch.int)
    return idx_next

torch._inductor.config.coordinate_descent_tuning = True
torch._inductor.config.triton.unique_kernel_names = True
torch._inductor.config.fx_graph_cache = True
torch._dynamo.config.cache_size_limit = 32
# TORCH_LOGS="perf_hints,recompiles,graph_breaks" python bench_new.py &> logs.txt
# TORCH_LOGS="recompiles"

torch.set_float32_matmul_precision('high')
torch.set_printoptions(linewidth=200)  # you can better see how the mask is shaped

device = "cuda:3"
attn_implementation = "sdpa"
all_dtype = torch.bfloat16
model_id = "meta-llama/Llama-2-7b-hf"
PREFILL_NUM_ITER = 3
NUM_ITER = 100
max_cache_length = 2048


tokenizer = AutoTokenizer.from_pretrained(model_id, padding_side="left", pad_token="<|end_of_text|>")
model = AutoModelForCausalLM.from_pretrained(model_id, torch_dtype=all_dtype, attn_implementation=attn_implementation).to(device,all_dtype)
model = model.eval()
BIG_INPUT = tokenizer([FRANCE_ARTICLE]*8, return_tensors="pt").to(device).input_ids

for BS in [8, 4, 2, 1]:
    for seq_length in [1, 512, 1024, 1800]:
        input_ids = BIG_INPUT[:BS, :seq_length]

        ##########  eager dynamic cache ######################################################################
        generated_ids = torch.zeros(BS, NUM_ITER, dtype = torch.long, device=device)
        past_key_values = DynamicCache()
        set_seed(123)
        with torch.no_grad():
            cache_position = torch.arange(seq_length , device=device)
            outputs = model(input_ids, past_key_values=past_key_values,cache_position=cache_position, return_dict=False, use_cache = True)
            new_token = sample(outputs[0])
            generated_ids[:,0] = new_token[:,0]

            ### DRY RUN
            cache_position = torch.tensor([seq_length] , device=device)
            for i in range(1, PREFILL_NUM_ITER):
                outputs = model(new_token, past_key_values=past_key_values,cache_position=cache_position, return_dict=False, use_cache = True)
                new_token = sample(outputs[0])
                cache_position+=1
                generated_ids[:,i] = new_token[:,0]

            start = time.perf_counter()
            for i in range(PREFILL_NUM_ITER, NUM_ITER):
                outputs = model(new_token, past_key_values=past_key_values,cache_position=cache_position, return_dict=False, use_cache = True)
                new_token = sample(outputs[0])
                generated_ids[:,i] = new_token[:,0]
                cache_position+=1
            torch.cuda.synchronize()
            total_time = time.perf_counter() - start
            forward_average_time =  total_time / NUM_ITER
            print(f"{BS},\t{seq_length},\t{max_cache_length}: eager dynamic cache\t - Took : {total_time:.4f} sec total, {forward_average_time*1000:.3f} ms per forward, {(NUM_ITER/total_time):.2f} tokens/s")
            # print(f"{tokenizer.batch_decode(generated_ids.tolist())}")

        def decode_one_tokens(model, cur_token, input_pos, cache_position, past_key_values):
            logits = model(cur_token, position_ids=input_pos, past_key_values=past_key_values, cache_position=cache_position, return_dict=False, use_cache = True)[0]
            new_token = sample(logits)
            return new_token

        past_key_values = StaticCache(model.config, BS, max_cache_length, device=device, dtype=torch.bfloat16)
        with torch.no_grad():
            ### PREFILL ##################
            cache_position = torch.arange(seq_length , device=device)
            generated_ids = torch.zeros(BS, seq_length+NUM_ITER+PREFILL_NUM_ITER+1, dtype = torch.int, device=device)
            generated_ids[:, cache_position] = input_ids.to(device).to(torch.int)

            logits = model(input_ids, past_key_values=past_key_values, cache_position=cache_position, return_dict=False, use_cache = True)[0]
            next_token = sample(logits)
            generated_ids[:, seq_length] = next_token[:,0]

            decode_one_tokens = torch.compile(decode_one_tokens, mode="reduce-overhead", fullgraph=True)

            ### DRY RUN ##################
            cache_position = torch.tensor([seq_length] , device=device)
            for i in range(1, PREFILL_NUM_ITER):
                next_token = decode_one_tokens(model, next_token.clone(), None, cache_position, past_key_values)
                generated_ids.index_copy_(1, cache_position, next_token)
                cache_position+=1

            ### RUN ##################
            start = time.perf_counter()
            for i in range(PREFILL_NUM_ITER, NUM_ITER+PREFILL_NUM_ITER):
                next_token = decode_one_tokens(model, next_token.clone(), None, cache_position, past_key_values)
                generated_ids.index_copy_(1, cache_position, next_token)
                cache_position+=1

            torch.cuda.synchronize()
            total_time = time.perf_counter() - start
            forward_average_time =  total_time / NUM_ITER
            print(f"{BS},\t{seq_length},\t{max_cache_length}: compiled static cache\t - Took : {total_time:.4f} sec total, {forward_average_time*1000:.3f} ms per forward, {(NUM_ITER/total_time):.2f} tokens/s")
            # print(f"{tokenizer.batch_decode(generated_ids[:,seq_length:].tolist())}")
            #####################################################################################
	FRANCE_ARTICLE = """<s>Marseille, France (CNN)The French prosecutor leading an investigation into the crash of Germanwings Flight 9525 insisted Wednesday that he was not aware of any video footage from on board the plane. Marseille prosecutor Brice Robin told CNN that "so far no videos were used in the crash investigation." He added, "A person who has such a video needs to immediately give it to the investigators." Robin\'s comments follow claims by two magazines, German daily Bild and French Paris Match, of a cell phone video showing the harrowing final seconds from on board Germanwings Flight 9525 as it crashed into the French Alps. All 150 on board were killed. Paris Match and Bild reported that the video was recovered from a phone at the wreckage site. The two publications described the supposed video, but did not post it on their websites. The publications said that they watched the video, which was found by a source close to the investigation. \"One can hear cries of 'My God' in several languages,\" Paris Match reported. "Metallic banging can also be heard more than three times, perhaps of the pilot trying to open the cockpit door with a heavy object. Towards the end, after a heavy shake, stronger than the others, the screaming intensifies. Then nothing." "It is a very disturbing scene," said Julian Reichelt, editor-in-chief of Bild online. An official with France's accident investigation agency, the BEA, said the agency is not aware of any such video. Lt. Col. Jean-Marc Menichini, a French Gendarmerie spokesman in charge of communications on rescue efforts around the Germanwings crash site, told CNN that the reports were "completely wrong" and "unwarranted." Cell phones have been collected at the site, he said, but that they "hadn\'t been exploited yet." Menichini said he believed the cell phones would need to be sent to the Criminal Research Institute in Rosny sous-Bois, near Paris, in order to be analyzed by specialized technicians working hand-in-hand with investigators. But none of the cell phones found so far have been sent to the institute, Menichini said. Asked whether staff involved in the search could have leaked a memory card to the media, Menichini answered with a categorical "no." Reichelt told "Erin Burnett: Outfront" that he had watched the video and stood by the report, saying Bild and Paris Match are "very confident" that the clip is real. He noted that investigators only revealed they\'d recovered cell phones from the crash site after Bild and Paris Match published their reports. "That is something we did not know before. ... Overall we can say many things of the investigation weren't revealed by the investigation at the beginning," he said. What was mental state of Germanwings co-pilot? German airline Lufthansa confirmed Tuesday that co-pilot Andreas Lubitz had battled depression years before he took the controls of Germanwings Flight 9525, which he's accused of deliberately crashing last week in the French Alps. Lubitz told his Lufthansa flight training school in 2009 that he had a "previous episode of severe depression," the airline said Tuesday. Email correspondence between Lubitz and the school discovered in an internal investigation, Lufthansa said, included medical documents he submitted in connection with resuming his flight training. The announcement indicates that Lufthansa, the parent company of Germanwings, knew of Lubitz's battle with depression, allowed him to continue training and ultimately put him in the cockpit. Lufthansa, whose CEO Carsten Spohr previously said Lubitz was 100% fit to fly, described its statement Tuesday as a "swift and seamless clarification" and said it was sharing the information and documents -- including training and medical records -- with public prosecutors. Spohr traveled to the crash site Wednesday, where recovery teams have been working for the past week to recover human remains and plane debris scattered across a steep mountainside. He saw the crisis center set up in Seyne-les-Alpes, laid a wreath in the village of Le Vernet, closer to the crash site, where grieving families have left flowers at a simple stone memorial. Menichini told CNN late Tuesday that no visible human remains were left at the site but recovery teams would keep searching. French President Francois Hollande, speaking Tuesday, said that it should be possible to identify all the victims using DNA analysis by the end of the week, sooner than authorities had previously suggested. In the meantime, the recovery of the victims' personal belongings will start Wednesday, Menichini said. Among those personal belongings could be more cell phones belonging to the 144 passengers and six crew on board. Check out the latest from our correspondents . The details about Lubitz's correspondence with the flight school during his training were among several developments as investigators continued to delve into what caused the crash and Lubitz's possible motive for downing the jet. A Lufthansa spokesperson told CNN on Tuesday that Lubitz had a valid medical certificate, had passed all his examinations and "held all the licenses required." Earlier, a spokesman for the prosecutor\'s office in Dusseldorf, Christoph Kumpa, said medical records reveal Lubitz suffered from suicidal tendencies at some point before his aviation career and underwent psychotherapy before he got his pilot's license. Kumpa emphasized there's no evidence suggesting Lubitz was suicidal or acting aggressively before the crash. Investigators are looking into whether Lubitz feared his medical condition would cause him to lose his pilot's license, a European government official briefed on the investigation told CNN on Tuesday. While flying was "a big part of his life," the source said, it\'s only one theory being considered. Another source, a law enforcement official briefed on the investigation, also told CNN that authorities believe the primary motive for Lubitz to bring down the plane was that he feared he would not be allowed to fly because of his medical problems. Lubitz's girlfriend told investigators he had seen an eye doctor and a neuropsychologist, both of whom deemed him unfit to work recently and concluded he had psychological issues, the European government official said. But no matter what details emerge about his previous mental health struggles, there's more to the story, said Brian Russell, a forensic psychologist. "Psychology can explain why somebody would turn rage inward on themselves about the fact that maybe they weren't going to keep doing their job and they're upset about that and so they're suicidal," he said. "But there is no mental illness that explains why somebody then feels entitled to also take that rage and turn it outward on 149 other people who had nothing to do with the person's problems." Germanwings crash compensation: What we know . Who was the captain of Germanwings Flight 9525? CNN's Margot Haddad reported from Marseille and Pamela Brown from Dusseldorf, while Laura Smith-Spark wrote from London. CNN's Frederik Pleitgen, Pamela Boykoff, Antonia Mortensen, Sandrine Amiel and Anna-Maja Rappard contributed to this report.""" * 40
	import torch, transformers, time, tokenizers, os
	os.environ["TOKENIZERS_PARALLELISM"] = "0"
	import torch._dynamo.config
	import torch._inductor.config
	from transformers import AutoModelForCausalLM, AutoTokenizer, DynamicCache, set_seed ,StaticCache
	from typing import Optional
	# use https://github.com/huggingface/transformers/pull/30476 branch of transformers
	print(tokenizers.__version__)
	print(transformers.__version__)

	def sample(logits, temperature: float = 0.6, top_k: Optional[int] = 5):
	logits = logits / max(temperature, 1e-5)
	if top_k is not None:
	v, _ = torch.topk(logits, min(top_k, logits.size(-1)))
	pivot = v.select(-1, -1).unsqueeze(-1)
	logits = torch.where(logits < pivot, -float("Inf"), logits)
	probs_sort = torch.nn.functional.softmax(logits, dim=-1)[:, -1]
	probs = torch.empty_like(probs_sort).exponential_(1)
	idx_next = torch.argmax(probs / probs_sort, dim=-1, keepdim=True).to(dtype=torch.int)
	return idx_next

	torch._inductor.config.coordinate_descent_tuning = True
	torch._inductor.config.triton.unique_kernel_names = True
	torch._inductor.config.fx_graph_cache = True
	torch._dynamo.config.cache_size_limit = 32
	# TORCH_LOGS="perf_hints,recompiles,graph_breaks" python bench_new.py &> logs.txt
	# TORCH_LOGS="recompiles"

	torch.set_float32_matmul_precision('high')
	torch.set_printoptions(linewidth=200) # you can better see how the mask is shaped

	device = "cuda:3"
	attn_implementation = "sdpa"
	all_dtype = torch.bfloat16
	model_id = "meta-llama/Llama-2-7b-hf"
	PREFILL_NUM_ITER = 3
	NUM_ITER = 100
	max_cache_length = 2048


	tokenizer = AutoTokenizer.from_pretrained(model_id, padding_side="left", pad_token="<\|end_of_text\|>")
	model = AutoModelForCausalLM.from_pretrained(model_id, torch_dtype=all_dtype, attn_implementation=attn_implementation).to(device,all_dtype)
	model = model.eval()
	BIG_INPUT = tokenizer([FRANCE_ARTICLE]*8, return_tensors="pt").to(device).input_ids

	for BS in [8, 4, 2, 1]:
	for seq_length in [1, 512, 1024, 1800]:
	input_ids = BIG_INPUT[:BS, :seq_length]

	########## eager dynamic cache ######################################################################
	generated_ids = torch.zeros(BS, NUM_ITER, dtype = torch.long, device=device)
	past_key_values = DynamicCache()
	set_seed(123)
	with torch.no_grad():
	cache_position = torch.arange(seq_length , device=device)
	outputs = model(input_ids, past_key_values=past_key_values,cache_position=cache_position, return_dict=False, use_cache = True)
	new_token = sample(outputs[0])
	generated_ids[:,0] = new_token[:,0]

	### DRY RUN
	cache_position = torch.tensor([seq_length] , device=device)
	for i in range(1, PREFILL_NUM_ITER):
	outputs = model(new_token, past_key_values=past_key_values,cache_position=cache_position, return_dict=False, use_cache = True)
	new_token = sample(outputs[0])
	cache_position+=1
	generated_ids[:,i] = new_token[:,0]

	start = time.perf_counter()
	for i in range(PREFILL_NUM_ITER, NUM_ITER):
	outputs = model(new_token, past_key_values=past_key_values,cache_position=cache_position, return_dict=False, use_cache = True)
	new_token = sample(outputs[0])
	generated_ids[:,i] = new_token[:,0]
	cache_position+=1
	torch.cuda.synchronize()
	total_time = time.perf_counter() - start
	forward_average_time = total_time / NUM_ITER
	print(f"{BS},\t{seq_length},\t{max_cache_length}: eager dynamic cache\t - Took : {total_time:.4f} sec total, {forward_average_time*1000:.3f} ms per forward, {(NUM_ITER/total_time):.2f} tokens/s")
	# print(f"{tokenizer.batch_decode(generated_ids.tolist())}")

	def decode_one_tokens(model, cur_token, input_pos, cache_position, past_key_values):
	logits = model(cur_token, position_ids=input_pos, past_key_values=past_key_values, cache_position=cache_position, return_dict=False, use_cache = True)[0]
	new_token = sample(logits)
	return new_token

	past_key_values = StaticCache(model.config, BS, max_cache_length, device=device, dtype=torch.bfloat16)
	with torch.no_grad():
	### PREFILL ##################
	cache_position = torch.arange(seq_length , device=device)
	generated_ids = torch.zeros(BS, seq_length+NUM_ITER+PREFILL_NUM_ITER+1, dtype = torch.int, device=device)
	generated_ids[:, cache_position] = input_ids.to(device).to(torch.int)

	logits = model(input_ids, past_key_values=past_key_values, cache_position=cache_position, return_dict=False, use_cache = True)[0]
	next_token = sample(logits)
	generated_ids[:, seq_length] = next_token[:,0]

	decode_one_tokens = torch.compile(decode_one_tokens, mode="reduce-overhead", fullgraph=True)

	### DRY RUN ##################
	cache_position = torch.tensor([seq_length] , device=device)
	for i in range(1, PREFILL_NUM_ITER):
	next_token = decode_one_tokens(model, next_token.clone(), None, cache_position, past_key_values)
	generated_ids.index_copy_(1, cache_position, next_token)
	cache_position+=1

	### RUN ##################
	start = time.perf_counter()
	for i in range(PREFILL_NUM_ITER, NUM_ITER+PREFILL_NUM_ITER):
	next_token = decode_one_tokens(model, next_token.clone(), None, cache_position, past_key_values)
	generated_ids.index_copy_(1, cache_position, next_token)
	cache_position+=1

	torch.cuda.synchronize()
	total_time = time.perf_counter() - start
	forward_average_time = total_time / NUM_ITER
	print(f"{BS},\t{seq_length},\t{max_cache_length}: compiled static cache\t - Took : {total_time:.4f} sec total, {forward_average_time*1000:.3f} ms per forward, {(NUM_ITER/total_time):.2f} tokens/s")
	# print(f"{tokenizer.batch_decode(generated_ids[:,seq_length:].tolist())}")
	#####################################################################################