all the input arrays must have same number of dimensions, but the array at index 0 has 2 dimension(s) and the array at index 1 has 1 dimension(s)

[athtareq]

Hello, I'm trying to solve

Here's the code:

from __future__ import absolute_import
from __future__ import division
from __future__ import print_function

import numpy as np
import tensorflow as tf

import deepxde as dde


def main():
    def pde (x , y ) :
        dy_x = tf . gradients (y , x ) [0]
        dy_x , dy_y = dy_x [: , 0:1] , dy_x [: , 1:]
        dy_xx = tf . gradients ( dy_x , x ) [0][: , 0:1]
        dy_yy = tf . gradients ( dy_y , x ) [0][: , 1:]
        return dy_xx + dy_yy - x[: , 0:1]*tf.sin(5.0*np.pi*x[: , 1:])+tf.exp(-((x[: , 0:1]-0.5)**2+(x[: , 1:]-0.5)**2)/0.02) 

    def func (x) :
        return np.zeros([len(x),1])

    def boundary(x, on_boundary):
        return on_boundary

    def boundary_l(x, on_boundary):
        return on_boundary and np.isclose(x[1], 0)

    def boundary_r(x, on_boundary):
        return on_boundary and np.isclose(x[1], 1)

    def boundary_t(x, on_boundary):
        return on_boundary and np.isclose(x[0], 1)

    def boundary_b(x, on_boundary):
        return on_boundary and np.isclose(x[0], 0)

    geom = dde.geometry.Rectangle([0,0],[1,1])
    geomtime = dde.geometry.Interval(0, 1)
    bc = dde.DirichletBC(geom, func, boundary)
    bc1 = dde.DirichletBC(geom, func, boundary_t)
    bc2 = dde.DirichletBC(geom, func, boundary_b)
    bc3 = dde.PeriodicBC(geomtime, 1, boundary_l, derivative_order=0, component=1)
    #bc4 = dde.PeriodicBC(geom, 0, boundary_b)


    data = dde.data.PDE(
        geom, pde, [bc1, bc2, bc3], num_domain=1200, num_boundary=120, num_test=1500
    )
    layer_size = [2] + [50] * 3 + [1]
    activation = "tanh"
    initializer = "Glorot uniform"
    net = dde.maps.FNN(layer_size, activation, initializer)

    model = dde.Model(data, net)
    model.compile("adam", lr=1e-3, loss="MSE")
    losshistory, train_state = model.train(epochs=50000)

    dde.saveplot(losshistory, train_state, issave=True, isplot=True)


if __name__ == "__main__":
    main()

I get the error in the title. Can anyone please tell me what's wrong ? Thanks !

[athtareq]

The error vanishes when I change this portion:

    bc3 = dde.PeriodicBC(geom, 0, boundary_t)
    bc4 = dde.PeriodicBC(geom, 0, boundary_b)

    data = dde.data.PDE(
        geom, pde, [bc1, bc2, bc3,bc4], num_domain=120, num_boundary=120, num_test=1500
    )
    metrics_=None
    net = dde.maps.FNN([2] + [50] * 4 + [1], "tanh", "Glorot uniform")
    model = dde.Model(data, net)

    model.compile("adam", lr=0.001, metrics=metrics_)
    model.train(epochs=20000)
    model.compile("L-BFGS-B", metrics=metrics_)
    losshistory, train_state = model.train()
    dde.saveplot(losshistory, train_state, issave=True, isplot=True)

if __name__ == "__main__":
    main()

Any ideas where it might have come from ? I get this error quite frequently

[forxltk]

geotime is an interval with only x, but bc3 have x[1]. And the output of your network is 1, while in bc3, component=1. I think you should have a check.

[lululxvi]

You also mixed the usage of geom and geomtime.