Light

kc-ml2 / meent Goto Github PK

Electromagnetic simulation (RCWA) & optimization package in Python

Home Page: https://meent.readthedocs.io/

License: MIT License

Python 100.00%

optimization jax pytorch optics rcwa fmm fourier-modal-method rigorous-coupled-wave-analysis electromagnetics simulation

meent's People

Contributors

Stargazers

Watchers

Forkers

chocopi2718 dalerxli jeongminshin kvnoct sprbull jackcheng1 eminsight

meent's Issues

modify 'setup.py' file

modify not to contain other folders than 'meent'

optimizer structure

meent/meent/on_torch/optimizer.py

Lines 52 to 80 in 598d828

    
           class Optimizer: 
        
               def __init__(self, optimizer, solver, target, lr = 0.001): 
        
                   self.solver = solver 
        
                   self.target = target 
        
                   if not isinstance(target, torch.Tensor): 
        
                       self.target = torch.Tensor(self.target) 
        
                   self.target.requires_grad = True 
        
                   if isinstance(optimizer,str): 
        
                       optimizer = getattr(torch.optim, optimizer) 
        
                   self.optimizer = optimizer([self.target], lr = lr) 
        
               def optimize(self, iterations = 1000): 
        
                   for iteration in range(iterations): 
        
                       E_conv_all = to_conv_mat(self.solver.ucell, fourier_order) 
        
                       o_E_conv_all = to_conv_mat(1 / self.solver.ucell, fourier_order) 
        
                       de_ri, de_ti, _, _, _ = self.solver.solve(self.solver.wavelength, E_conv_all, o_E_conv_all) 
        
                       self.optimizer.zero_grad() 
        
                       loss = self.loss(de_ti) 
        
                       loss.backward() 
        
                       self.optimizer.step() 
        
                       print(loss) 
        
               def loss(self, de_ti): 
        
                   return -de_ti[3, 2]

@yonghakim

Check and enforce Fourier_order type as integer

Even though Fourier_order was set to integer (e.g. 40) by user, it's changed to float (e.g. 40.0) and raises TypeError.

force wavelength to be float

Recommend Projects

React

A declarative, efficient, and flexible JavaScript library for building user interfaces.
Vue.js

🖖 Vue.js is a progressive, incrementally-adoptable JavaScript framework for building UI on the web.
Typescript

TypeScript is a superset of JavaScript that compiles to clean JavaScript output.
TensorFlow

An Open Source Machine Learning Framework for Everyone
Django

The Web framework for perfectionists with deadlines.
Laravel

A PHP framework for web artisans
D3

Bring data to life with SVG, Canvas and HTML. 📊📈🎉

Recommend Topics

javascript

JavaScript (JS) is a lightweight interpreted programming language with first-class functions.
web

Some thing interesting about web. New door for the world.
server

A server is a program made to process requests and deliver data to clients.
Machine learning

Machine learning is a way of modeling and interpreting data that allows a piece of software to respond intelligently.
Visualization

Some thing interesting about visualization, use data art
Game

Some thing interesting about game, make everyone happy.

Recommend Org

Facebook

We are working to build community through open source technology. NB: members must have two-factor auth.
Microsoft

Open source projects and samples from Microsoft.
Google

Google ❤️ Open Source for everyone.
Alibaba

Alibaba Open Source for everyone
D3

Data-Driven Documents codes.
Tencent

China tencent open source team.

	class Optimizer:
	def __init__(self, optimizer, solver, target, lr = 0.001):
	self.solver = solver

	self.target = target
	if not isinstance(target, torch.Tensor):
	self.target = torch.Tensor(self.target)
	self.target.requires_grad = True

	if isinstance(optimizer,str):
	optimizer = getattr(torch.optim, optimizer)
	self.optimizer = optimizer([self.target], lr = lr)

	def optimize(self, iterations = 1000):
	for iteration in range(iterations):
	E_conv_all = to_conv_mat(self.solver.ucell, fourier_order)
	o_E_conv_all = to_conv_mat(1 / self.solver.ucell, fourier_order)

	de_ri, de_ti, _, _, _ = self.solver.solve(self.solver.wavelength, E_conv_all, o_E_conv_all)

	self.optimizer.zero_grad()
	loss = self.loss(de_ti)
	loss.backward()
	self.optimizer.step()

	print(loss)

	def loss(self, de_ti):
	return -de_ti[3, 2]