| from pathlib import Path |
| import sys,os,torch |
| import torch.distributed as dist |
| from torch.nn.parallel import DistributedDataParallel as DDP |
| ROOT=Path(__file__).resolve().parents[1];sys.path.insert(0,str(ROOT)) |
| from model.wrf_ml import * |
| c=load_config(ROOT);rank=int(os.getenv("RANK",0));world=int(os.getenv("WORLD_SIZE",1));distributed=world>1 |
| if distributed:dist.init_process_group("gloo") |
| torch.manual_seed(c["seed"]);torch.set_num_threads(2);base=RadiationBiLSTM(**c["model"]);m=DDP(base) if distributed else base;opt=torch.optim.Adam(m.parameters(),lr=c["train"]["learning_rate"]);losses=[] |
| for i in range(rank,c["data"]["samples"],world):x,p,b=synthetic_columns(i);ph,bh=m(x);loss=((ph-p)/torch.tensor([500,300,5,5])) .pow(2).mean()+((bh-b)/500).pow(2).mean();opt.zero_grad();loss.backward();opt.step();losses.append(float(loss)) |
| v=torch.tensor([sum(losses),len(losses)],dtype=torch.float64) |
| if distributed:dist.all_reduce(v) |
| path=ROOT/c["paths"]["checkpoint"] |
| if rank==0:path.parent.mkdir(parents=True,exist_ok=True);torch.save({"model":base.state_dict(),"model_config":c["model"]},path);write_json(ROOT/c["paths"]["training_metrics"],{"normalized_mse":float(v[0]/v[1]),"world_size":world});print(path) |
| if distributed:dist.destroy_process_group() |
|
|
