X-Git-Url: https://www.fleuret.org/cgi-bin/gitweb/gitweb.cgi?a=blobdiff_plain;f=vignette_set.py;h=0b6de7e17d7db8fc4897f844c889acf549cda562;hb=c5660d2adc37f116088867ae71f66a89b34c94f1;hp=695fed3442b862faee2984adb406a946163aad5e;hpb=0b891219d91e981f96e5321bcf0db6c3beea0017;p=pysvrt.git diff --git a/vignette_set.py b/vignette_set.py index 695fed3..0b6de7e 100755 --- a/vignette_set.py +++ b/vignette_set.py @@ -22,6 +22,7 @@ import torch from math import sqrt +from torch import multiprocessing from torch import Tensor from torch.autograd import Variable @@ -30,38 +31,54 @@ import svrt ###################################################################### +def generate_one_batch(s): + problem_number, batch_size, random_seed = s + svrt.seed(random_seed) + target = torch.LongTensor(batch_size).bernoulli_(0.5) + input = svrt.generate_vignettes(problem_number, target) + input = input.float().view(input.size(0), 1, input.size(1), input.size(2)) + return [ input, target ] + class VignetteSet: + def __init__(self, problem_number, nb_batches, batch_size, cuda = False): self.cuda = cuda self.batch_size = batch_size self.problem_number = problem_number self.nb_batches = nb_batches self.nb_samples = self.nb_batches * self.batch_size - self.targets = [] - self.inputs = [] + + seeds = torch.LongTensor(self.nb_batches).random_() + mp_args = [] + for b in range(0, self.nb_batches): + mp_args.append( [ problem_number, batch_size, seeds[b] ]) + + self.data = [] + for b in range(0, self.nb_batches): + self.data.append(generate_one_batch(mp_args[b])) + + # Weird thing going on with the multi-processing, waiting for more info + + # pool = multiprocessing.Pool(multiprocessing.cpu_count()) + # self.data = pool.map(generate_one_batch, mp_args) acc = 0.0 acc_sq = 0.0 - for b in range(0, self.nb_batches): - target = torch.LongTensor(self.batch_size).bernoulli_(0.5) - input = svrt.generate_vignettes(problem_number, target) - input = input.float().view(input.size(0), 1, input.size(1), input.size(2)) - if self.cuda: - input = input.cuda() - target = target.cuda() + input = self.data[b][0] acc += input.sum() / input.numel() acc_sq += input.pow(2).sum() / input.numel() - self.targets.append(target) - self.inputs.append(input) mean = acc / self.nb_batches std = sqrt(acc_sq / self.nb_batches - mean * mean) for b in range(0, self.nb_batches): - self.inputs[b].sub_(mean).div_(std) + self.data[b][0].sub_(mean).div_(std) + if cuda: + self.data[b][0] = self.data[b][0].cuda() + self.data[b][1] = self.data[b][1].cuda() def get_batch(self, b): - return self.inputs[b], self.targets[b] + return self.data[b] ######################################################################