xautodl/exps/NAS-Bench-201/statistics-v2.py

#####################################################
# Copyright (c) Xuanyi Dong [GitHub D-X-Y], 2019.08 #
#####################################################
import os, sys, time, argparse, collections
import numpy as np
import torch
from pathlib import Path
from collections import defaultdict, OrderedDict
from typing import Dict, Any, Text, List

lib_dir = (Path(__file__).parent / ".." / ".." / "lib").resolve()
if str(lib_dir) not in sys.path:
    sys.path.insert(0, str(lib_dir))
from log_utils import AverageMeter, time_string, convert_secs2time
from config_utils import dict2config

# NAS-Bench-201 related module or function
from models import CellStructure, get_cell_based_tiny_net
from nas_201_api import NASBench201API, ArchResults, ResultsCount
from procedures import bench_pure_evaluate as pure_evaluate, get_nas_bench_loaders

api = NASBench201API(
    "{:}/.torch/NAS-Bench-201-v1_0-e61699.pth".format(os.environ["HOME"])
)


def create_result_count(
    used_seed: int,
    dataset: Text,
    arch_config: Dict[Text, Any],
    results: Dict[Text, Any],
    dataloader_dict: Dict[Text, Any],
) -> ResultsCount:
    xresult = ResultsCount(
        dataset,
        results["net_state_dict"],
        results["train_acc1es"],
        results["train_losses"],
        results["param"],
        results["flop"],
        arch_config,
        used_seed,
        results["total_epoch"],
        None,
    )
    net_config = dict2config(
        {
            "name": "infer.tiny",
            "C": arch_config["channel"],
            "N": arch_config["num_cells"],
            "genotype": CellStructure.str2structure(arch_config["arch_str"]),
            "num_classes": arch_config["class_num"],
        },
        None,
    )
    network = get_cell_based_tiny_net(net_config)
    network.load_state_dict(xresult.get_net_param())
    if "train_times" in results:  # new version
        xresult.update_train_info(
            results["train_acc1es"],
            results["train_acc5es"],
            results["train_losses"],
            results["train_times"],
        )
        xresult.update_eval(
            results["valid_acc1es"], results["valid_losses"], results["valid_times"]
        )
    else:
        if dataset == "cifar10-valid":
            xresult.update_OLD_eval(
                "x-valid", results["valid_acc1es"], results["valid_losses"]
            )
            loss, top1, top5, latencies = pure_evaluate(
                dataloader_dict["{:}@{:}".format("cifar10", "test")], network.cuda()
            )
            xresult.update_OLD_eval(
                "ori-test",
                {results["total_epoch"] - 1: top1},
                {results["total_epoch"] - 1: loss},
            )
            xresult.update_latency(latencies)
        elif dataset == "cifar10":
            xresult.update_OLD_eval(
                "ori-test", results["valid_acc1es"], results["valid_losses"]
            )
            loss, top1, top5, latencies = pure_evaluate(
                dataloader_dict["{:}@{:}".format(dataset, "test")], network.cuda()
            )
            xresult.update_latency(latencies)
        elif dataset == "cifar100" or dataset == "ImageNet16-120":
            xresult.update_OLD_eval(
                "ori-test", results["valid_acc1es"], results["valid_losses"]
            )
            loss, top1, top5, latencies = pure_evaluate(
                dataloader_dict["{:}@{:}".format(dataset, "valid")], network.cuda()
            )
            xresult.update_OLD_eval(
                "x-valid",
                {results["total_epoch"] - 1: top1},
                {results["total_epoch"] - 1: loss},
            )
            loss, top1, top5, latencies = pure_evaluate(
                dataloader_dict["{:}@{:}".format(dataset, "test")], network.cuda()
            )
            xresult.update_OLD_eval(
                "x-test",
                {results["total_epoch"] - 1: top1},
                {results["total_epoch"] - 1: loss},
            )
            xresult.update_latency(latencies)
        else:
            raise ValueError("invalid dataset name : {:}".format(dataset))
    return xresult


def account_one_arch(
    arch_index: int,
    arch_str: Text,
    checkpoints: List[Text],
    datasets: List[Text],
    dataloader_dict: Dict[Text, Any],
) -> ArchResults:
    information = ArchResults(arch_index, arch_str)

    for checkpoint_path in checkpoints:
        checkpoint = torch.load(checkpoint_path, map_location="cpu")
        used_seed = checkpoint_path.name.split("-")[-1].split(".")[0]
        ok_dataset = 0
        for dataset in datasets:
            if dataset not in checkpoint:
                print(
                    "Can not find {:} in arch-{:} from {:}".format(
                        dataset, arch_index, checkpoint_path
                    )
                )
                continue
            else:
                ok_dataset += 1
            results = checkpoint[dataset]
            assert results[
                "finish-train"
            ], "This {:} arch seed={:} does not finish train on {:} ::: {:}".format(
                arch_index, used_seed, dataset, checkpoint_path
            )
            arch_config = {
                "channel": results["channel"],
                "num_cells": results["num_cells"],
                "arch_str": arch_str,
                "class_num": results["config"]["class_num"],
            }

            xresult = create_result_count(
                used_seed, dataset, arch_config, results, dataloader_dict
            )
            information.update(dataset, int(used_seed), xresult)
        if ok_dataset == 0:
            raise ValueError("{:} does not find any data".format(checkpoint_path))
    return information


def correct_time_related_info(
    arch_index: int, arch_info_full: ArchResults, arch_info_less: ArchResults
):
    # calibrate the latency based on NAS-Bench-201-v1_0-e61699.pth
    cifar010_latency = (
        api.get_latency(arch_index, "cifar10-valid", hp="200")
        + api.get_latency(arch_index, "cifar10", hp="200")
    ) / 2
    arch_info_full.reset_latency("cifar10-valid", None, cifar010_latency)
    arch_info_full.reset_latency("cifar10", None, cifar010_latency)
    arch_info_less.reset_latency("cifar10-valid", None, cifar010_latency)
    arch_info_less.reset_latency("cifar10", None, cifar010_latency)

    cifar100_latency = api.get_latency(arch_index, "cifar100", hp="200")
    arch_info_full.reset_latency("cifar100", None, cifar100_latency)
    arch_info_less.reset_latency("cifar100", None, cifar100_latency)

    image_latency = api.get_latency(arch_index, "ImageNet16-120", hp="200")
    arch_info_full.reset_latency("ImageNet16-120", None, image_latency)
    arch_info_less.reset_latency("ImageNet16-120", None, image_latency)

    train_per_epoch_time = list(
        arch_info_less.query("cifar10-valid", 777).train_times.values()
    )
    train_per_epoch_time = sum(train_per_epoch_time) / len(train_per_epoch_time)
    eval_ori_test_time, eval_x_valid_time = [], []
    for key, value in arch_info_less.query("cifar10-valid", 777).eval_times.items():
        if key.startswith("ori-test@"):
            eval_ori_test_time.append(value)
        elif key.startswith("x-valid@"):
            eval_x_valid_time.append(value)
        else:
            raise ValueError("-- {:} --".format(key))
    eval_ori_test_time, eval_x_valid_time = float(np.mean(eval_ori_test_time)), float(
        np.mean(eval_x_valid_time)
    )
    nums = {
        "ImageNet16-120-train": 151700,
        "ImageNet16-120-valid": 3000,
        "ImageNet16-120-test": 6000,
        "cifar10-valid-train": 25000,
        "cifar10-valid-valid": 25000,
        "cifar10-train": 50000,
        "cifar10-test": 10000,
        "cifar100-train": 50000,
        "cifar100-test": 10000,
        "cifar100-valid": 5000,
    }
    eval_per_sample = (eval_ori_test_time + eval_x_valid_time) / (
        nums["cifar10-valid-valid"] + nums["cifar10-test"]
    )
    for arch_info in [arch_info_less, arch_info_full]:
        arch_info.reset_pseudo_train_times(
            "cifar10-valid",
            None,
            train_per_epoch_time
            / nums["cifar10-valid-train"]
            * nums["cifar10-valid-train"],
        )
        arch_info.reset_pseudo_train_times(
            "cifar10",
            None,
            train_per_epoch_time / nums["cifar10-valid-train"] * nums["cifar10-train"],
        )
        arch_info.reset_pseudo_train_times(
            "cifar100",
            None,
            train_per_epoch_time / nums["cifar10-valid-train"] * nums["cifar100-train"],
        )
        arch_info.reset_pseudo_train_times(
            "ImageNet16-120",
            None,
            train_per_epoch_time
            / nums["cifar10-valid-train"]
            * nums["ImageNet16-120-train"],
        )
        arch_info.reset_pseudo_eval_times(
            "cifar10-valid",
            None,
            "x-valid",
            eval_per_sample * nums["cifar10-valid-valid"],
        )
        arch_info.reset_pseudo_eval_times(
            "cifar10-valid", None, "ori-test", eval_per_sample * nums["cifar10-test"]
        )
        arch_info.reset_pseudo_eval_times(
            "cifar10", None, "ori-test", eval_per_sample * nums["cifar10-test"]
        )
        arch_info.reset_pseudo_eval_times(
            "cifar100", None, "x-valid", eval_per_sample * nums["cifar100-valid"]
        )
        arch_info.reset_pseudo_eval_times(
            "cifar100", None, "x-test", eval_per_sample * nums["cifar100-valid"]
        )
        arch_info.reset_pseudo_eval_times(
            "cifar100", None, "ori-test", eval_per_sample * nums["cifar100-test"]
        )
        arch_info.reset_pseudo_eval_times(
            "ImageNet16-120",
            None,
            "x-valid",
            eval_per_sample * nums["ImageNet16-120-valid"],
        )
        arch_info.reset_pseudo_eval_times(
            "ImageNet16-120",
            None,
            "x-test",
            eval_per_sample * nums["ImageNet16-120-valid"],
        )
        arch_info.reset_pseudo_eval_times(
            "ImageNet16-120",
            None,
            "ori-test",
            eval_per_sample * nums["ImageNet16-120-test"],
        )
    # arch_info_full.debug_test()
    # arch_info_less.debug_test()
    return arch_info_full, arch_info_less


def simplify(save_dir, meta_file, basestr, target_dir):
    meta_infos = torch.load(meta_file, map_location="cpu")
    meta_archs = meta_infos["archs"]  # a list of architecture strings
    meta_num_archs = meta_infos["total"]
    assert meta_num_archs == len(
        meta_archs
    ), "invalid number of archs : {:} vs {:}".format(meta_num_archs, len(meta_archs))

    sub_model_dirs = sorted(list(save_dir.glob("*-*-{:}".format(basestr))))
    print(
        "{:} find {:} directories used to save checkpoints".format(
            time_string(), len(sub_model_dirs)
        )
    )

    subdir2archs, num_evaluated_arch = collections.OrderedDict(), 0
    num_seeds = defaultdict(lambda: 0)
    for index, sub_dir in enumerate(sub_model_dirs):
        xcheckpoints = list(sub_dir.glob("arch-*-seed-*.pth"))
        arch_indexes = set()
        for checkpoint in xcheckpoints:
            temp_names = checkpoint.name.split("-")
            assert (
                len(temp_names) == 4
                and temp_names[0] == "arch"
                and temp_names[2] == "seed"
            ), "invalid checkpoint name : {:}".format(checkpoint.name)
            arch_indexes.add(temp_names[1])
        subdir2archs[sub_dir] = sorted(list(arch_indexes))
        num_evaluated_arch += len(arch_indexes)
        # count number of seeds for each architecture
        for arch_index in arch_indexes:
            num_seeds[
                len(list(sub_dir.glob("arch-{:}-seed-*.pth".format(arch_index))))
            ] += 1
    print(
        "{:} There are {:5d} architectures that have been evaluated ({:} in total).".format(
            time_string(), num_evaluated_arch, meta_num_archs
        )
    )
    for key in sorted(list(num_seeds.keys())):
        print(
            "{:} There are {:5d} architectures that are evaluated {:} times.".format(
                time_string(), num_seeds[key], key
            )
        )

    dataloader_dict = get_nas_bench_loaders(6)
    to_save_simply = save_dir / "simplifies"
    to_save_allarc = save_dir / "simplifies" / "architectures"
    if not to_save_simply.exists():
        to_save_simply.mkdir(parents=True, exist_ok=True)
    if not to_save_allarc.exists():
        to_save_allarc.mkdir(parents=True, exist_ok=True)

    assert (save_dir / target_dir) in subdir2archs, "can not find {:}".format(
        target_dir
    )
    arch2infos, datasets = {}, (
        "cifar10-valid",
        "cifar10",
        "cifar100",
        "ImageNet16-120",
    )
    evaluated_indexes = set()
    target_full_dir = save_dir / target_dir
    target_less_dir = save_dir / "{:}-LESS".format(target_dir)
    arch_indexes = subdir2archs[target_full_dir]
    num_seeds = defaultdict(lambda: 0)
    end_time = time.time()
    arch_time = AverageMeter()
    for idx, arch_index in enumerate(arch_indexes):
        checkpoints = list(
            target_full_dir.glob("arch-{:}-seed-*.pth".format(arch_index))
        )
        ckps_less = list(target_less_dir.glob("arch-{:}-seed-*.pth".format(arch_index)))
        # create the arch info for each architecture
        try:
            arch_info_full = account_one_arch(
                arch_index,
                meta_archs[int(arch_index)],
                checkpoints,
                datasets,
                dataloader_dict,
            )
            arch_info_less = account_one_arch(
                arch_index,
                meta_archs[int(arch_index)],
                ckps_less,
                datasets,
                dataloader_dict,
            )
            num_seeds[len(checkpoints)] += 1
        except:
            print("Loading {:} failed, : {:}".format(arch_index, checkpoints))
            continue
        assert (
            int(arch_index) not in evaluated_indexes
        ), "conflict arch-index : {:}".format(arch_index)
        assert (
            0 <= int(arch_index) < len(meta_archs)
        ), "invalid arch-index {:} (not found in meta_archs)".format(arch_index)
        arch_info = {"full": arch_info_full, "less": arch_info_less}
        evaluated_indexes.add(int(arch_index))
        arch2infos[int(arch_index)] = arch_info
        # to correct the latency and training_time info.
        arch_info_full, arch_info_less = correct_time_related_info(
            int(arch_index), arch_info_full, arch_info_less
        )
        to_save_data = OrderedDict(
            full=arch_info_full.state_dict(), less=arch_info_less.state_dict()
        )
        torch.save(to_save_data, to_save_allarc / "{:}-FULL.pth".format(arch_index))
        arch_info["full"].clear_params()
        arch_info["less"].clear_params()
        torch.save(to_save_data, to_save_allarc / "{:}-SIMPLE.pth".format(arch_index))
        # measure elapsed time
        arch_time.update(time.time() - end_time)
        end_time = time.time()
        need_time = "{:}".format(
            convert_secs2time(arch_time.avg * (len(arch_indexes) - idx - 1), True)
        )
        print(
            "{:} {:} [{:03d}/{:03d}] : {:} still need {:}".format(
                time_string(), target_dir, idx, len(arch_indexes), arch_index, need_time
            )
        )
    # measure time
    xstrs = [
        "{:}:{:03d}".format(key, num_seeds[key])
        for key in sorted(list(num_seeds.keys()))
    ]
    print("{:} {:} done : {:}".format(time_string(), target_dir, xstrs))
    final_infos = {
        "meta_archs": meta_archs,
        "total_archs": meta_num_archs,
        "basestr": basestr,
        "arch2infos": arch2infos,
        "evaluated_indexes": evaluated_indexes,
    }
    save_file_name = to_save_simply / "{:}.pth".format(target_dir)
    torch.save(final_infos, save_file_name)
    print(
        "Save {:} / {:} architecture results into {:}.".format(
            len(evaluated_indexes), meta_num_archs, save_file_name
        )
    )


def merge_all(save_dir, meta_file, basestr):
    meta_infos = torch.load(meta_file, map_location="cpu")
    meta_archs = meta_infos["archs"]
    meta_num_archs = meta_infos["total"]
    assert meta_num_archs == len(
        meta_archs
    ), "invalid number of archs : {:} vs {:}".format(meta_num_archs, len(meta_archs))

    sub_model_dirs = sorted(list(save_dir.glob("*-*-{:}".format(basestr))))
    print(
        "{:} find {:} directories used to save checkpoints".format(
            time_string(), len(sub_model_dirs)
        )
    )
    for index, sub_dir in enumerate(sub_model_dirs):
        arch_info_files = sorted(list(sub_dir.glob("arch-*-seed-*.pth")))
        print(
            "The {:02d}/{:02d}-th directory : {:} : {:} runs.".format(
                index, len(sub_model_dirs), sub_dir, len(arch_info_files)
            )
        )

    arch2infos, evaluated_indexes = dict(), set()
    for IDX, sub_dir in enumerate(sub_model_dirs):
        ckp_path = sub_dir.parent / "simplifies" / "{:}.pth".format(sub_dir.name)
        if ckp_path.exists():
            sub_ckps = torch.load(ckp_path, map_location="cpu")
            assert (
                sub_ckps["total_archs"] == meta_num_archs
                and sub_ckps["basestr"] == basestr
            )
            xarch2infos = sub_ckps["arch2infos"]
            xevalindexs = sub_ckps["evaluated_indexes"]
            for eval_index in xevalindexs:
                assert (
                    eval_index not in evaluated_indexes and eval_index not in arch2infos
                )
                # arch2infos[eval_index] = xarch2infos[eval_index].state_dict()
                arch2infos[eval_index] = {
                    "full": xarch2infos[eval_index]["full"].state_dict(),
                    "less": xarch2infos[eval_index]["less"].state_dict(),
                }
                evaluated_indexes.add(eval_index)
            print(
                "{:} [{:03d}/{:03d}] merge data from {:} with {:} models.".format(
                    time_string(), IDX, len(sub_model_dirs), ckp_path, len(xevalindexs)
                )
            )
        else:
            raise ValueError("Can not find {:}".format(ckp_path))
            # print ('{:} [{:03d}/{:03d}] can not find {:}, skip.'.format(time_string(), IDX, len(subdir2archs), ckp_path))

    evaluated_indexes = sorted(list(evaluated_indexes))
    print(
        "Finally, there are {:} architectures that have been trained and evaluated.".format(
            len(evaluated_indexes)
        )
    )

    to_save_simply = save_dir / "simplifies"
    if not to_save_simply.exists():
        to_save_simply.mkdir(parents=True, exist_ok=True)
    final_infos = {
        "meta_archs": meta_archs,
        "total_archs": meta_num_archs,
        "arch2infos": arch2infos,
        "evaluated_indexes": evaluated_indexes,
    }
    save_file_name = to_save_simply / "{:}-final-infos.pth".format(basestr)
    torch.save(final_infos, save_file_name)
    print(
        "Save {:} / {:} architecture results into {:}.".format(
            len(evaluated_indexes), meta_num_archs, save_file_name
        )
    )


if __name__ == "__main__":

    parser = argparse.ArgumentParser(
        description="NAS-BENCH-201",
        formatter_class=argparse.ArgumentDefaultsHelpFormatter,
    )
    parser.add_argument(
        "--mode",
        type=str,
        choices=["cal", "merge"],
        help="The running mode for this script.",
    )
    parser.add_argument(
        "--base_save_dir",
        type=str,
        default="./output/NAS-BENCH-201-4",
        help="The base-name of folder to save checkpoints and log.",
    )
    parser.add_argument("--target_dir", type=str, help="The target directory.")
    parser.add_argument(
        "--max_node", type=int, default=4, help="The maximum node in a cell."
    )
    parser.add_argument(
        "--channel", type=int, default=16, help="The number of channels."
    )
    parser.add_argument(
        "--num_cells", type=int, default=5, help="The number of cells in one stage."
    )
    args = parser.parse_args()

    save_dir = Path(args.base_save_dir)
    meta_path = save_dir / "meta-node-{:}.pth".format(args.max_node)
    assert save_dir.exists(), "invalid save dir path : {:}".format(save_dir)
    assert meta_path.exists(), "invalid saved meta path : {:}".format(meta_path)
    print(
        "start the statistics of our nas-benchmark from {:} using {:}.".format(
            save_dir, args.target_dir
        )
    )
    basestr = "C{:}-N{:}".format(args.channel, args.num_cells)

    if args.mode == "cal":
        simplify(save_dir, meta_path, basestr, args.target_dir)
    elif args.mode == "merge":
        merge_all(save_dir, meta_path, basestr)
    else:
        raise ValueError("invalid mode : {:}".format(args.mode))