Skip to main content

Write a function with a for loop with tensorflow.compat.v1 to implement a integral

Hi! I'm working with deepxde using tensorflow.compat.v1 as a backend and I need to approximate an integral using tensorflow but I haven't able to figure it out. As context the integral is

integral

The difficulty here is that the variable t is the upper limit of the integral but also is in the function q. So for each t I need to create a tensor and iterate it.

This is my code so far

def tf_integral(t_tensor, a):
    dt = tf.constant(0.1)
    ta = tf.TensorArray(tf.float32, size=0, dynamic_size=True)
    ta = ta.unstack([0., 1.])
    for i in range(n_days):
        t = t_tensor[i]
        # breakpoint()
        eta_tensor = tf.range(dt, t + dt, dt)
        value = tf.reduce_sum(tf_q(t - eta_tensor, b, c) * tf.exp(eta_tensor, d)) * dx
        ta = ta.write(i, value)
    return ta.stack()

where t_tensor is a tensor of shape [None, 1] (a column), n_days is a integer since I wasn't able to get something like len(t_tensor). I'm getting the following error.

ValueError: in user code:

    File "/mnt/d/Documents/git/insight/notebooks/china project/china-v-fitted.py", line 130, in tf_integral_dv_dphi  *
        nu_tensor = tf.range(dx, t + dx, dx)

    ValueError: Shape must be rank 0 but is rank 1
         for 'limit' for ' = Range[Tidx=DT_FLOAT](Const, add, Const)' with input shapes: [], [1], [].

I'm not sure if I can use @tf.function decorator with tensorflow.compat.v1 backend since it is my first time working with tensorflow.

I've been looking for alternatives and it seems I could use tf.compat.v1.while_loop and/or tf.compat.v1.wrap_function but honestly I'm not sure.

As a reference, the code with scipy works well

def integral(t, a):
    dt = 0.1
    eta_array = np.arange(dt, t + dt, dt)
    value = np.sum(q(t - eta_array) * np.exp(a * eta_array)) * dx
    return value

Thanks in advance!



source https://stackoverflow.com/questions/75961673/write-a-function-with-a-for-loop-with-tensorflow-compat-v1-to-implement-a-integr

Comments

Popular posts from this blog

Where and how is this Laravel kernel constructor called? [closed]

Where and how is this Laravel kernel constructor called? public fucntion __construct(Application $app, $Router $roouter) { } I have read the documentation and some online tutorial but I can find any clear explanation. I am learning Laravel and I am wondering where does this kernel constructor receives its arguments from. "POSTMOTERM" CLARIFICATION: Here is more clarity.I have checked the boostrap/app.php and it is only used for boostrapping the interfaces into the container class. What is not clear to me is where and how the Kernel class is instatiated and the arguments passed to the object calling the constructor.Something similar to; obj = new kernel(arg1,arg2) or, is the framework using some magic functions somewhere? Special gratitude to those who burn their eyeballs and brain cells on this trivia before it goes into a full blown menopause alias "MARKED AS DUPLICATE". To some of the itchy-finger keyboard warriors, a.k.a The mods,because I believe in th...

Why is my reports service not connecting?

I am trying to pull some data from a Postgres database using Node.js and node-postures but I can't figure out why my service isn't connecting. my routes/index.js file: const express = require('express'); const router = express.Router(); const ordersCountController = require('../controllers/ordersCountController'); const ordersController = require('../controllers/ordersController'); const weeklyReportsController = require('../controllers/weeklyReportsController'); router.get('/orders_count', ordersCountController); router.get('/orders', ordersController); router.get('/weekly_reports', weeklyReportsController); module.exports = router; My controllers/weeklyReportsController.js file: const weeklyReportsService = require('../services/weeklyReportsService'); const weeklyReportsController = async (req, res) => { try { const data = await weeklyReportsService; res.json({data}) console...

Confusion between commands.Bot and discord.Client | Which one should I use?

Whenever you look at YouTube tutorials or code from this website there is a real variation. Some developers use client = discord.Client(intents=intents) while the others use bot = commands.Bot(command_prefix="something", intents=intents) . Now I know slightly about the difference but I get errors from different places from my code when I use either of them and its confusing. Especially since there has a few changes over the years in discord.py it is hard to find the real difference. I tried sticking to discord.Client then I found that there are more features in commands.Bot . Then I found errors when using commands.Bot . An example of this is: When I try to use commands.Bot client = commands.Bot(command_prefix=">",intents=intents) async def load(): for filename in os.listdir("./Cogs"): if filename.endswith(".py"): client.load_extension(f"Cogs.{filename[:-3]}") The above doesnt giveany response from my Cogs ...